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FOREWORD 


In  recent  years  many  engineers  concerned  with  dynamic 
motion  environments  have  been  exposed,  in  one  way  or  an¬ 
other,  to  the  term  “mechanical  impedance."  Dependent  upon 
his  available  time  and  inclination  he  may  have  looked  in 
usual  texts  and  literature  and  determined  that  mechanical 
impedance  is  the  “analog"  of  electirical  impedance.  He  may 
also  have  seen  other  terms,  e.g.,  “mobility"  or  “admittance* 
which  appear  to  be  very  similar.  However,  as  Judged  by 
collective  remarks  of  shock  and  vibration  engineers,  what 
he  probably  has  not  seen  is  an  explanation  of  how  this  quan¬ 
tity  can  be  useful  to  him  in  his  day-to-day  problems.  It  was 
just  this  situation  which  prompted  the  scheduling  of  a  session 
of  arranged  and  invited  papers  for  the  30th  Symposium  on 
Shock,  Vibration,  and  Associated  Elnvironments.  The  intent 
was  to  present  the  concepts  of  mechanical  impedance  in 
“tutorial"  fashion,  to  describe  some  applications,  and  to 
encourage  discussion  of  the  subject.  The  technical  content 
of  the  session  and  the  invited  papers  were  arranged  by 
R.  O.  Belsheim  of  the  U.S.  Naval  Research  Laboratory. 

The  papers  reprinted  here  contain  the  material  pre¬ 
sented  at  the  Symposium  together  with  the  discussion  which 
followed  the  presentations  and  an  additional  paper,  not  pre¬ 
sented  at  the  meeting,  on  the  same  subject. 

January  1962 
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Section  2 

MECHANICAL  IMPEDANCE 


INTRODUCTORY  REMARKS 


R.  O.  Belsheim 

U.  S.  Naval  Research  Laboratory 


These  papers  have  been  arranged  to  present,  in  a  "tutorial"  sense,  the 
concept,  determination,  and  application  of  what  may  be  referred  to  as 
mechanical  impedance  procedures.  Some  historical  and  introductory 
material  is  presented  to  prepare  the  way  for  the  papers. 


During  the  last  decade  or  so,  a  marked 
change  in  the  speed  with  which  operational  en> 
gi..eering  events  occur  has  greatly  increased 
the  scope  and  number  of  problems  to  be  solved 
by  the  structural  dynamiclst,  test  engineer,  de¬ 
signer,  and  others.  This  increase  is  caused  by 
the  ne^  for  greater  performance  reliability, 
but  with  reduced  weight  of  structures,  and  with 
more  complex  and  indeterminate  motional  en¬ 
vironments.  To  meet  this  need  dynamic  anal¬ 
ysis,  design,  and  testing  must  be  done  more 
exactly,  i.e.,  with  less  "factor  of  ignorance" 
margin.  In  a  historical  sense,  the  electrical 
engineers  faced  a  similar  situation  some  sixty 
years  ago,  with  respect  to  circuit  analyses. 
Their  solution  was  to  develop  methods,  tech¬ 
niques,  and  theorems  which  permit  considera¬ 
tion  of  complex  circuits  by  use  of  a  quantitative 
parameter  called  impedance  [1,2];"'  this  pa¬ 
rameter  is  the  ratio  of  sinusoidal  voltage  and 
sinusoidal  current.  Because  many  circuits  are 
linear,  this  ratio  is  independent  of  amplitude 
and  is  a  function  of  sinusoidal  frequency  only. 
However,  because  voltage  and  current  are  not 
in  phase  with  each  other,  their  ratio,  z,  can  be 
expressed  as  a  complex  number  with  a  magni¬ 
tude  and  a  phase  angle,  or,  alternatively,  with 
real  and  imaginary  parts.  Over  the  years  de¬ 
velopment  and  evaluation  of  a  great  many  cir¬ 
cuit  theorems,  design  rules,  etc.,  has  almost 


eliminated  the  classical  method  of  solving  linear 
circuit  problems  by  use  of  the  pertinent  differ¬ 
ential  equations. 

Workers  in  the  field  of  acoustics,  particu¬ 
larly  in  areas  where  electromechanical  devices 
are  important,  have  also  used  similar  concepts 
for  some  30  to  40  years  [3, 4].  Because  energy 
is  often  converted  from  electrical  to  mechani¬ 
cal  forms  (or  vice  versa)  extensive  study  of 
electrical- mechanical  analogs  [5]  has  proved 
useful  and  impedance  of  mixed  units,  e.g.,  output 
force  to  input  current,  is  used. 

The  field  of  control  by  servomechanisms, 
which  has  mushroomed  during  the  last  20  to  30 
years,  also  uses  similar  parameters  to  facili¬ 
tate  analysis  of  control  systems  [6,7].  Because 
these  systems  are  specially  designed  it  has  been 
possible  to  develop  qiecialized  (and  simple) 
techniques,  since  certain  restrictions  are  im¬ 
posed  in  the  design  stages. 


MEANING  OF  MECHANICAL 
IMPEDANCE 

With  regard  to  linear  and  elastic  struc¬ 
tures  or  mechanical  systems,  the  literature 
during  the  past  25  years  contains  a  number  of 


♦Numbers  in  brackets  refer  to  references  at  end. 
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papers  which  discuss  and/or  utilize  mechanical 
impedance  concepts.  Of  the  earlier  papers 
one  should  mention  the  works  of  Carter  [8] , 
Firestone  [9],  Biot  [10],  Duncan  [11],  and 
Manley  [12],  all  of  which  were  concerned  with 
systems  subjected  to  free  or  forced  vibration. 
Much  of  the  early  work  utilized  the  analogy  be¬ 
tween  the  equations  governing  mechanical  sys¬ 
tems  and  electrical  circuits  and  did  not  really 
deal  with  impedance  quantities  as  a  unique  pa¬ 
rameter.  This  brings  up  a  misconception  which 
some  engineers  may  have,  i.e.,  that  impedance 
is  something  new  and  unique  which  has  become 
available  suddenly.  Of  course  this  is  not  the 
case;  mechanical  impedance  can  be  an  extremely 
useful  "tool”  for  solving  problems  in  structural 
dynamics  but,  as  with  any  tool  it  should  be  used 
with  proper  regard  for  its  limitations  and  in 
conjunction  with  other  methods  and  techniques. 
Its  unique  quality  is  that  it  permits  a  convenient 
yet  precise  description  of  the  dynamic  charac¬ 
teristics  of  a  structure,  either  at  a  point  or  be¬ 
tween  two  points;  this  descr4>tion  is  equally 
meaningful  regardless  of  the  complexity  of  the 
structure,  and  allows  "equivalent  black  box" 
manipulations. 

Although  Impedance  is  a  parameter  derived 
from  ratios  of  sinusoidal  quantities,  a  structure 
properly  described  by  its  impedance  may  be 
analyzed  for  response  to  sinusoidal,  random, 
or  transient  excitations  [11, 13-17]. 

The  limitations  on  use  of  impedance  may 
be  grouped  to  Include:  (a)  the  assumption  that 
the  system  treated  is  linear  and  elastic,  (b)  the 
requirement  that  the  desired  impedance  quan¬ 
tity  be  determined  (analytically  or  experimen¬ 
tally)  in  a  reasonable  time,  yet  with  sufficient 
precision  [18],  and  (c)  the  necessity  for  exten¬ 
sive  computations  which  may  require  electronic 
computation  [19].  None  of  these  limitations  are 
really  precise  so,  as  with  electrical  network 
theory,  the  supporting  theorems,  rules,  and 
procedures  must  be  developed  and  proved  over 
a  period  of  time.  Because  the  mechanical- 
electrical  analogy  is  applicable,  a  considerable 
fund  of  such  theorems  is  already  available,  al¬ 
though  each  should  be  reconsidered  in  terms  of 
the  aforementioned  limitations  on  mechanical 
impedance.  Also,  even  though  it  may  be 


convenient  to  solve  mechanical  problems  by 
using  an  analog,  it  is  desirable  to  present 
mechanical  impedance  in  its  own  language,  so 
that  the  required  limitations  are  correctly 
understood  and  applied. 


DEHNITION  OF  MECHANICAL 
IMPEDANCE 

Mechanical  impedance  has  been  defined 
loosely  as  the  complex  ratio  of  sinusoidal  force 
to  sinusoidal  velocity.  More  precisely,  if  the 
sinusoidal  force  F  and  the  resulting  sinusoidal 
velocity  v  are  given  by 


F  =  F,j  cos  (ait  -f  0) 

(1) 

V  =  Vjj  cos  cat  , 

(2) 

the  impedance  Z  by  definition  is 

(3) 

where  6 ,  the  phase  angle  between  the  force  and 
velocity,  is  a  function  of  the  frequency  u.  This 
is  not  a  convenient  form  for  the  impedance  ex¬ 
pression.  If  we  use  complex  notation  (where 
Eqs.  (1)  and  (2)  are  the  real  parts  only)  then 

F  =  (4) 

V  =  (5) 

and 


z 


(6) 


where  the  impedance  z  is  now  a  complex  func¬ 
tion  of  the  exciting  frequency  and  is  expressed 
in  polar  form.  The  modulus  or  absolute  value 
Z^  is  the  ratio  of  the  force  magnitude  to  the  ve¬ 
locity  magnitude.  The  argument  le  gives  the 
phase  relationship  between  the  two  quantities, 
with  force  leading  velocity. 

It  is  often  convenient  to  express  the  imped¬ 
ance  in  terms  of  the  real  and  imaginary  parts: 


•  It  should  be  emphasized  that,  although  mechani¬ 
cal  impedance  is  defined  as  the  complex  ratio 
of  sinusoidal  force  to  sinvisoidal  velocity,  five 
other  ratios  between  sinusoidal  force  and  sinus¬ 
oidal  displacement,  velocity,  and  acceleration 
are  equally  useful  and  used.  The  word  imped¬ 
ance  is  often  used  loosely  to  refer  to  all  six 
ratios.  It  is  sometimes  convenient  to  do  this 
in  this  paper. 


Z  =  Zj  +  iZ,  . 

(7) 

Note  that 

Z=  Z„+  iZ,  =  ' 

(Z,/Zr)  .  W 

^  ' 

(9) 

2 


and 

^  "  ‘""’‘(ft)' 

The  real  part  Zg  is  often  called  the  resistive 
part  of  the  impedance  and  the  imaginary  part  Z, 
the  reactive. 

Two  types  of  mechanical  impedance  can  be 
considered:  driving-point  impedance  and  trans¬ 
fer  impedance.  If  the  velocity  is  obtained  at  the 
point  of  application  of  the  exciting  force,  and  in 
the  same  direction  as  the  force,  the  resulting 
ratio  is  called  driving-point  impedance.  If  the 
velocity  is  obtained  at  some  other  point  of  in¬ 
terest  on  the  structure,  or  in  a  direction  unlike 
that  of  the  driving  force,  the  ratio  is  referred 
to  as  transfer  impedance.  The  impedance  is 
usually  expressed  as  Z.^  where  the  second  sub¬ 
script  ic  refers  to  the  point  and  direction  of  ap¬ 
plication  of  the  exciting  force  and  the  first  sub¬ 
script  j ,  to  the  point  and  direction  where  the 
velocity  is  measured.  For  driving-point  im¬ 
pedance  j  =  k ,  and  the  driving-point  impedance 
is  often  expressed  with  the  single  subscript  Zj . 

Impedance  can  be  thought  of  as  a  measure 
of  the  resistance  of  the  system  to  motion.  If 
the  system  is  inherently  difficult  to  move,  stiff 
and/or  massive,  its  impedance  will  be  high. 
Conversely,  if  the  system  is  easy  to  move,  light 
and/or  flexible,  its  impedance  will  be  low. 
However,  since  impedance  is  a  frequency- 
dependent  quantity,  these  generalizations  are 
most  reliable  for  low  frequencies.  In  general  a 
plot  of  impedance  versus  frequency  will  show 
large  fluctuations  in  impedance  magnitude,  with 
a  number  of  peaks  and  valleys.  With  driving- 
point  impedance  the  valleys,  points  of  low- 
impedance,  indicate  ease  of  motion  and  corre¬ 
spond  to  a  resonant  condition  in  the  structure. 
Similarly  the  peaks,  points  of  high  impedance, 
correspond  to  antiresonant  conditions.  With 
transfer  impedance,  these  observations  are 
sometimes  valid. 


TERMINOLOGY  AND  UNITS 

As  mentioned  earlier,  six  ratios  between 
applied  sinusoidal  force  and  resultant  sinusoidal 
displacement,  velocity,  and  acceleration  may 
be  used  to  express  the  same  information.  Be¬ 
cause  certain  usages  are  more  convenient  with 
one  or  the  other  of  these  six  ratios,  examples 
of  most  of  these  may  be  found  in  the  literature 
[9,  11,  12,  16,20-27].  Unfortunately,  the  ter¬ 
minology  has  not  yet  become  standardized,  since 
individual  researchers  have  generated  individual 


terminology.  This  is  even  illustrated  in  the 
other  papers  in  this  group,  where  both  "mobil¬ 
ity"  and  "admittance"  are  used  to  mean  the 
complex  ratio  of  velocity  to  force.  The  Ameri¬ 
can  Standards  Association  recently  has  pub¬ 
lished  recommended  terminology  [28],  which 
should  improve  this  confusing  situation.  Tabic  1 
lists  the  ASA-recommended  terminology,  along 
with  some  alternative  terms  which  have  been 
used  for  the  same  parameters. 

This  table  also  shows  the  relationships  be¬ 
tween  some  of  the  parameters;  the  relation¬ 
ships  which  are  not  shown  (e.g.,  between  me¬ 
chanical  impedance  and  mobility),  might  appear 
to  be  simple  inversions,  but  this  is  not  gener¬ 
ally  true.  For  most  cases,  except  for  very 
simple  ones  such  as  a  pure  mass  system,  the 
inversion  is  complex,  in  that  it  requires  in¬ 
verting  a  matrix  expression  [29].  Because  of 
this,  and  also  because  measurement  or  ana¬ 
lytical  procedures  influence  the  choice  of  pa¬ 
rameter,  there  will  undoubtedly  continue  to  be 
different  parameters  in  use  by  different  inves¬ 
tigators.  However,  Table  1  shows  that,  for  the 
three  parameters  in  which  force  is  always  in 
the  numerator  (or,  conversely,  in  the  denom¬ 
inator),  simple  relationships  exist.  Indeed,  a 
single  plot  on  a  log-log  graph  can  display  all 
three  parameters.  Various  "standard"  graph 
forms  have  been  proposed,  such  as  that  shown 
in  Fig.  1,  as  suggested  by  D.  V.  Wright  [30]  for 
use  in  work  of  the  Structural  Impedance  Panel, 
Ship  Noise  Committee  of  the  Bureau  of  Ships. 
This  form  also  standardizes  on  the  British 
(in- lb- sec)  system  of  units  although  others 
have  used  the  cgs  (centimeter-gram- second) 
and  the  mks  (meter-kilogram- second)  systems. 
While  terminology  and  units  are  technically 
uninteresting,  a  wide  divergence  in  usage  can 
be  very  time-consuming  for  others  trying  to 
use  the  data. 

While  the  preceding  discussion  has  been 
concerned  with  linear  force  and  motion,  pure 
rotational  impedances,  where  force  becomes 
torque  and  linear  motion  becomes  angular  mo¬ 
tion,  or  "mixed"  impedances,  where  one  linear 
quantity  and  one  rotational  quantity  are  used, 
may  be  used  also.  Of  course,  the  complexity 
increases  when  more  degrees  of  freedom  are 
considered. 


SUMMARIZING  REMARKS 

By  this  time,  it  must  be  clear  to  the  reader 
that  mechanical  impedance  is  definitely  not  a 
panacea  for  the  relief  of  overworked  engineers, 
but  is  instead  an  additional  tool  which  requires 
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Sample  graph  showing  Magnitude  on  F/a  and  F/d  Scales 
and  Phase  Angle  on  Degree  Scales  [30] 


TABLE  1 

Factors*  Relating  Impedance  Parameters,  as  Defined  by 
American  Standards  Association  [28] 


Sym¬ 

bol 

Defini¬ 

tion 

V/F 

F/d 

Alternative 

Names 

Dynamic 

Stiffness 

K 

F/d 

x  =  (iai)^ 

X  =  ia> 

x=  1 

Dynamic 

Modulus; 

Mechanical 

Impedance 

Mechanical 

Impedance 

Z 

F/V 

X  =  io) 

X  =  1 

X  =  l/icii 

Immobility 

Effective 

Mass 

M 

F/a 

X  =  1 

X  =  l/io) 

X  =  l/(  ifj)* 

Apparent 

Mass 

??  ? 

a/F 

X  =  1(1^.))^ 

X  =  1/iiti 

X-  1 

Mobility 

Y 

V/F 

X  =  1/iw 

x=  1 

X  =  ici) 

Mechanical 

Admittance 

Dynamic 

Flexibility 

A 

d/F 

X  =  1 

X  =  io; 

x  =  (ia))2 

Mechanical 

Admittance 

■i-Sample  from  line  2  and  column  4:  (F/V)  -  (i(^)  *  (F/a). 


proper  usage  to  yield  good  engineering  results. 
The  other  four  papers  in  this  session  were 
prepared  by  competent  engineers  to  present 

(1)  the  fundamentals  needed  for  analytical  deter¬ 
mination  and  understanding  of  impedance  pa¬ 
rameters,  (2)  the  fundamentals  needed  for  meas¬ 
urement  of  impedance,  (3)  a  review  of  practical 
problems  which  have  used  impedance  in  obtain¬ 
ing  solutions,  and  (4)  a  detailed  procedure  for 
obtaining  rational  design  information  through 
use  of  impedance  information.  In  the  writer's 


opinion,  this  set  of  papers  will  have  accom¬ 
plished  its  purpose  if  it  whets  the  interest  of 
the  reader  sufficiently  so  that  he  (or  she)  will 
consult  other  referenced  works  and  consider 
where  this  tool  can  be  of  use.  Although  incom¬ 
plete,  the  references  and  bibliographies  herein 
contain  much  additional  information.  No  attempt 
has  been  made  to  include  laboratory  reports  on 
particular  projects;  during  the  last  several 
years,  there  have  been  a  considerable  number 
of  these. 
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ANALYTICAL  DETERMINATION  OF  MECHANICAL  IMPEDANCE 


Robert  Plunkett 
University  of  Minnesota 
Minneapolis,  Minnesota 


This  pape  r  discusses  the  impedance,  dynamic  stiffness,  effective  mass ,  mo¬ 
bility,  and  receptance  of  simple  lumped  elements  and  shows  how  they  may  be 
combined  to  find  the  impedance  of  systems  of  moderate  complexity.  It  also 
indicates  methods  for  finding  the  impedance  of  simple,  uniform,  continuous 
systems  and  lists  references  dealing  with  more  complex  structures. 


INTRODUCTION 

Impedance  methods  cannot  replace  other 
types  of  vibration  analysis  but  may  supplement 
them  to  solve  certain  kinds  of  problems  more 
easily.  In  its  simplest  form,  the  impedance 
concept  applies  only  to  steady-state,  sinusoidal 
vibrations  of  linear  structures,  those  in  which 
the  vibration  level  at  any  given  frequency  is 
directly  proportional  to  the  applied  force.  The 
key  to  the  viewpoint  introduced  by  the  impedance 
concept  is  that  it  concentrates  on  the  vibration 
response  of  one  point  of  a  structure  as  a  func¬ 
tion  of  frequency,  whereas  the  usual  types  of 
analysis  determine  the  way  in  which  the  whole 
structure  responds  and  let  frequency  be  a  sec¬ 
ondary  parameter. 

The  major  advantage  one  derives  from  this 
new  viewpoint  is  the  ability  to  find  the  effect  on 
vibratory  response  of  interconnecting  systems 
of  known  characteristics  or  of  modifying  some 
of  the  elements  of  a  given  system.  In  the  stand¬ 
ard  methods  of  analysis,  if  any  part  of  the  sys¬ 
tem  is  changed  in  any  way,  the  complete  anal¬ 
ysis  must  be  repeated  and  each  time  we  find  the 
behavior  of  the  whole  structure;  using  impedance 
methods,  we  must  make  a  separate  analysis  for 
the  ressponse  of  each  point  of  a  complex  struc¬ 
ture  but  a  change  in  one  of  the  elements  re¬ 
quires  a  new  analysis  only  for  that  element. 

This  paper  discusses  some  of  the  methods 
used  for  finding  the  impedance  of  various  kinds 
of  structures  but  does  not  have  space  to  give  all 
the  details  or  many  examples.  A  recent  sym¬ 
posium  [1]  includes  12  papers,  many  of  which 
demonstrate  specific  methods  for  finding  im¬ 
pedances  and  also  gives  an  extensive  bibliogra¬ 
phy  of  other  specific  cases.  Another  recent  set 


of  magazine  articles  [2]  goes  into  the  analysis 
of  some  simple  systems  in  great  detail  with 
numerical  examples.  One  of  the  current  texts 
[3]  shows  the  way  in  which  impedance  methods 
fit  into  the  whole  field  of  vibration  analysis  and 
gives  extensive  tables  of  mode  shapes  and  fre¬ 
quencies  useful  in  this  approach. 


SIMPLE  ELEMENTS 

The  easiest  place  to  start  is  with  the  clas¬ 
sic  building  blocks  of  linear  vibration  theory- 
spring,  viscous  damper,  and  mass— connected 
into  a  damped  single-degree-of-freedom  sys¬ 
tem  (Fig.  la). 


o  bed 


Figure  1 


From  Fig.  la  we  see  that  all  three  ele¬ 
ments  have  the  same  steady- state  amplitude  of 
vibration,  yj  cos  (<i)t  *(p) ,  and  that  P,  cos  is 
the  total  force  applied  to  all  three  of  them;  it 
will  be  remembered  that  we  are  planning  to  con¬ 
sider  only  a  steady- state  sinusoidal  vibration 
of  constant  amplitude  and  to  ignore  any  transients. 
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Tlie  governing  equation  of  kinematic  constraint 
is  that  each  of  the  elements  has  the  same  mo¬ 
tion.  The  equation  of  dynamic  constraint  is: 

Pj  cos  uJt  =  F,  +Fc  +  F„  .  (1) 

The  next  step  is  to  find  the  F's  in  terms 
of  y.  By  definition, 

Fj  =  ky  -  kyj  cos  (2) 

Likewise, 

Fjj  =  cv  =  cy  =  -Ci--yj  sin  (u>t+rf:)  (3) 

and 

“  ma  =  my  =  cos  (wt  +^)  .  (4) 

Now  we  may  expand  Eqs.  (2),  (3),  and  (4)  by  the 
trigonometric  identities: 

sin  (A  +  B)  =  sin  A  cos  B  +  cos  A  sin  B 

and 

cos  (A  +  B)  =  cos  A  cos  B  -  sin  A  sin  B  , 
and  substitute  into  Eq.  (1)  to  get 
Pj  cos  wt  =  (k  cos  <i>  -  c  sin  ^ 

-  mci)^  cos  <}:)  yj  cos 
+  (-k  sin  4  ~  c  cos  >4 


The  most  difficult  and  time-consuming  part  of 
this  development  has  been  the  algebraic  manip¬ 
ulation.  The  labor  can  be  reduced  by  using 
phasors  which  rely  on  Euler's  formula, 

=  cos  i/.  +  j  sin  v  .  (3) 

where  =  -l.  If  we  now  let  the  sinusoidal 
force  be 

P  =  P,  , 

y  becomes 


where  y  is  a  complex  number,  including  the 
phase  angle  and  yj  is  a  real  number,  the 
modulus  or  absolute  value  of  y  If  these  are 
now  substituted  into  the  above  relationships. 


F,  =  ky  =  kye'"',  (2a) 

F|3  =  cv  =  cy  =  jcyicc*"^  ,  and  (3a) 

Ff^  =  ma  =  my  =  -my  e  ^ .  (4a) 

If  these  are  substituted  into  Eq.  (1), 

Pj  e*"’  =  (k  +  jccc  -  yc’'"'  .  (5a) 


The  tangent  of  the  phase  angle  is  the  negative 
of  the  ratio  of  the  real  to  imaginary  parts  of 
the  coefficient, 


+  mc<;^  sin  4)  Vj  sin  «Jt  .  (5) 


tan 


(6a) 


When  sin  wt  =  1  ,  then  cos  <jt  =  0  ,  so  that 
(-k  sin  v'  -  c.»J  cos  4  +  mij^  sin  <^)  =  0  . 


or. 


and  the  absolute  value  is  the  square  root  of  the 
sum  of  the  squares,  or 


P,  =  [(k 


Vl 


(7a) 


t.m 


4 


cai 

mij^  -  k 


(6) 


at  a  great  saving  in  algebra. 


Since 


It  is  also  convenient  to  nondimensionalize 
Eqs.  (6)  and  (7)  by  letting 


2  ,  tan^  </  2 

5  111-“./  = - s—  and  cos  ■‘4 


1  +  4 


1  +  tan^ 4 


Pj  cos  .ct  -  \4k  -  ^  y  j  cos  ,  (7) 


u>^  =  k/m 

4'  = 

C,.  =  2/v.m  =  2m  =  2k  =  2Z^  (9) 


after  a  certain  amount  of  algebraic  manipulation. 


,o 


C  C 


c 


Of  necessity.  Eqs.  (6)  and  (7)  are  exactly 
those  one  finds  in  any  text  in  mechanical  vibra¬ 
tions  for  the  steady-state  response  and  phase 
angle  of  a  single-degree-of-freedom  system. 


Vo  =  Pi  'k  . 


Then, 


tan  <y 


2pr 

1  -  .:2 


(6b) 
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and 


We  may  rewrite  Eq.  (7b)  either  as 


yi  =  Vo  +  (2p/3)2]  .  (7b) 

It  will  be  seen  that  both  4-  and  are  func¬ 

tions  of  /3  onlyi  if  p  is  considered  as  a  param¬ 
eter,  the  dimensionless  damping.  We  may  thus 
make  two  plots— ^  versus  yS  and  yj/y^  versus 
/3— with  p  as  a  parameter  and  have  the  re¬ 
sponse  curves  for  all  single-degree-of-freedom 
systems. 


PRESENTATION  OF  RESULTS 

Very  often  for  engineering  work,  however, 
we  are  interested  in  actual  frequencies  and  the 
numerical  values  of  spring,  mass,  and  damping; 
in  this  case,  we  must  multiply  /3  by  (or  ) 
and  p  by  to  get  dimensional  values.  It  has 
recently  become  customary  to  plot  log  Vl/Vo 
and  4  versus  log  /J  so  that  dimensional  values 
may  be  found  by  relabeling  the  graphs  without 
changing  the  shape  of  the  curves;  Fig.  2  shows 
nondimensional  curves  and  the  results  for 
f|,  =  10  cps. 


V 


or 


(7c) 


(7d) 


While  any  of  the  forms  (7b),  (c),  or  (d),  are 
equally  useful,  the  author  prefers  (7c)  because 
it  is  symmetrical  in  /3  and  l,  p  and  it  has 
some  other  technical  advantages.  Therefore 
this  paper  will  henceforth  present  all  develop¬ 
ments  and  results  in  terms  of  velocity  as  in 
Fig.  3  which  is  a  replot  of  Fig.  2. 

Heretofore,  all  results  have  been  given  in 
terms  of  reference  displacement  y„,  velocity 


Figure  Z 


Figure  3 


The  development  just  presented  and  Fig.  2 
give  the  result  for  vibration  displacement.  Al¬ 
though  it  was  customary  for  many  years  to  con¬ 
sider  mainly  vibration  displacement,  it  is  not 
the  only  quantity  of  interest  since  vibration 
velocity  controls  both  stress  and  noise  [4]  and 
vibration  acceleration  is  controlled  by  unbal¬ 
ance  [5].  It  is  easy  to  get  one  from  another 
since  if 


j  -jt 

y=yif 

V  “  y  =  j  ■  y  y. 

(I  =  y  =  j.v  -  - 


(10) 


,  or  acceleration  as  was  pointed  out  be¬ 
fore,  dimensionless  charts  are  not  directly  use¬ 
ful  for  engineering  applications  where  we  usu¬ 
ally  want  to  know  how  much  vibration  we  get 
per  unit  force.  The  impedance  is  the  quantita¬ 
tive  measure  of  how  well  a  vibratory  system 
resists  a  force;  it  is  the  force  per  unit  velocity. 
Thus  for  this  system. 


Pi 


V 

o 


yiin 


Zo 


[Eq.  (9)1.  (12) 
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Borrowing  still  further  from  circuit  theory, 
l/2p  is  often  known  as  the  quality  factor  Q. 
Thus,  finally  the  impedance  for  a  single- 
degree-of-freedom  system  is 


CHAINS  OF  SIMPLE  ELEMENTS 

We  have  defined  impedance  as  the  ratio 
(complex)  of  force  to  velocity,  or 


Z  = 


(13) 


(15) 


The  mobility  is  a  measure  of  how  mobile 
the  system  is,  how  easily  it  moves,  or  the  ratio 
of  velocity  to  force 


M  = 


V 

P 


z‘ 


(14) 


Other  ratios  which  have  been  used  are:  re- 
ceptance  [3],  ratio  of  displacement  to  force; 
dynamic  stiffness  [6],  ratio  of  force  to  dis¬ 
placement;  effective  mass,  ratio  of  force  to 
acceleration;  and  acceleration  mobility  [5], 
ratio  of  acceleration  to  force.  Figure  4  is  a 
plot  of  the  mobility  of  a  single-degree-of- 
freedom  system  with  a  natural  frequency  of  10 
cps,  a  Q  of  10  or  p  of  0.05  and  a  mass  which 
weighs  10  pounds.  Figure  5  is  a  plot  of  the 
impedance  of  the  same  system. 


Figure  4 


Figure 


Referring  to  Fig.  1  and  Eq.  (1),  we  see  that  all 
the  elements  of  a  single-degree-of-freedom 
system  have  the  same  velocity  and  the  total 
force  is  the  sum  of  the  forces.  Thus  if  we  di¬ 
vide  Eq.  (1)  by  the  common  velocity, 


V  V  V  V 


or 

=  z.  +  Zd  +  Z.  •  (16) 

Equation  (16)  says  that  the  impedance  of  this 
system  is  the  sum  of  the  impedances  of  the  ele¬ 
ments  because  the  velocity  is  common.  In  Fig. 
6  the  force  on  the  mass  and  the  spring  are  the 
same  and  the  velocity  of  the  end  point  of  the 
spring  is  the  sum  of  that  due  to  the  mass  and 
that  due  to  the  spring,  or 


and 

M,  =  M,  +  ,  (17) 


Figure  6 


Whether  we  add  impedances  or  mobilities 
depends  on  whether  the  forces  add  and  the  ve¬ 
locity  is  the  same  (mechanical  parallel)  or  the 
velocities  add  and  the  force  is  the  same  (me¬ 
chanical  series).  It  is  possible  to  construct  the 
impedance  of  a  lumped  system  by  successively 
adding  the  elementary  components. 

Let  us  consider  the  case  of  torsional  vibra¬ 
tion  of  a  shaft  with  disks  distributed  along  it. 

Fig.  7.  We  may  start  from  the  left  end  and  build 
up  the  impedances  until  we  get  to  the  right  end. 
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Thus, 


s,  e, 
/> '  r*'  />’  /I  ^ 

'n  'm  m 


U  k, 

Ji 


''23  '  ’'3 

Jj  J 


Figure  7 


ZjU  -  Zj  +  ZjL 


=  j-jj  +  -j: 


1 


‘12 


(23) 


We  may  proceed  as  above  to  the  left  side  of  Jj, 


Let  Z,  be  the  impedance  of  the  disk  Jj  which 
is,  by  definition,  the  ratio  of  the  applied  torque 
to  the  resulting  vibratory  angular  velocity:  then, 

^IR  ”  j“'Ji  •  (18) 

The  vibratory  velocity  of  is  the  vibratory 
velocity  of  the  left  end  of  kj^,  which  is  the  right 
end  of  Zj,  plus  the  increase  in  vibratory  veloc¬ 
ity  across  kjj. 


Vj  =  Vj  +  Vi2  .  (19) 


''.IL  =  +  ■'2R 


-  J± 


^^23 


1 

i.  ’ 


i  1 


''3L 


JlJ  1 

- 1 


^  JiiL  A 


(24) 


j^Ji 


However,  the  torque  applied  to  the  right  end  of 
k,2  is  the  torque  applied  to  the  right  end  of  j,  , 
or 


f  12 

^12 

Ti2 

"2L 

■  ^*1R 

+  M,2 

But, 


M 


IR 


1 


M 


12 


Thus, 


Now, 


M 


2L 


j.o 


k 


1  2 


+ 


1 


1  1 


M 


2L 


j  '■ 


+ 


i",!  1 


(20) 


(21) 


(22) 


The  velocity  of  is  the  same  as  that  of  the 
right  end  of  k,^  and  the  torque  applied  to  the 
right  end  of  of  11*^1  applied  to 

k,i  and  that  applied  to  ,l2,or 


and  to  the  right  side  of  jj, 

^3R  ■ 


=  + 


j<'^J2  + 


I 


-ill  +  1 

'<12  j  J 1 


(25) 


The  pattern  is  now  obvious,  and  we  may  write 
for  the  impedance,  looking  at  the  right  end  of 
J4  > 


Z4R  =  i“il4  + 


<34 


j-Jj  +  i 


<23  ,  1 

+  -p;;  —T 


(26) 


‘12 


j“tli 


The  approach  is  obviously  the  same  as  that  of 
the  Holzer  method,  except  that  frequency  is 
carried  as  a  parameter  instead  of  being  evalu¬ 
ated  numerically.  The  use  one  makes  of  the 
final  expression  depends  on  the  application  and 
relative  magnitude  of  the  terms;  one  might  use 
the  theorems  of  continued  fractions  |7]  to  find 
the  characteristic  equation  containing  the  four 
natural  frequencies  (one  of  them  is  zero)  or  one 
may  use  various  approximations  to  find  the  im¬ 
pedance  over  a  limited  frequency  range. 

The  development  just  presented  shows  that 
impedance  and  mobility  are  equally  important 
in  our  analysis.  Which  we  use,  depends  on 
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whether  we  are  joining  elements  so  that  they 
see  the  same  velocity  and  add  forces,  or  see 
the  same  force  and  add  velocities.  In  the  first 
case,  we  speak  of  mechanical  parallel  and  add 
impedances  (Fig.  8a);  in  the  second,  we  speak 
of  mechanical  series  and  add  mobilities  (Fig.  8b). 


^ — 1^1  I* — *  F 


Strips 


b 


Figure  8 


COMBINING  SUBSYSTEMS 

It  was  stated  in  the  Introduction  that  im¬ 
pedance  methods  differ  from  previous,  more 
classical  methods  in  focussing  attention  on  the 
vibration  response  of  one  point  as  a  function  of 
frequency.  It  is  this  aspect  of  the  approach  that 
is  used  for  combining  subsystems.  This  may 
be  illustrated  by  considering  the  classical  case 
of  the  damped  dynamic  vibration  absorber  at¬ 
tached  to  a  single-degree-of-freedom  system 
[8]  (Fig.  9).  This  system  may  be  considered  in 
two  parts  (Fig.  10). 


^WWVWvj'n  i  I 


lAWWVWi 


Cz 


Figure  9 


The  governing  equations  are  clearly, 


V  =  Vj  ^  Vj  , 

and 

F  =  F,  +  Fj  . 


Thus, 


or 

Z  =  Zi  +  Zj  ,  (27) 

since  the  velocity  is  the  same  for  both  pieces 
and  the  forces  add  (mechanical  parallel,  Fig. 

8a).  We  have  already  seen  that 

=  j  (^^1  -  ~)  ■ 

And  it  is  not  difficult  to  show  that 

^2  ■  71  2  rn  ■ 

I  K  2  "  cj*  mj  )  +  ^*2 

If  Eqs.  (28)  and  (29)  are  substituted  into  Eq.  (27), 


z 


1  + 


)u!  Cj 


1 


(30) 


This  development  may  be  compared  with  that 
used  in  the  classic  approacli  [8],  One  might 
point  out  some  conclusions  that  may  be  more 
obvious  in  this  formulation  than  in  the  usual 
one.  If  there  is  no  absorber,  Zj  becomes  zero 
at  =  kj/mj  and  amplitude  becomes  exces¬ 
sive;  from  Eq.  (30)  we  see  that  z  can  never  be 
zero  unless  cj  is,  in  which  case  there  are  two 
frequencies  of  resonance,  one  lower  than  either 
/kj/m,  or  v/k^ntj  and  one  higher  than  either. 

The  larger  the  mass  ratio,  m^/mj ,  the  more 
these  two  new  resonant  frequencies  differ  from 
the  old  ones.  Such  a  device  is  used  for  a  system 
driven  at  a  constant  frequency  to  increase  the 
impedance  and  reduce  the  response  at  that  fre¬ 
quency:  this  means  that  Zj  should  be  as  large 
as  possible  at  the  frequency  it  is  desired  to  sup¬ 
press.  It  is  clear  from  Eq.  (29)  that  this  fre¬ 
quency  is  very  close  to  \/k  j/mj  and  that  the 
maximum  value  of  Zj  decreases  as  in¬ 
creases;  cj  is  useful  only  to  control  amplitude 
at  the  two  new  resonances  and  not  at  v^kj.  . 

If  these  conclusions  were  made  quantitative, 
they  would  confirm  the  major  results  of  previ¬ 
ous  investigations.  What  is  more  interesting, 
these  conclusions  are  not  confined  to  vibration 
absorbers  used  on  single-degree-of-freedom 
systems.  The  major  point  of  interest  is  that 
Eqs.  (27)  and  (29)  still  hold,  but  Eq.  (28)  no 
longer  gives  the  impedance  of  the  primary  sys¬ 
tem.  If  we  can  find  some  expression  for  Z, , 
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M,j  F,  +  MijFj 


analytic  or  measured,  we  may  adjust  kj/mj, 
cticj/k]  and  ni]  for  adequate  control  at  any  de¬ 
sired  frequency. 


V,  = 

V2  =  Mj,F,  +  MjjFj  (31) 


MULTIMASS  SYSTEMS 

Up  to  this  point,  we  have  discussed  the  vi¬ 
bration  response  of  one  point  of  a  system  when 
it  is  excited  by  a  sinusoidal  force  of  the  same 
character  as  the  motion  acting  at  the  same 
point.  In  the  case  of  a  rotation  (Fig.  7),  the 
force  must  be  a  torque  about  the  same  axis;  in 
the  case  of  linear  motion,  the  force  must  be  in 
the  same  vectorial  direction  as  the  motion.  To 
state  it  more  rigorously,  the  force  must  be  that 
generalized  force  corresponding  to  the  gener¬ 
alized  coordinate  of  the  Lagrangian.  When  we 
deal  with  cases  where  more  than  one  point,  or 
motions  at  one  point  in  more  than  one  direction, 
are  involved,  we  must  use  somewhat  more  com¬ 
plicated  expressions.  Suppose  we  are  concerned 
with  a  system  of  many  masses,  springs,  and 
dampers  and  that  the  reactions  or  constraints 
are  fully  defined.  Then,  if  the  system  be  linear, 
the  vibration  response  to  any  one  sinusoidal  vi¬ 
bratory  force  is  sinusoidal,  of  the  same  fre¬ 
quency  and  proportional  to  the  force.  The  am¬ 
plitude  of  any  one  point  is,  in  general,  different 
from  that  of  any  other  so  that  we  may  write 


vj  =  MjjF, 

^3  =  MjiFj.  etc. 

In  these  equations,  Mj , ,  Mji ,  M31 ,  and  the  like, 
are  complex  numbers  giving  the  ratio  of  the  vi¬ 
bration  velocity  at  points  1,  2,  3,  etc.,  to  the  vi¬ 
bratory  force  at  point  1.  Referring  to  our  pre¬ 
vious  discussion,  they  are  functions  of  frequency, 
and  may  be  called  mobilities;  is  called  a 
driving  point,  or  input  mobility  and  Mj,  or  Mj, 
are  called  transfer  mobilities.  If  we  now  have 
a  force  operating  on  point  2  alone  and  none  on 
point  1,  we  may  write 

v^  =  MjjFj 

Vj  “  MjjFj 

V  3  =  M32  ^2  •  etc. 

Now,  if  both  forces  act  at  the  same  time  and 
with  the  same  frequency,  the  vibration  ampli¬ 
tudes  will  be  sum  of  those  due  to  each  force 
separately: 


''3  "  **31  *’'1  **32*^2  ■ 

In  expressions  of  this  type,  Mij  is  the  velocity 
response  of  point  i  due  to  a  unit  vibratory  force 
at  point  j  with  all  other  forces  except  the  stated 
reactions  or  constraints  being  zero;  it  may  be 
shown  that  for  all  such  passive  mechanical  sys¬ 
tems,  Mij  =  Mji  (dynamic  reciprocity)  [9,10]. 
When  we  start  defining  the  more  general  form 
of  impedance,  we  find  the  situation  somewhat 
more  difficult.  Since  the  vibration  velocities 
are  caused  by  the  vibratory  forces,  we  may  not 
use  more  independently  defined  velocities  than 
forces.  Thus  in  the  equations  similar  to  Eqs. 
(31),  we  must  have  as  many  equations  for  force 
as  we  have  given  velocities. 


*^1  '  ^11  ''1  ^  ^12  ''2 
*^2  =  ^21  ''1  *  ^22  ''2 


(32) 


if  there  are  only  two  forces.  In  this  case,  Zij 
is  the  force  at  point  i  due  to  a  unit  velocity  of 
point  j  when  all  other  velocities  corresponding 
to  forces  being  considered  are  held  zero;  that 
is,  all  other  forces  must  be  replaced  by  rigid 
constraints.  Under  these  conditions,  it  may  be 
shown  that  reciprocity  still  holds,  Zij  =  Zj  i ; 
however,  it  may  no  longer  be  said  that  the  im¬ 
pedances,  Z ,  are  the  reciprocals  of  the  mobili¬ 
ties  M.  Instead,  one  must  say  that  the  square 
matrix  of  impedances  is  the  inverse  of  the  cor¬ 
responding  square  matrix  of  the  mobilities  [11]. 
It  should  be  noted  that,  in  general,  the  values  of 
the  Z’s  will  depend  on  how  many  forces  are 
being  considered,  whereas  the  values  of  the  M’s 
do  not.  As  for  simpler  systems,  the  impedance 
or  stiffness  formulation  must  be  used  for  some 
analyses  [6, 12]  and  mobility  or  receptance  for¬ 
mulation  for  others  [3, 13]. 

Starting  with  these  formulations,  we  may 
set  up  general  expressions  for  the  Mobilities 
or  Impedances  in  terms  of  the  resonant  fre¬ 
quencies  and  mode  shapes.  While  these  proce¬ 
dures  have  been  known  to  be  valid  since  the 
time  of  Fourier  [9],  it  has  been  only  since  the 
development  of  high-speed  digital  computers 
that  they  have  become  practical  [14, 15].  Since 
these  methods  (Given  or  Jacobi)  are  now  read¬ 
ily  available,  we  will  assume  that  all  of  the 
normal  modes  and  frequencies  of  an  undamped 
multimass  system  are  known.  It  may  be  shown 
that  the  mobilities  can  be  expanded  as  a  series 
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where  n  is  the  number  of  masses,  natural  fre¬ 
quencies,  and  mode  shapes.  The  a's  may  be 
found  from  the  mode  shape  [14];  the  easiest  ap¬ 
proach  is  by  matrix  multiplication,  an  operation 
well  suited  to  digital  calculators.  If  Eq.  (33)  is 
converted  into  one  fraction,  the  denominator  is 
the  product  of  the  individual  denominators  and 
the  numerator  is  a  polynomial  of  degree  n  -  1 
in  01^,  or 


(34) 


where  'U] ,  , ,  and  the  like,  are  the  antireso¬ 

nant  frequencies  of  the  point.  This  type  of  ex¬ 
pression  is  known  as  Foster's  Theorem  (16). 

In  this  case,  the  denominator  is  the  same  for 
all  points  and  depends  only  on  the  resonant  fre¬ 
quencies  of  the  system;  the  constant  and  the 
numerator  are  different  for  all  points  since  the 
antiresonant  frequencies  are  point  and  not  sys¬ 
tem  parameters. 

The  impedance  functions  of  Eqs.  (32)  may 
be  expressed  in  a  similar  way,  except  that  now 
the  zeros  are  system  parameters  and  the  poles 
are  the  point  parameters.  One  of  the  more  use¬ 
ful  aspects  of  Foster's  theorem  is  that  the  mo¬ 
bility  or  impedance  may  be  expressed  quite  ac¬ 
curately  over  a  limited  frequency  range  by 
using  only  the  resonant  or  antiresonant  fre¬ 
quencies  that  fall  within  the  range  [6].  Figure 
11  is  a  schematic  representation  of  mobility  or 
impedance  over  such  a  limited  range  [2, 5, 14]. 

CONTINUOUS  SYSTEMS 

Foster's  theorem  or  the  expansion  in  nor¬ 
mal  modes,  Eq.  (33),  lets  us  express  the  mobility 


Figure  1 1 


or  impedance  of  continuous  systems  in  terms 
of  a  finite  number  of  frequencies.  In  fact,  when 
we  are  talking  of  real  systems,  this  is  what  we 
always  do;  all  real  systems  have  an  infinite  num¬ 
ber  of  natural  frequencies  so  that  only  finite  ex¬ 
pression  is  only  on  approximation.  We  may 
also  solve  such  systems  by  approximating  them 
with  lumped  springs  and  masses  and  solving 
the  reduced  physical  problem. 

The  arithmetical  labor  of  this  approach  is 
reduced  by  matrix  formulation  [17, 18]  and  by 
extensive  tables  of  the  mode  sh^es  and  fre¬ 
quencies  of  simple  shapes  [3].  Impedance  meth¬ 
ods  may  also  be  useful  for  analyzing  nonlinear 
vibration  problems  where  most  of  the  system 
is  linear  and  the  nonlinearity  involves  only  a 
few  elements  [19]. 


DAMPING 

The  phase  and  amplitude  effects  of  damping 
are  discussed  at  some  length  in  several  of  the 
papers  of  Ref.  1,  notably  Refs.  12  and  16,  and 
these  precise  formulations  are  often  useful  in 
servomechanism  or  feedback  problems.  For 
most  vibration  problems,  small  damping  is 
most  easily  handled  by  assuming  an  amplifica¬ 
tion  factor  for  each  mode  shape  as  in  Figs.  4 
and  5.  It  can  be  shown  that  this  will  also  con¬ 
trol  the  antiresonances  in  a  similar  fashion  [20]. 

If  the  damping  is  not  linear  but  still  small, 
it  may  be  handled  as  Jacobsen's  equivalent 
viscous  damping  in  the  manner  described  in 
most  vibration  texts.  (Also  on  damping,  "Anal¬ 
ysis  of  Dynamic  Systems  Using  the  Mechanical 
Impedance  Concept,"  W.  C.  Sperry,  Noise  Con¬ 
trol,  7,  2,  March-April  1961,  pp.  13-21.) 
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DISCUSSION 


A  question  was  asked  about  the  electrical 
analogy. 

Dr.  Plunkett:  I  was  afraid  somebody  was 
going  to  ask  that  question.  The  answer  is  yes 
and  no.  I  do  not  like  to  think  of  this  tool  in 
terms  of  an  analogy.  Because  if  you  are  going 
to  try  and  set  up  an  analogy  you  can  set  it  up  in 
a  number  of  different  ways  depending  on  the  ap¬ 
plication.  From  my  standpoint,  I  find  it  easier 
to  think  of  it  in  terms  of  its  primitive  defini¬ 
tion,  in  terms  of  a  ratio  between  a  sinusoidal 
force  and  the  resulting  sinusoidal  motion.  Now 
it  depends,  if  you  are  going  to  call  it  admittance. 


it  depends  on  whether  you  call  your  analogy  be¬ 
tween  motion  and  voltage,  or  between  motion 
and  current.  You  can  set  up  your  analogy  either 
way  and  this  then  determines  whether  you  call 
it  admittance  or  not,  which  is  why  I  personally 
dislike  the  use  of  the  word  admittance.  I  think 
it  is  unfortunate  that  we  have  used  the  word 
impedance,  that's  bad  enough.  Does  that  suc¬ 
cessfully  evade  your  question? 

Mr.  Saunders  (General  Electric):  One 
thing  t  would  like  to  mention.  I  think  Dr. 
Plunkett  has  done  a  great  and  very  interesting 
task  on  this.  One  thing  he  has  neglected  to 


I 


16 


mention  and  that  is  the  application  of  matrix 
methods.  Now,  if  you  have  a  very  large  and 
complex  system,  it's  very  easy  to  use  the 
matrix  method;  in  fact,  you  let  the  computing 
machine  do  all  the  work  for  you.  When  you 
have  the  matrix  method,  if  you  have  to  break 
into  it  as  Dr.  Plunkett  mentioned,  why  you  can 
actually  store  the  matrix  in  the  machine,  you 
can  actually  start  from  there  and  put  your  new 
system  on,  or  add  anything  else  on  to  it.  Ac¬ 
tually  the  machine  will  do  most  of  your  work 
for  you.  One  of  the  other  things  I  would  like  to 
mention  also  is  about  the  Myklestaad  method. 
Now  if  you  have,  say,  a  non-uniform  beam  like 
a  missile,  you  can  actually  determine  your 
transfer  impedances  by  using  the  Myklestadd 
procedure  very  easily,  and  simply.  Therefore, 
by  using  a  computer  and  using  some  of  these 
other  tools  you  can  actually  reduce  your  work 
down  to  a  bare  minimum.  Let  the  machine  do 
all  the  work  for  you. 

Dr.  Plunkett:  I  can  only  agree  with  what 
he  says,  and  1  can  only  reiterate  what  I  said 
earlier— that  impedance  methods  do  not  re¬ 
place  any  other  method  that  is  useful  in 
vibration— and  you  can  make  use  of  any  of  them 
for  certain  applications. 

Mr.  Stern  (General  Electric):  As  you  prob¬ 
ably  know,  many  times  when  items  are  being 
evaluated  to  meet  a  particular  test  their  failure 
rates  can  be  rather  excessive  due  to  the  specs 
calling  for  a  rather  rigid  fixture.  Whereas,  in 
the  actual  application,  the  supporting  struc¬ 
tures  are  usually  quite  soft,  and  even  for  the 
same  input,  you  usually  don't  get  the  same 
severity  of  damage.  I  wondered  how  you  felt 


about  the  possibility  of  having  impedances,  fix¬ 
ture  impedances  called  down  in  specifications? 
Particularly  on  missile  equipment  or  satellite 
equipment  where  weight  is  important. 

Dr.  Plunkett:  Well,  I  think  you  may  be  ask¬ 
ing  the  wrong  person  this  because  Dr.  Belsheim 
is  the  authority.  Let  me  say  it  though  because 
he  may  not  want  to.  The  people  at  the  U.  S. 
Naval  Research  Laboratory  in  recent  years  have 
demonstrated  quite  conclusively  that  this  lack 
of  appreciation  of  the  influence  of  support  im¬ 
pedance  on  shock  level,  in  particular,  has  led 
to  shock  specifications  being  considerably  more 
severe  than  they  ought  to  be.  I  am  not  sure  he 
would  want  to  be  quoted  on  a  particular  number. 
My  personal  impression  is  that  they  are  at  least 
five  times  too  high,  but  that  is  not  an  official 
opinion,  fortunately.  I  would  think  that  one  of 
the  places  where  we  could  get  most  mileage  out 
of  this  kind  of  approach  is  by  making  sure  that 
our  specification  does  indeed  match  the  environ¬ 
ments  to  which  we  think  that  we  are  trying  to 
subject  the  equipment.  Very  often,  we  have  got¬ 
ten  in  exactly  that  kind  of  a  trap.  This  is  the 
converse  problem  of  course  of  the  isolator  ef¬ 
fectiveness,  where  isolators  turn  out  not  to  be 
effective  because  they  are  not  supported  prop¬ 
erly.  Likewise,  if  we  make  a  measurement  and 
then  don't  support  things  properly,  why  we  will 
get  a  lower  vibration  or  shock  environment  than 
we  expected.  If,  on  the  other  hand,  our  ma¬ 
chines  are  carefully  ^ecified  to  be  extremely 
rigid,  our  environment  will  turn  out  to  be  much 
higher  than  we  thought.  The  only  answer  I  can 
give  to  that  is  I'm  also  in  favor  of  motherhood 
and  against  sin.  .  .firmly.  .  .oh!  I'm  an  anti¬ 
communist  too. 
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INSTRUMENTS  AND  METHODS  FOR  MEASURING 
MECHANICAL  IMPEDANCE 


R.  R.  Bouche 
Endevco  Corporation 
Pasadena,  California 


After  mentioning  some  of  the  various  methods  used  for  experimentally 
measuring  mechanical  impedance,  this  paper  describes  the  performance 
characteristics  that  affect  the  accuracy  of  recently  developed  meciiani- 
cal  impedance  heads.  The  results  of  calibrations  and  evaluation  tests 
performed  on  an  impedance  head  are  described.  Tests  on  a  simple 
structure  are  made  to  illustrate  the  suitability  of  the  liead  for  malting 
impedance  measurements. 


INTRODUCTION 

Various  methods  have  been  employed  for 
the  experimental  determination  of  the  mechani¬ 
cal  impedance  of  structures. 

An  indirect  method  [l]  employs  an  electro¬ 
dynamic  vibration  exciter  equipped  with  both  a 
driving  coil  and  velocity  sensing  coil.  The 
transfer  electrical  impedance  of  the  coils  is 
measured  with  known  mechanical  impedances 
attached  to  the  exciter  and  with  the  structure  of 
unknown  impedance  attached.  The  mechanical 
impcuvmce  of  the  structure  is  computed  from 
equations  relating  it  to  the  measured  electrical 
impedances.  Other  indirect  methods  for  meas¬ 
uring  mechanical  impedance  include  (a)  meas¬ 
uring  the  terminal  electrical  impedance  of  a 
variable  reluctance  transducer  under  blocked 
and  free  conditions  and  (b)  measuring  the 
change  in  response  of  two  identical  reeds  of 
known  mechanical  impedance  when  one  is  con¬ 
nected  to  the  structure  of  unknown  impedance. 
These  methods  and  others  are  briefly  discussed 
by  Coleman  [2]. 

A  transient  method  also  has  been  used  for 
measuring  mechanical  impedance  [3].  A  tran¬ 
sient  force  is  applied  and  the  resulting  motion 
of  the  structure  is  measured  with  an  acceler¬ 
ometer.  The  impedance  is  determined  from 
"the  ratio  of  the  Fourier  spectra  of  the  force 
and  response  histories." 

The  most  direct  technique  to  determine 
mechanicai  impedance  is  to  utilize  force  and 


motion  measuring  transducers.  Recently  a 
number  of  experimentists  [2-10]  have  employed 
this  technique  on  structures.  Since  this  tech¬ 
nique  is  quite  widely  accepted,  the  remainder 
of  this  paper  is  devoted  to  a  calibration  and 
evaluation  of  an  impedance  head  and  its  use  in 
measuring  the  impedance  of  a  simple  structure. 
A  mechanical  impedance  head  is  a  transducer 
with  built-in  force  and  motion  measuring  de¬ 
vices.  Although  one  type  of  impedance  head  is 
herein  described  in  detail,  the  calibration  pro¬ 
cedures  and  discussion  of  accuracies  are  ap¬ 
plicable  to  all  transducers  and  their  use  in 
making  mechanical  impedance  measurements. 


DESCRIPTION  OF  IMPEDANCE  HEAD 

Mechanical  impedance  is  defined  as  the 
ratio  of  the  force  exciting  a  structure  to  the 
resulting  velocity.  The  motion  transducer  in 
the  impedance  head  may  be  a  displacement,  ve¬ 
locity,  or  acceleration  transducer.  Plotting  the 
apparent  stiffness  (ratio  of  force  to  displace¬ 
ment)  or  apparent  weight  (ratio  of  force  to  ac¬ 
celeration)  on  suitable  graph  paper  provides 
directly  the  magnitude  of  the  mechanical  im¬ 
pedance  when  the  excitation  is  sinusoidal. 
Sinusoidal  excitation  also  makes  it  possible  to 
utilize  a  variety  of  instruments  for  accurately 
measuring  the  phase  angle  of  the  impedance. 

Most  impedance  measurements  are  made 
in  the  frequency  range  from  10  cps  to  5000  cps. 
In  order  to  have  wide  application,  the  imped¬ 
ance  head  transducers  should  be  flat  (free  of 
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frequency  distortion)  and  have  zero  phase  angle 
between  the  force  and  motion  transducers 
throughout  this  frequency  range.  If  the  phase 
angle  is  not  zero,  the  data  may  be  corrected,  if 
the  phase  shift  is  known  and  the  excitation  si¬ 
nusoidal.  The  transducers  must  also  be  ampli¬ 
tude  linear  in  order  to  measure  a  wide  range  of 
impedances.  Mechanical  impedance  measure¬ 
ments  usually  are  made  at  forces  ranging  from 
less  than  1  pound  up  to  about  100  pounds,  and 
accelerations  up  to  several  hundred  times  the 
acceleration  of  gravity.  Even  very  heavy  struc¬ 
tures  can  be  measured  in  the  low  force  ranges 
since  the  frequencies  of  primary  interest  are 
those  above  which  the  structure  acts  as  a  rigid 
body.  Significant  motions  result  from  small 
forces  near  resonances  (i.e.,  regions  of  low 
impedance).  The  impedance  head  must  also 
measure  motions  corre^onding  to  accelera¬ 
tions  much  less  than  1  g  (386  in.  per/'sec^)  in 
order  to  measure  accurately  high  impedances 
at  antiresonances  in  the  structure.  For  this 
reason,  accelerometers  are  preferred  for  mak¬ 
ing  impedance  measurements  at  high  frequen¬ 
cies,  since  the  outputs  from  velocity  and  dis¬ 
placement  pickups  are  quite  small  for  small 
amplitude  motions  at  frequencies  much  above 
1000  cps. 

Another  requirement  for  accurately  meas¬ 
uring  large  impedances  is  that  the  motion 
transducer  should  not  respond  to  force  excita¬ 
tion  particularly  when  the  acceleration  tends  to 
zero.  If  false  outputs  result  from  the  motion 
transducer  due  to  dynamic  stresses  applied  to 
the  case  of  the  impedance  head,  the  indicated 
impedance  would  less  than  the  actual  imped¬ 
ance.  The  impedance  head  should  be  designed 
so  this  interaction  does  not  occur. 

On  the  other  hand,  the  interaction  of  the 
force  transducer  responding  to  applied  acceler¬ 
ations  normally  occurs  in  an  imp^ance  head. 
This  output  from  the  force  tran^ucer  equals 
the  mass  of  the  end  of  the  impedance  head  plus 
any  externally  attached  mass  times  their  accel¬ 
eration.  This  characteristic  is  used  to  dynam¬ 
ically  calibrate  the  force  transducer  in  the 
impedance  head. 

The  mass  of  the  impedance  head  should  be 
small  and  its  stiffness  large  so  that  the  effect 
of  the  impedance  head  on  the  re^nse  of  the 
structure  is  not  significant.  If  the  head  has 
these  characteristics,  the  measured  data  on 
many  structures  need  not  be  corrected  for  the 
head  mass  and  stiffness. 

One  type  of  impedance  head  and  an  accom¬ 
panying  force  generator  are  shown  in  Fig.  1. 
The  impedance  head  is  1  inch  high  and  nearly 


Fig.  1  -  Electrodynamic  force  generator 
and  impedance  head  used  for  measuring 
the  impedance  of  structures 


2  inches  in  diameter.  It  has  three  piezoelectric 
force  pickups  and  three  piezoelectric  acceler- 
ometsrs  built  into  the  head.  A  bolt  is  attached 
to  the  force  generator  and  passed  through  the 
hole  in  the  center  of  the  head  to  rigidly  connect 
the  head  and  outer  casing  of  the  impedance  head 
to  a  structure.  The  force  generator  consists  of 
an  electrodynamic  vibration  exciter  with  its 
driving  coil  rigidly  connected  to  the  outer  cas¬ 
ing.  The  permanent  magnet  of  the  exciter 
which  serves  as  a  reaction  mass  is  flexibly 
attached  to  the  casing  by  thin  leaf  springs.  It 
is  rated  at  1-pound  driving  force  which  is  suffi¬ 
cient  for  testing  many  structures. 


CALIBRATION 

Sensitivity 

The  sensitivities  of  the  force  pickup  and 
accelerometer  parts  of  the  impedance  head 
were  performed  on  an  electro^namic  vibration 
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Fig.  Z  -  Setup  used  to  calibrate  the  force  pickup  and  accelerometer 
in  the  impedance  head 


exciter  using  the  test  setup  shown  in  Fig.  2. 

The  sensitivity  of  the  force  pickup  was  deter¬ 
mined  by  measuring  its  output  with  two  differ¬ 
ent  weights  attached  to  the  impedance  head. 

The  ratio  of  the  outputs  from  the  force  pickup 
to  standard  accelerometer  mounted  on  top  of 
the  impedance  head  was  measured.  The  meas¬ 
urement  was  repeated  with  a  2-pound  weight 
inserted  between  the  standard  pickup  and  ta¬ 
pered  adapter  block  on  top  of  the  impedance 
head.  The  force  sensitivity  was  the  change  in 
force  pickup  output  divided  by  2  pounds  and  the 
applied  acceleration.  From  this  force  sensitiv¬ 
ity  determination  performed  at  one  frequency, 
the  total  weight  acting  on  the  force  pickup  in¬ 
cluding  the  effective  portion  of  the  mounting 
bolt  through  the  head  and  the  effective  end 
weight  of  the  head  was  determined.  Knowing 
the  total  weight  the  frequency  response  was  de¬ 
termined  by  measuring  the  force  pickup  to 
standard  accelerometer  outputs  from  10  cps  to 
5000  cps  without  the  2-pound  weight  attached. 
This  part  of  the  calibration  was  repeated, 
measuring  the  ratio  of  the  output  of  the  accel¬ 
erometer  in  the  head  to  the  standard  acceler¬ 
ometer  output,  to  determine  the  accelerometer 
sensitivity  over  the  same  frequency  range.  The 
2-pound  weight  was  not  used  at  the  higher  fre¬ 
quencies  in  order  to  avoid  relative  motion 


between  the  standard  accelerometer  and  imped¬ 
ance  head.  A  voltage  ratio  circuit  [ll]  was  used 
to  reduce  calibration  errors. 

The  results  of  the  calibration  are  shown  in 
Fig.  3.  The  sensitivity  of  the  accelerometer  in 
the  impedance  head  is  nearly  constant  through¬ 
out  the  frequency  range,  increasing  slightly 
near  5000  cps.  Also,  the  sensitivities  of  the 
force  pickup  are  nearly  constant  in  the  same 
frequency  range.  The  sensitivity  for  use  with 
an  aluminum  1/2-inch-diameter  bolt  is  5.8 
mv/lb.  The  sensitivity  is  reduced  to  5.1  mv/lb 
when  a  1/2-inch-diameter  steel  bolt  is  used. 
This  change  in  sensitivity  occurs  since  a  small 
portion  of  the  total  force  is  applied  to  the  bolt 
while  the  rest  of  the  force  is  applied  to  the 
head.  Of  course,  the  force  sensitivity  is  the 
voltage  output  divided  by  the  total  force.  The 
portion  of  the  total  force  applied  to  the  bolt  is 
the  ratio  of  the  bolt  stiffness  to  head  stiffness. 
The  impedance  head  stiffness  computed  from 
the  modulus  of  elasticity  and  geometric  dimen¬ 
sions  of  the  various  parts  is  2.7  x  10^  Ib/in. 
Therefore,  the  change  in  sensitivity  resulting 
from  using  different  size  bolts  can  be  computed. 
The  computed  change  in  sensitivity  for  the  steel 
and  aluminum  1/2-inch-diameter  bolts  is  ap¬ 
proximately  12  percent  which  agrees  closely 
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Fig.  3  -  Calibration  results  on  an  impedance  head.  The  force  pickup  and  acceler¬ 
ometer  were  calibrated  by  comparison  witii  a  standard  accelerometer.  Tlie 
phase-angle  measurements  were  made  with  a  dual-beam  oscilloscope. 


with  the  results  shown  in  Fig.  3.  When  no  bolt 
is  used,  the  increase  in  sensitivity  is  7  percent 
compared  to  the  value  obtained  with  the  alumi¬ 
num  bolt. 

The  phase  angle  between  the  force  and  ac¬ 
celeration  outputs  is  zero  degrees  throughout 
the  frequency  range.  This  indicates  that  the 
acceleration  motion  and  force  applied  to  a 
structure  are  faithfully  reproduced  in  the  output 
signals  from  the  head  and  no  correctiem  to  the 
measured  phase  angle  of  the  impedance  need  be 
made. 

The  effective  end  weight  of  the  impedance 
head  was  measured  by  repeating  the  atove 
described  force  calibration  with  several  differ¬ 
ent  weights  attached.  The  output  from  the  force 
piclnip  corresponding  to  1-g  acceleration  was 
measured  for  each  weight  applied  externally  to 
the  impedance  head.  These  data  were  plotted 
and  the  points  fitted  with  a  least  squares  line. 
The  results  for  two  impedance  heads  are  shown 
in  Fig.  4.  The  intercepts  of  the  lines  with  the 
abscissa  indicates  the  effective  end  weight  of 
the  impedance  head  is  about  0.3  pound. 


Amplitude  Linearity. 

The  amplitude  linearities  of  the  impedance 
head  were  verified  on  electrodynamic  vibration 
exciters  equipped  with  resonant  beams  and  on  a 
shoclc  motion  calibrator. 

Calibrations  up  to  112  pounds  and  100  g 
were  made  on  vibration  exciters.  With  a  small 
weight  on  top  of  the  impedance  head,  sufficiently 
high  accelerations  were  achieved  by  mounting 
it  at  the  center  of  a  free- free  beam.  In  order 
to  achieve  sufficiently  high  forces,  it  was  nec¬ 
essary  to  mount  the  liead  on  the  end  of  the  beam. 
A  standard  accelerometer  was  mounted  on  top 
of  the  head  and  the  calibration  was  made  using 
the  same  method  used  for  the  above  sensitivity 
calibrations. 

Accelerations  up  to  502  g  were  applied  on 
shoclc  motion  calibrators  similar  to  that  de¬ 
scribed  in  Ref.  12.  The  impedance  head  and 
standard  accelerometer  were  attached  to  an 
anvil.  Two- inch-diameter  and  4-inch-diameter 
hardened  steel  balls  were  dropped  on  the  anvil. 
Rubber  paddings  were  used  on  the  impact 
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Applied  WtigM,  Pounds 

Fig.  4  -  Normalized  force  output  measured  with  different 
applied  weights.  The  impedance  head  apparent  end  weight  is 
indicated  by  the  abscissa  intercept  (0.3  pound). 


surface  to  obtain  the  desired  accelerations.  The 
amplitude  linearity  was  determined  by  measur¬ 
ing  the  ratio  of  the  standard  acceierometer  to 
impedance  head  accelerometer  outputs  which 
were  photographed  on  the  screen  of  an  oscillo¬ 
scope.  The  results  are  given  in  Table  1  and 
Fig.  5.  The  deviations  from  amplitude  linearity 
are  less  than  the  measurement  errors.  The 
Model  2225  standard  accelerometer  (Fig.  5) 
was  previously  calibrated  throughout  its  useful 
acceleration  range  by  the  absolute  method  de¬ 
scribed  in  Ref.  [12]. 


TABLE  1 

Amplitude  Linearity  of  Endevco 
Impedance  Head  Model  2110 


Force  Gage 

Accelerometer 

Applied 

Force 

(lb) 

Deviation 

(%) 

Applied 

Acceleration 

(g) 

Deviation 

(%) 

0.03 

+  1.0 

0.03 

0.0 

0.1 

+0.5 

0.1 

0.0 

1 

0.0 

1 

0.0 

17 

-0.5 

10 

0.0 

28 

0.0 

50 

0.0 

56 

0.0 

100 

+0.4 

70 

0.0 

95 

0.0 

84 

0.0 

192 

+0.1 

112 

+0.5 

412 

+0.9 

502 

0.0 

Resonance  Frequency 

Resonance  frequency  measurements  were 
also  made  on  the  shock  machine.  Less  padding 
was  used  on  the  anvil  to  reduce  the  pulse  dura¬ 
tion  in  order  to  produce  resonance  frequency 
excitation  in  the  impedance  head.  The  results 
are  shown  in  Fig.  6.  The  oscillogram,  Fig.  6A, 
indicates  the  resonance  frequency  of  the  accel¬ 
erometer  part  of  the  impedance  head  is  about 
28,000  cps.  The  theoretical  response  for  an 
ideal  pickup,  Eq.  10  in  Ref.  13,  with  a  28,000- 
cps  resonance  frequency  is  in  close  agreement 
with  the  slight  rise  near  5000  cps  noted  in  the 
sinusoidal  calibration  (Fig.  3)  of  the  impedance 
head.  The  beating  present  in  Fig.  6A  apparently 
is  due  to  the  individual  accelerometers  in  the 
head  having  slightly  different  resonance  fre¬ 
quencies.  The  accelerometer  part  of  the  im¬ 
pedance  head  consists  of  three  individual  accel¬ 
erometers  connected  in  parallel  with  a  single 
output. 

The  re^X)nse  was  measured  also  with  the 
force  pickup  in  the  impedance  head  subjected  to 
shock  excitation.  Measurements  were  made  on 
two  different  structures.  The  first  structure 
consisted  of  four  parts  including  the  impedance 
head  for  a  total  length  of  about  4  inches.  The 
second  structure  had  five  parts  for  a  total 
length  of  7  inches.  The  corresponding  re¬ 
sponses  from  the  force  pickup  indicates  a 
structural  resonance  of  20,400  cps  and  10,300 
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Fig.  5  -  Comparison  calibrations  performed  on  a  shock  machine  to 
verify  the  amplitude  linearity  of  the  accelerometer  in  an  impedance 
head.  The  standard  accelerometer  (2225)  and  impedance  head  (2110) 
were  mounted  on  an  anvil  in  a  drop*ball  shock  machine  and  the  ratios 
of  the  outputs  compared. 


cps,  Figs.  6B  and  C,  respectively.  These  fre¬ 
quencies  correspond  to  the  resonance  frequen¬ 
cies  for  the  above  two  structures.  This  indi¬ 
cates  that  the  impedance  head  is  capable  of 
making  measurements  at  very  high  frequencies. 
Considering  the  mass- stiffness  characteristics 
of  the  impedance  head  only,  computations  would 
indicate  its  resonance  frequencies  would  occur 
at  frequencies  much  higher  than  those  present 
in  Figs.  6B  and  C. 


Other  Characteristics 

In  addition  to  the  calibrations  just  pre¬ 
sented,  tests  were  performed  to  establish  the 
other  characteristics  of  the  impedance  head. 

The  force  sensitivity  of  the  accelerometer 
part  of  the  impedance  head  was  determined  on  a 
hydraulic  testing  machine.  The  applied  load 


was  released  instantaneously  with  a  quick-acting 
valve.  An  oscilloscope  was  used  to  measure 
the  output  from  the  accelerometer.  No  signifi¬ 
cant  output  from  the  accelerometer  was  meas¬ 
ured  with  an  applied  load  of  5000  pounds  which 
was  the  capacity  of  the  testing  machine.  This 
indicated  the  response  of  the  accelerometer  is 
of  the  order  of  0.0001  g/lb  or  less  which  cor¬ 
responded  to  the  resolution  obtainable  with  the 
test  equipment.  Since  this  interaction  in  the 
impedance  head  is  so  small,  it  is  capable  of 
measuring  very  large  impedances  accurately. 

The  other  characteristics  including  trans¬ 
verse  sensitivity  of  the  accelerometer,  capaci¬ 
tances,  internal  resistances,  and  the  like,  are 
typical  to  that  of  other  piezoelectric  transduc¬ 
ers.  Both  the  force  pickup  and  accelerometer 
are  insulated  from  the  case  of  the  impedance 
head  in  order  to  eliminate  electrical  noise  in 
auxiliary  electronic  instruments. 
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Fig.  6  -  Resonance  response  of  an 
impedance  head.  (A)  Resonance 
frequency  of  the  accelerometer 
part  of  the  impedance  head  is  Z8,000 
cps.  (B)  and  (C)  Force  response 
with  the  impedance  head  mounted 
in  two  structures  resonate  at  Z0,400 
cps  and  10,300  cps,  respectively. 


contact  area  of  the  block  with  the  beam  is  an 
annulus  with  an  area  of  0.23  inch^.  The  force 
and  acceleration  outputs  from  the  impedance 
head  were  connected  to  a  pair  of  voltage  ampli¬ 
fiers  and  directly  to  voltmeters  and  a  dual 
beam  oscilloscope.  No  filtering  circuits  were 
used.  The  power  amplifier  used  with  the  force 
generator  was  operated  at  sufficiently  low  lev¬ 
els  so  no  distortion  was  present,  except  through 
very  narrow  frequency  ranges  at  resonance 
where  the  force  or  acceleration  tends  to  zero. 
The  magnitude  of  the  impedance  was  computed 
from  the  measured  voltage  outputs  times  the 
force  pickup  and  accelerometer  sensitivities 
determined  from  the  above  calibrations.  No 
corrections  were  made  to  the  data  for  the  ap¬ 
parent  weight  of  the  impedance  head,  bolt,  or 
tapered  adapter  block.  The  phase  angle  of  the 
apparent  weight  was  measured  directly  on  a 
dual-beam  oscilloscope  and  by  Lissajous 
figures. 

The  results  are  given  in  Fig.  8.  The  ex¬ 
perimentally.  measured  resonances  and  the  the¬ 
oretically  computed  resonances  for  the  beam 
are  given  in  Table  2.  The  cantilever  beam  res¬ 
onance  at  78  cps  was  characterized  by  a  single 
nodal  point  at  the  center  of  the  beam  and  large 


TABLE  2 

Resonance  Frequencies  of  a  Beam 
Measured  with  an  Impedance  Head 


Measured 

Resonance 

Frequency 

(cps) 

Theoretical 

Resonance 

Frequency 

(cps) 

- - - 

Normal  Mode 

Number 

Type 

78 

74 

1st 

Fixed- free 

121 

117 

1st 

Free-free 

482 

460 

2nd 

Fixed-free 

619 

631 

3rd 

Free-free 

1340 

1289 

3rd 

Fixed- free 

1513 

1560 

5th 

Free-free 

2560 

2526 

4th 

Fixed- free 

2790 

2900 

7th 

Free- free 

4080 

4176 

5th 

Fixed- free 

4400 

4652 

9th 

Free -free 

IMPED.\NCE  OF  A  BEAM 

In  order  to  illustrate  typical  ranges  of  im¬ 
pedance  and  frequency,  measurements  were 
made  on  an  uniform  l^am.  The  test  setup  is 
illustrated  in  Fig.  7.  The  force  generator  and 
impedance  head  were  connected  rigidly  to  a 
3/4-  by  3-  by  36-inch  aluminum  beam  (2024-T4) 
with  a  1/2- inch-diameter  steel  bolt.  The  bolt 
passed  through  a  steel  tapered  adapter  block 
provided  with  a  1/2-inch  clearance  hole.  The 


di  placements  at  both  beam  ends.  The  free- 
free  beam  frequencies  excited  were  odd  num¬ 
bered  modes  characterized  by  an  antinodal  point 
at  the  center  of  the  beam.  It  was  difficult  to 
excite  the  even-mode  frequencies  because  of 
the  effect  of  the  bolt  and  adapter  block  when  the 
center  of  the  beam  tends  to  form  an  inflection 
point.  It  is  interesting  to  note  that  an  attempt 
to  correct  for  the  apparent  mass  of  the  imped¬ 
ance  head,  bolt,  and  adapter  block  would  not 
alter  the  results  at  or  near  the  antiresonances 
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Fig.  7  -  Mechanical  impedance  head  and  force  generator  connected 
to  a  free-free  beam.  The  beam  is  suspended  from  overhead  by 
rubbe  r bands. 


(regions  of  high  impedances).  Also,  an  attempt 
to  correct  for  this  apparent  mass  at  resonances 
(regions  of  low  impedances)  would  produce  only 
small  differences  between  the  measured  results 
and  theoretical  resonances.  In  the  case  of  this 
beam,  good  results  were  obtained  without  cor¬ 
recting  for  the  impedance  head  and  fixture  ap¬ 
parent  weight. 


CONCLUSIONS 

Although  several  different  methods  have 
been  used  to  measure  mechanical  impedance, 
the  direct  measurement  method  with  force  and 
motion  transducers  is  becoming  attractive. 


Impedance  heads  are  available  for  making 
measurements  over  wide-frequency  and  imped¬ 
ance  ranges.  Accurate  impedance  measure¬ 
ments  are  easily  made  with  small  heads  that 
are  sufficiently  stiff.  Measurements  made  on 
an  aluminum  beam  indicate  good  accuracy  is 
maintained  without  correcting  for  the  effective 
end- mass  of  the  impedance  head. 
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Fig.  8  -  Mechanical  impedance  of  an  aluminum  beam  measured  with  an  impedance  head: 
(A)  The  magnitude  of  the  impedance,  apparent  weight,  or  apparent  stiffness  is  indicated 
by  the  appropriate  scales.  (B)  This  plot  is  the  phase  angle  of  the  apparent  weight. 


26 


Occibqls 


Soils,"  Journal  of  Acoustical  Society  of 
America,  Vol.  28,  No.  1,  January  1956, 
pp.  73-76. 

[10]  E.  L.  R.  Corliss  and  W.  Koidan,  "Mechani¬ 
cal  Impedance  of  the  Forehead  and  Mas¬ 
toid,"  Journal  of  Acoustical  Society  of 
America,  Vol.  27,  No.  6,  November  1955, 
pp.  1164-1172. 

[11]  R.  R.  Bouche,  "High  Frequency  Response  and 
Transient  Motion  Performance  Character¬ 
istics  of  Piezoelectric  Accelerometers," 


Proc.  Instrument  Society  of  America  (to  be 
published). 

[12]  R.  R.  Bouche,  "The  Absolute  Calibration  of 
Pickups  on  a  Drop-Ball  Shock  Machine  of 
the  B^listic  Type,"  Proc.  Institute  of  En¬ 
vironmental  Sciences,  1961,  pp.  115-121. 

[13]  "American  Standard  Methods  for  the  Cali¬ 
bration  of  Shock  and  Vibration  Pickups," 
S2.2- 1959,  American  Standards  Associa¬ 
tion,  Inc.,  10  East  40th  Street,  New  York  16, 
New  York. 


DISCUSSION 


Mr.  Murfin  (Sandia  Corp.):  In  your  table  of 
resonances  (Table  2),  all  those  with  even  modes 
under  free-free,  as  you  say,  would  be  difficult 
to  excite  since  they  have  an  inflection  point  at 
the  point  of  support.  If  those  are  replaced  by 
the  higher  modes  of  fixed-free  beams,  the 
agreement  between  the  theoretical  frequencies 
and  the  measured  frequencies  is  really  aston¬ 
ishing.  I  get  the  first  eight  to  match  admost 
perfectly  and  the  last  two  not  quite  so  good. 

Mr.  Bouche:  I  see  you  had  the  time  to  do  a 
little  bit  more  work  than  I  did. 

Mr.  Murfin:  It  is  also  interesting  that  in 
the  fixed-free,  there  is  rigid  body  motion,  hence 
you  would  expect  to  get  a  high  impedance  which 
matches  up  quite  nicely. 

Mr.  Bouche;  I  take  it  then  that  you  have 
continued  the  theoretical  analysis  and  computed 
the  resonances  for  the  four  that  I  did  not  indi¬ 
cate  in  the  paper.  Is  that  correct? 

Mr.  Murfin:  Yes,  I  have. 

NOTE:  As  a  result  of  Mr.  Murfin's  comments, 
additional  computations  have  been  made 
by  Mr.  Bouche.  The  results  are  in¬ 
cluded  in  Table  2  and  in  the  text. 

Mr.  Trummel  (Jet  Propulsion  Laboratory): 

I  wonder  if  you  could  comment  on  the  effect  of 
the  stiffness  of  the  gage  on  frequency  limita¬ 
tions? 

Mr.  Bouche:  I  can  comment  on  it  to  a  cer¬ 
tain  extent.  As  I  indicated,  the  stiffness  of  the 
total  gage  is  2.7  x  lO’  pounds  per  inch.  But,  if 
one  considers  only  the  portion  of  the  impedance 
head  between  the  top  of  the  force  transducers 
and  the  point  where  the  impedance  head  is  at¬ 
tached  to  the  structure,  we  have  a  much  higher 
stiffness  because  there  is  only  a  thin  steel  plate 


in  this  region.  Its  stiffness  is  about  a  factor  of 
10  greater,  about  3  x  10^  pounds  per  inch. 

Now,  from  the  practical  point  of  view  my  ex¬ 
perimental  results  on  the  beam  that  I  measured 
showed  stiffnesses  which  were  quite  large. 

This  may  demonstrate  an  approximate  figure 
for  the  maximum  stiffness  that  could  be  meas¬ 
ured.  In  that  case  I  obtained  readings  which 
were  in  excess  of  10’  pounds  per  inch,  nearly 
10^  pounds  per  inch. 

Mr.  Trummel:  Well,  in  some  of  our  expe¬ 
rience  in  measuring  impedance  we  found  that 
we  couldn't  quite  use  that  large  an  area  for  at¬ 
tachment.  In  some  cases  you  have  to  use  an 
adapter  which  would  effectively  give  you  a  much 
less  stiff  attaching  spring  constant.  You  have 
never  talked  about  the  effect  of  that  on  how  far 
out  in  the  frequency  range  you  can  go.  If  you 
use  that  large  an  area  you  stiffen  your  structure 
much  more  than  you  should. 

Mr.  Bouche:  This  would  reduce  the  range, 
however,  I  think,  in  most  practical  structures 
that  occur  in  the  field,  in  most  cases  we  don't 
reach  stifbiesses  quite  this  large. 

Voice:  What  are  the  low-frequency  meas¬ 
uring  limits  of  your  impedance  head  ?  How 
far  is  it  flat? 

Mr.  Bouche:  I  performed  by  calibrations 
down  to  10  cps,  but  I  expect  the  impedance  head 
itself  will  be  flat  to  much  lower  frequencies, 
say  of  the  order  of  2  cps.  The  primary  limita¬ 
tion  here  is  the  electronic  equipment  that  is 
used  with  the  impedance  head.  Piezoelectric 
transducers  inherently  have  a  relatively  large 
internal  electrical  impedance  and  the  impedance 
is  a  function  of  its  capacity.  So,  at  very  low 
frequencies,  one  must  use  amplifiers  or  cathode 
followers  that  have  high  impedances.  If  a  high 
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enoui^h  input  impedance  is  maintained  in  the 
cathode  follower,  there  is  no  reason  to  ejqject 
that  the  transducer  will  not  be  flat  down  to  2 


cycles.  It  is  a  little  bit  more  difficult  to  cali¬ 
brate  at  2  cycles  than  at  10  cycles. 


28 


APPLICATIONS  OF  IMPEDANCE  INFORMATION 


Ralph  E.  Blake 

Lockheed  Missiles  and  Space  Company 
Sunnyvale,  California 


Only  a  few  applications  of  impedance  information  liave  yet  been  made 
to  the  solution  of  engineering  problems.  However,  engineering  research 
is  being  carried  on  by  several  groups  to  develop  methods,  apparatus, 
and  theorems  for  application  to  some  important  problems  in  shock  and 
vibration.  The  greatest  current  effort  is  on  developing  more  effective 
sound  and  vibration  isolation  systems.  Improvements  in  methods  of 
making  field  measurements  of  dynamic  environment  and  methods  of 
simulation  in  the  testing  laboratory  are  also  being  studied.  Ultimately, 
any  area  of  shock  and  vibration  work  which  deals  in  complicated  linear 
systems  should  benefit  from  applications  of  the  techniques  and  knowl¬ 
edge  being  developed. 


INTRODUCTION 

This  paper  is  concerned  with  practical,  en¬ 
gineering  application  of  impedance  techniques 
and  theories  such  as  those  discussed  in  the  two 
previous  papers.  Altliough  applications  to  engi¬ 
neering  problems  are  just  beginning,  engineer¬ 
ing  research  is  making  rapid  progress  toward 
further  uses,  particularly  in  some  rather  im¬ 
portant  areas  of  shock  and  vibration.  In  general, 
we  might  hoi)e  to  apply  impedance  ideas  to  any 
problem  inwiiich  we  know  the  dynamic-response 
characteristics  of  some  parts  and  need  to  know 
the  response  when  the  parts  are  joined  together. 
Complexity  of  the  details  is  not  a  problem  so 
long  as  the  parts  are  linear. 

The  engineers  working  on  structural  im¬ 
pedance  have  had  to  "write  the  book"  as  they 
went  along,  developing  equipment  and  techniques 
before  running  the  experiments  that  they  really 
wanted.  Progress  can  be  expected  to  become 
more  evident  as  tliey  build  up  from  these  foun¬ 
dations.  One  of  the  most  important  a^ects  of 
this  progress  is  tiie  least  predictable.  As 
masses  of  data  are  gathered  and  studied  and  as 
more  cases  are  worked  out  theoretically,  we 
can  hope  that  new  generalizations,  or  laws,  will 
be  developed.  Such  general  laws  of  behavior  of 
complex  mechanical  systems  may  someday  en¬ 
able  us  to  tackle  problems  which  look  too  com¬ 
plicated  and  unwieldy  today. 

In  .addition  to  describing  some  of  the  cur¬ 
rent  .oiiplications.  I  will  discuss  the  present 


work  toward  future  applications.  Then  I  wish 
to  illustrate  some  generalizations  that  have  be¬ 
gun  to  show  up  in  the  results  already  reported. 

This  paper  does  not  pretend  to  survey  all 
of  the  activities  in  the  field.  First,  I  am  sure 
that  I  have  not  seen  all  the  pertinent  reports. 
But  those  that  I  have  seen  have  too  much  mate¬ 
rial  to  permit  me  in  the  time  allowed  to  do 
more  than  select  a  few  examples.  I  asked  Dr. 
Belsheim  to  send  me  reports  from  his  files  to 
help  me  brush  up  on  recent  work.  But  I  was 
staggered  to  get  a  carton-full;  it  was  a  stack, 

6  inches  high,  with  a  grand  total  of  over  1000 
pages.  If  structural  impedance  work  is  still  in 
its  infancy,  it's  a  big  baby. 

THE  ISOLATION  MOUNTING  PROBLEM 

Most  of  the  work  on  structural  impedance 
thus  far  has  been  aimed  at  reducing  the  machin¬ 
ery  noise  emitted  by  a  Navy  ship  or  submarine 
which  might  be  picked  up  by  an  enemy  hydro¬ 
phone.  The  hull  and  its  sound  radiation  char¬ 
acteristics  can  be  studied  only  in  expensive 
field  tests.  Machinery  noise  generation  and 
transmission  can  be  measured  in  engineering 
laboratories.  Impedances  of  machinery  and 
foundations  can  be  measured.  Then,  the  effect 
of  various  possible  flexible  mounts  and  founda¬ 
tions  can  be  studied  theoretically  and  exjjeri- 
mentally.  The  goal  is  to  find  design  features 
and  techniques  to  minimize  the  transmitted 
sound.  All  of  these  steps  are  being  actively 
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legend: 

1.  ELECTROMAGNETIC  SHAKER  i.  IMPEDANCE  HEAD 

2.  BED  PLATE  6.  MASSIVE  BLOCK 

3.  rubber  isolation  mounts  7  MACHINE 
4  TYPICAL  WEIGHT  LOCATIONS 


Fig.  1  -  Machine-bed  plate  isolation  system  foundation 


Fig.  Z  -  Test  barge  "Dezi-Bellc''  in  air 


carried  out;  the  following  few  enmples  may  in¬ 
dicate  the  scope  and  nature  of  the  work. 

Figure  1  shows  an  experimental  setup  [l] 
at  the  Electric  Boat  Division.  Impedance 
measurements  are  being  made  on  a  special 
foundation  structure  while  it  supports  a  piece 
of  naval  machinery  on  some  rubter  mounts.  In 
this  test,  the  foundation  is  on  a  massive  block. 


Figure  2  is  a  view  of  a  barge  [2]  at  the 
Electric  Boat  Company  which  is  structurally 
very  much  like  a  section  of  submarine  hull. 

By  adding  foundations  inside,  the  sound  trans¬ 
mission  path  from  foundation  to  hull  to  water 
can  be  studied.  Figure  3  shows  a  point  imped¬ 
ance  measurement  [2]  that  was  made  on  the 
foundation  in  the  barge.  The  complexity  of  the 
structure  can  be  inferred  from  the  number  of 
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Measured  point  impedance,  foundation  G  in  barge  in  water,  X  direction 


Fig.  4  -  Experimental  setup  for  external  excitation  tests 


resonant  peaks  on  the  record.  I  would  guess 
that  over  a  hundred  similar  impedance  meas¬ 
urements  have  been  made. 

Figure  4  shows  the  apparatus  used  at  the 
Arthur  D.  Little  Company  [3]  for  some  imped¬ 
ance  measurements  on  the  mounting  feet  of  an 
electric  motor.  One  of  the  transfer  impedances 
measured  between  one  foot  and  another  is  shown 
on  Fig.  5.  Now,  the  Reciprocity  Theorem  says 
that  the  transfer  impedance  will  be  the  same  if 
one  interchanges  the  force- input  and  velocity- 
output  points.  But  the  interchange  in  this  ex¬ 
periment  gave  the  result  shown  on  Fig.  6. 
Perhaps  the  unsatisfactory  comparison  between 
these  curves  was  due  to  nonlinearity  of  the 
structure.  But  the  experimenter  hesitated  to 
state  this  as  a  conclusion  without  more  direct 
evidence.  This  action  is  in  keeping  with  the 
need  to  proceed  cautiously  and  check  everything 
in  this  unfamiliar  game  until  we  learn  to  avoid 
the  pitfalls.  We  expect  most  structures  to  be 
linear  but  should  not  assume  it  for  any  struc¬ 
ture. 

Theoretical  and  laboratory  studies  of  vari¬ 
ous  resilient  mounts  are  being  carried  on  by 
several  investigators  [4,  5,6).  The  dynamic 
modulus  and  damping  of  many  rubber-like  mate¬ 
rials  have  been  measured  up  to  several  thousand 
cps.  These  have  then  been  used  in  theoretical 
analyses  and  digital  computer  studies  to  show 
the  transmission  with  various  types  of  machine 
iind  foundation  impedances. 

So  far,  the  conclusions  that  have  been  drawn 
are  concerned  with  parts  of  the  over-all  sound 
tr.ansmission  path.  For  e>ample,  it  has  been 
reported  that  wave  transmission  effects  become 
important  in  rubber  mounts  at  high  frequencies 


[4],  that  foundation  resonances  seriously  in¬ 
crease  sound  radiated  into  the  water  [l],  and 
the  like.  It  is  hoped  that  such  generalizations 
will  grow  in  number  and  scope  until  engineering 
design  improvements  can  be  made  in  the  hard¬ 
ware.  Of  course  much  that  is  learned  will  apply 
to  any  isolator  problem. 


SHOCK  AND  VIBRATION  MEASUREMENT 
AND  TESTING 

Most  of  us  at  this  symposium  dealt  with 
shock  and  vibration  as  an  environment  which 
might  damage  some  structure  or  equipment. 

We  go  out  in  the  field  and  measure  acceleration¬ 
time  histories  which  we  hope  are  typical.  Then 
we  set  up  laboratory  test  machines  to  duplicate 
the  worst  features  of  the  field  environment. 
These  machines  are  then  used  to  test  the  struc¬ 
ture.  It  is  usually  true  that  the  dynamic  sys¬ 
tem  on  which  the  measurement  was  made  is 
somewhat  different  from  the  one  on  which  the 
equipment  will  be  used.  Most  often  the  equip¬ 
ment  itself  is  new  and  will  supplant  the  one 
used  during  the  field  trial. 

Now,  Norton's  Theorem  [7)  shows  that  a 
change  in  the  equipment  impedance  will  change 
the  input  motion,  so  that  our  measurement  does 
not  tell  us  what  the  test  machine  should  do.  It 
is  common  practice  to  assume  that  the  meas¬ 
urement  is  "typical"  and  that  several  measure¬ 
ments  will  somehow  envelope  the  one  we  need. 

It  is  true  that  a  small  equipment  will  not  influ¬ 
ence  the  motion  of  a  heavy  foundation  at  low 
frequencies;  but  this  assumption,  in  practice, 
becomes  increasingly  untenable  as  frequencies 
get  higher  and  supporting  structures  get 
trimmed  to  minimum  weights.  The  best  defense 
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Transfer  mobility  and  the  phase  versus  frequency  fig.  6  -  Transfer  mobility  and  the  phase  versus  frequency 

(G.  E.  3-hp  motor)  (G.  E.  induction  motor) 


of  the  assumption  is  that  there  is  no  real  alter¬ 
native  available,  and  it  has  not  led  to  disaster, 
so  far  as  we  know. 

Although  we  may  have  no  alternative  today, 
we  can  hardly  be  comfortable  with  the  knowl¬ 
edge  that  our  methods  do  not  have  a  firm  basis. 
This  problem  of  what  we  should  measure  in  the 
field,  how  we  should  present  data,  and  what  the 
testing  machine  should  duplicate  has  interested 
me  for  many  years. 

Our  first  real  evidence  that  something  was 
wrong  came  from  comparing  the  actual  strength 
of  some  heavy  naval  equipment  with  that  re¬ 
quired  by  an  envelope  through  the  peaks  of  shock 
spectra  measured  at  the  equipment  mounting 
points  [8].  Although  the  equipment  easily  sur¬ 
vived  the  shock,  our  calculations  called  for  con¬ 
siderable  strengthening.  Further  study  showed 
that  the  heavy  equipment  had  produced  valleys 
(or  "dips")  in  the  shock  spectra  at  the  equipment 
design  frequencies.  In  this  instance  we  had 
overestimated  the  design  acceleration  by  a  fac¬ 
tor  of  five.  The  overdesign  is  much  less  for 
equipment  that  is  lighter  relative  to  its  founda¬ 
tion  [9, 10|. 

It  would  seem  that  a  quantitative  e;q>res- 
sion  of  what  we  mean  by  "lighter  relative  to  its 


foundation"  would  require  measurements  of  im¬ 
pedance  of  equipment  and  foundation  through  the 
frequency  range  of  interest.  Testing  would  be 
more  true-to-life  if  the  output  impedance  of  the 
test  machine  duplicated  the  important  charac¬ 
teristics  of  a  foundation  [ll]. 

So,  our  present  position  is  that  we  can  see 
that  what  we  are  doing  is  inadequate.  Theory 
tells  us  what  to  do,  but  this  would  not  be  practi¬ 
cal  in  most  cases  [7].  What  is  needed  is  more 
experience  and  data  to  show  us  what  parts  of 
the  theoretically  perfect  approach  can  be  safely 
neglected.  In  short,  further  research  is  needed 
to  seek  out  practical  approximations. 

AN  APPUCATION  TO  TESTING 

Dr.  Belsheim  and  his  colleagues  at  U.  S. 
Naval  Research  Laboratory  (USNRL)  have  made 
an  interesting  application  of  impedance  infor¬ 
mation  to  the  testing  of  some  rocket  payloads 
[12].  The  only  measurement  of  the  vibration  of 
the  rocket  motor  to  which  the  satellite  was  to  be 
attached  was  made  during  a  static  firing.  Then 
the  impedance  of  the  motor  attachment  point 
was  also  measured  with  the  setup  shown  in  Fig. 
7.  This  impedance  was  recorded  as  an  "appar¬ 
ent  weight"  (force  divided  by  acceleration)  as 


shown  in  Fig.  8.  The  apparent  weight  of  the 
static-firing  fixture  was  measured,  as  well  as 
that  of  the  payload  (Fig.  9).  From  this  group  of 
measurements,  it  became  possible  to  estimate 
the  input  vibration  which  would  have  occurred 
if  the  payload  had  been  attached  in  the  place  of 
the  test  stand. 


AN  APPLICATION  TO  VIBRATION 
MEASUREMENT 

I  have  strongly  implied  that  vibration 
measurements  should  be  combined  with  a 
measurement  of  the  impedance,  but  we  are  not 
yet  able  to  make  use  of  these  on  a  routine  ba¬ 
sis.  At  Lockheed,  we  have  been  using  a  very 
simple  type  of  impedance  measurement  to 
check  on  the  installation  of  vibration  pickups  in 
the  Polaris  missile.  A  simple  spring  scale  is 
used  to  apply  a  known  force  to  the  pickup  loca¬ 
tion  through  a  piece  of  fine  music  wire,  as 
shown  in  Fig.  10.  Cutting  the  wire  causes  a 
sudden  step  change  in  the  force.  The  resulting 
transient  oscillations,  such  as  Fig.  11,  are  re¬ 
corded  with  an  oscilloscope  and  camera. 

Suppose  a  10-pound  step  force  excited  a 
10-g  oscillation  at  1500  cps,  as  it  did  in  Fig. 

11.  We  conclude  that  a  1500-cps  resonance 
involves  the  motion  of  only  a  pound  or  so  of 
structure.  Hence  a  ISOO-cps  peak  in  a  flight 
record  is  due  to  a  "local  resonance"  and  is  not 
representative  of  the  hundreds  of  pounds  of 
structure  that  we  are  interested  in.  Although 
the  lightweight  and  high  frequency  of  modern 
pickups  have  largely  done  away  with  resonance 
of  the  pickup  or  its  bracket,  resonance  of  the 
local  structure  can  produce  just  as  misleading 
a  signal. 

With  this  simple  step  force  system,  we 
have  been  able  to  show  whether  a  pickup  mount¬ 
ing  is  stiff  enough;  we  have  also  been  able  to 
edit  some  high  peaks  out  of  our  flight  vibration 
spectra.  If  we  had  needed  it,  the  transient  re¬ 
sponse  could  have  been  analyzed  into  an  imped¬ 
ance  curve. 


DEVELOPMENT  OF  GENERALIZATIONS 

In  this  brief  and  incomplete  survey,  I  have 
tried  to  show  the  scope  and  promise  of  work 
toward  application  of  impedance  information. 

In  order  to  be  able  to  think  about  such  complex 
structures  and  intricate  impedance  data,  it  is 
clear  that  we  must  discover  some  simplifying 
principles.  In  the  reports  already  published 
there  are  many  examples  of  summing  up  the 
important  trends  noted  in  the  data. 


In  a  Westinghouse  report  [5],  we  read  "The 
terminating  impedance  of  the  heavier  structures 
usually  has  a  resistive  component  about  equal 
to  the  reactive  at  frequencies  above  about  200 
cps." 

An  Electric  Boat  Division  report  says, 
"Damping  and  the  use  of  concentrated  masses 
are  effective  in  reducing  the  transmission  effi¬ 
ciency  of  foundations." 

Every  report  has  several  conclusions  of 
this  general  type.  It  is  this  sort  of  boiling  down 
of  the  data,  together  with  some  additional 
sweeping  generalizations,  which  will  enable  us 
to  solve  problems  which  seem  too  complex  now. 

Although  the  laws  which  we  hope  to  find  are 
largely  unpredictable,  we  can  search  for  them 
by  trying  to  ask  the  right  questions.  I  will 
close  with  a  discussion  of  some  studies  which 
we  have  made  at  Lockheed  to  find  possibly  use¬ 
ful  trends. 

In  the  August  1960  session  [ll]  of  this  sym¬ 
posium,  a  theoretical  study  of  a  very  simple 
three-mass  model  of  a  vehicle-equipment  sys¬ 
tem  was  presented.  Results  showed  how  the 
response  of  the  "equipment"  varied  with  natural 
frequency,  mass,  and  damping.  The  behavior 
of  this  simple  system  suggested  some  interest¬ 
ing  effects  which  should  be  watched  for  in  real 
complex  systems. 

More  recently,  I  have  been  interested  in 
the  implication  of  Norton's  Theorem  that  a 
foundation  vibration  can  become  much  more 
severe  if  the  equipment  impedance  is  suitably 
changed.  There  is  a  value  of  equipment  imp^- 
ance  which  will  make  the  resonant  frequency  of 
the  total  system  coincide  with  any  given  forcing 
frequency.  Such  a  resonant  condition  defines 
the  worst  that  could  happen  to  an  equipment  at 
that  frequency.  It  would  be  of  some  interest  to 
see  how  a  vibration  measurement  which  one 
might  obtain  in  the  field  compares  with  this 
worst  possible  case. 

A  computer  program  which  had  been  used 
at  Lockheed  for  studying  a  50-degree  of-freedom 
representation  of  a  missile  was  readily  adapted 
to  this  new  investigation.  The  vibration  re¬ 
sponse  at  an  output  point  near  one  end  of  the 
missile  was  calculated  for  a  unit  sinusoidal 
force  near  the  other  end.  This  is  the  transfer 
admittance  plotted  on  Fig.  12  for  no  equipment 
at  the  output  point.  Then,  it  was  assumed  that 
an  undamped  equipment  was  fastened  to  the 
output  point  of  the  missile  (which  was  damped 
to  a  Q  of  25).  The  reactive  admittance  of  the 
equipment  was  assumed  to  be  equal  and  opposite 
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I'ig.  12  -  Theoretical  transfer  admittances  through  a 
missile  structure 


to  that  of  the  output  point  of  the  missile.  This 
produced  a  resonant  response  of  the  system 
with  amplitude  limited  only  by  the  damping  of 
the  missile.  The  foundation  response  on  Fig. 

12  for  this  case  of  "Equipment  T\ined  to  Reso¬ 
nance"  shows  an  increase  by  a  factor  which 
often  exceeds  10. 

The  acceleration  input  to  the  equipment  is 
only  half  the  story;  the  force  input  is  also  of 
interest.  As  a  basis  for  comparison,  we  first 
computed  the  force  which  would  act  against  an 
equipment  which  was  an  infinite  rigid  mass. 

The  ratio  of  this  "blocked  force"  to  the  input 
force  is  shown  in  Fig.  13.  The  force  for  an 
"Equipment  Tuned  to  Resonance"  is  also  plotted. 
The  latter  exceeds  the  blocked  force  by  the 
same  ratio  as  that  for  the  admittance  curves 


of  Fig.  12.  The  resonant  output  force  was  often 
many  times  larger  than  the  input  force. 

Curves  of  the  type  shown  on  Figs.  12  and 
13  could  be  plotted  from  any  measurement  of 
vibration  of  a  structure,  if  one  also  measured 
the  resistive  and  the  total  admittances  of  the 
pickup  location.  The  curves  suggest  how  un¬ 
conservative  a  given  measurement  could  be  in 
the  worst  circumstances. 

This  study  by  itself  should  not  be  re¬ 
garded  as  conclusive.  It  was  cited  as  an  il¬ 
lustration  of  the  sort  of  probing  and  question¬ 
ing  which  we  should  be  doing  to  establish  a 
more  reliable  basis  for  shock  and  vibration 
practices. 
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DISCUSSION 


Mr.  Gertel  (Allied  Research):  I  just  want 
to  probe  a  little  bit  on  the  data  that  you  pre¬ 
sented  on  the  experiments  perform^  by  the 
Arthur  D.  Little  Company,  showing  transfer 
impedances  which  didn't  match  going  forwards 
and  backwards.  Is  it  possible  that  if,  perhaps, 
rotational  impedance,  say  torques  and  rotational 
accelerations,  were  measured  rather  than  the 
forces  and  translational  accelerations,  there 
might  have  been  a  better  correlation?  Was  that 
the  problem  ? 

Mj.  Blake:  It  seems  possible  to  me.  I 
haven't  worked  that  out.  But  it  is  conceivable 
that  there  was  some  sort  of  rotational  constraint 
offered  by  the  impedance  head,  although  I  am 
not  familiar  with  the  apparatus. 


Mr.  Gertel:  What  I  was  thinking  is  that  the 
picture  showed  the  item  of  equipment  mounted 
from  bungee  cords.  Probably  the  elastic  center 
was  fairly  central  to  the  equipment,  yet  the 
forces  were  applied  on  a  mounting  foot  which  wa.s 
eccentric  to  the  elastic  center  of  this  bungee 
system.  This  suggests  that  the  rotational  mo¬ 
tions  were  more  strongly  involved,  or  were  at 
any  rate  present,  in  addition  to  the  translational 
observations. 

Mr.  Blake:  I  don't  think  that  I  can  add  to 
your  comment.  I  suspect  that  the  flexible  cord 
didn't  play  any  part  at  this  frequency.  Theo¬ 
retically,  this  theorem  is  supposed  to  be  good 
regardless  of  whether  there  is  translational 
motion  present. 
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Mr.  Gertel:  It  was  just  a  thought.  Now  I 
had  one  other  experience  I  wanted  to  relay  about 
this  go  slow  on  impedance.  Like  Bob  Plunkett, 
I'm  in  favor  of  impedance  and  I'm  all  for  women 
and  motherhood.  We  recently  had  a  very  strange 
experience.  We  had  some  damped  elastomer  on 
which  we  were  making  some  impedance  meaS' 
urements.  In  plotting  the  forces  and  accelera¬ 
tions  on  impedance  paper,  we  came  out  with  a 
beautiful  spring-like  characteristic  indicating 
no  damping  present  whatsoever  in  this  particu¬ 
lar  elastomer.  As  an  additional  check,  we  ran 
a  conventional  transmissibility  test,  where  a 
mass  was  mounted  on  top  of  this  elastomer  and 
the  input- output  relationships  were  plotted.  The 
mass  involv^  was  about  20  pounds  and  the  fre¬ 
quency  was  about  60  cycles.  We  wound  up  with 
a  transmissibility  of  4-1/2,  indicating  that  there 
was  a  considerable  amount  of  damping  present 
in  this  elastomer.  The  impedance  equipment  or 
the  impedance-measuring  experiment  would 
have  led  us  to  believe  that  there  was  absolutely 
no  damping  in  this  particular  item.  It  gave  us 
a  very  weird,  strange,  misleading  kind  of  re¬ 
sult  as  the  additional  esqieriment  showed.  I  just 
mention  that  to  indicate  that  there  are  some 
subtle  things,  which  I  don't  understand,  which 
happened  in  this  particular  e^qieriment.  There 
are  lots  of  subtle  things  happening  in  impedance 
work  which  some  of  us  would  like  to  learn  more 
about. 

Mr.  Nutt  (Whitworth  Gloster  Aircraft  Ltd., 
England):  In  describing  the  experiment  to  de¬ 
termine  the  transmissibility  from  one  point  to 
another,  you  remarked  that  the  possible  ejqila- 
nation  of  the  discrepancy  was  nonlinearity  in 
the  system  and  then  went  on  to  say  that  this  was 
unusual.  This  is  not  our  experience.  We  find 
that  practicaily  every  structure,  well  in  fact,  if 
one  wants  to  be  a  purist,  every  structure  is 
nonlinear  up  to  a  point,  but  that  most  missile 
structures  are  nonlinear  to  a  significant  degree. 
Would  you  agree  with  this? 

Mr.  Blake:  Yes,  I'm  afraid  so.  I'm  opti¬ 
mistic  about  structures  being  linear  but  one 
experience  of  mine  is  that  perhaps  they  are 
not.  Our  attemiits  to  measure  Polaris  vibration 
in  flight  show  a  very  b.ad  scatter  of  successive 


measurements  in  successive  flights  for  pickups 
at  the  same  point.  One  possible  explanation  is 
that  these  missiles,  off  the  same  production 
line,  are  not  in  fact  dynamically  the  same. 

Maybe  there  are  subtle  differences  in  the 
joints— I  wish  I  knew.  But  that  the  missile  is 
not  linear  is  our  best  contender  for  an  expla¬ 
nation. 

Voice:  About  one  year  ago,  we  ran  some 
impedance  tests  on  an  old  vehicle.  We  were 
interested  in  the  apparent  mass,  and  I  was 
really  surprised  that  after  we  went  through  the 
first  natural  frequency,  our  apparent  mass  in¬ 
stead  of  190  pounds  with  the  fixture  was  only 
about  20  pounds.  It  really  amazed  me.  The 
other  thing  I  would  like  to  ask  you  is  if  you  have 
say  a  spent  rocket  case  with  some  of  the  vis- 
colastic  material  in  there,  would  your  curves 
on  the  impedance  plot  be  shifted  to  the  right  or 
left  as  you  would  have  in  a  nonlinear  type  vibra¬ 
tion  or  would  there  be  no  effect  at  all  on  this? 

Mr.  Blake:  You  are  asking  the  influence  of 
some  propellant  inside  ? 

Voice:  That's  right,  because  you  have  the 
viscolastic  material  acting  as  a  nonlinear-type 
damping  system. 

Mr.  Blake:  I  don't  think  I  can  answer  your 
question. 

Voice:  Perhaps  Dr.  Belsheim  may  be  able 
to  answer  that. 

Mr.  Blake:  I  believe  his  test  was  with  a 
motor  empty  and  full. 

Dr.  Belsheim:  I  had  one  empty  and  one  full 
and  there  was  quite  a  difference  in  the  apparent 
weight  looking  into  the  motor  in  each  of  these 
two  cases.  Considerably  more  damping  showed 
up  apparently  in  the  loaded  motor,  but  this  is 
rather  inconclusive  and  I  wouldn't  want  to  say 
much  more. 

Mr.  Blake:  Surely  if  you  were  above  what 
I  will  call  the  first  natural  frequency  of  the  pro¬ 
pellant,  it  could  go  either  way. 


♦ 
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APPLICATION  OF  MECHANICAL  ADMITTANCE  DATA 
TO  THE  SOLUTION  OF  A  PRACTICAL  PROBLEM 


Roy  W.  Muslain 

Norlronics  Division,  Northrop  Corporation 
Hawthorne,  California 


This  paper  describes  the  practical  application  of  mechanical  admit¬ 
tance  (reciprocal  of  impedance)  data  to  the  numerical  solution  of  a  vi¬ 
bration  isolation  problem.  The  dynamic  solution  follows  the  methodol¬ 
ogy  of  the  design  procedure  which  was  formulated  by  the  SAE  G-5 
Aerospace  Shock  and  Vibration  Committee  under  the  chairmanship  of 
Dr.  C.  T.  Molloy.  The  solution  results  in  the  selection  of  the  most 
favorable  vibration  isolator  to  protect  rigid  equipment  from  dynamic 
e.'icitalions  of  a  noinTgid  supporting  structure.  The  equipment  mass, 
the  fragility  of  the  equipment,  the  free  vibration  of  the  supporting  struc¬ 
ture,  and  tlio  mechanical  admittance  of  the  supporting  structure  are 
processed  by  techniques  given  in  the  G-5  Document,  "Design  of  Vibra¬ 
tion  Isolation  Systems."  The  solution  is  presented  in  steps  and  pro¬ 
duces  optimti  values  of  stiffness  and  damping  for  the  selected  vibration 
isolator.  An  envelope  of  transmissibility  curves  for  all  permitted  sys¬ 
tem  designs  is  established. 


INTRODUCTION 

The  concepts  and  tundanientals  ol  mechani¬ 
cal  impedance  which  are  essential  to  an  under¬ 
standing  ol  this  paper  have  been  discussed  and 
described  adequately  in  the  previous  presenta¬ 
tions  ol  this  session.  The  primary  goal  ol  the 
learned  SAE  G-5  Aerospace  Shock  and  Vibra¬ 
tion  Committee  has  been  the  preparation  of  the 
document,  "Design  of  Vibration  Isolation  Sys¬ 
tems."  The  intent  of  this  worthwhile  G-5  Docu¬ 
ment  is  to  provide  the  design  engineer  with  the 
methodology  necessary  to  select  the  optimum 
vibration  isolation  design  through  the  applica¬ 
tion  of  mechanical  impedance  techniques. 

Whereas  the  formulation  of  the  design 
procedure  and  its  supporting  mathematical 
equations  will  be  covered  in  the  G-5  Docu¬ 
ment,  the  presentation  of  the  procedure  proper 
and  the  mathematics  (especially  derivations) 
will  be  minimized  in  this  paper.  However, 
formulations  and  procedural  explanations  will 
be  given  wherever  needed  for  clarity  of 
presentation. 

This  numerical  solution  will  serve  as  a 
guide  for  engineers  desiring  to  apply  the  design 


procedure  of  the  document.  This  solution  pre¬ 
sents  a  practical  application  which  determines 
whether  vibration  isolators  are  necessary  to 
protect  delicate  airlwrne  equipment.  Then,  if 
isolators  are  required,  the  solution  provides 
data  needed  to  select  a  satisfactory  isolation 
system. 


THE  SAMPLE  PROBLEM 

A  specific  source  receiver  configuration 
has  been  selected  for  treatment  in  this  numeri¬ 
cal  illustration.  This  practical  application  as¬ 
sumes  (1)  that  the  supporting  structure  is  a 
nonrigid  mass  which  is  dynamically  linear  with 
resonant  vibrational  characteristics,  (2)  that 
the  equipment  behaves  as  rigid  mass,  (3)  that 
the  vibration  isolators  arc  linear  massless 
spring-dashpot  elements,  (4)  that  there  is 
mainly  translational  vibrations  in  each  of 
three  mutually  perpendicular  directions,  and 
(5)  that  the  susceptibility  to  failure  of  the  equip¬ 
ment  can  be  described  by  either  random  or 
sinusoidal  fragility  curves.  The  procedure  is 
limited  to  the  selection  of  center-of-gravity 
isolation  systems  [ij. 


43 


Required  Data 

The  weight  of  the  equipment,  the  fragility 
levels  of  the  equipment,  the  free  vibration  of 
the  supporting  structure,  and  the  mechanical 
admittance  of  the  supporting  structure  are  re¬ 
quired  to  adequately  accomplish  the  solution  of 
the  problem.  The  weight  of  the  equipment  (W) 
in  this  example  is  100  pounds.  The  fragility  of 
the  equipment  iu,|  is  given  in  the  accompany¬ 
ing  graph  (Fig.  1).  Levels  lUfl  are  di.splayed 
which  can  be  expected  to  cause  equipment  fail¬ 
ures.  This  graph  also  shows  the  free  vibration 
of  the  supporting  structure  |U,1  ,  i.e.,  the  vi¬ 
bration  amplitudes  which  exist  at  the  support 
points  without  the  load  of  the  equipment.  The 
mechanical  admittance  of  the  supporting  struc¬ 
ture  is  given  in  Fig.  2.  The  mechanical  admit¬ 
tance  of  the  source  is  the  ratio  of  the  complex 
amplitude  of  velocity  to  the  complex  amplitude 
of  sinusoidal  force  producing  that  velocity,  both 


quantities  being  measured  at  the  structural 
support  of  the  equipment  |2].  The  mechanical 
admittance  of  the  source,  ,  is  equal  to  A,  i  jn^ 
where  and  are  the  real  and  imaginary 
parts  of  the  admittance,  respectively.  The 
reciprocal  of  admittance  is  the  mechanical  im¬ 
pedance  which  is  equal  to  Rs+i^s  where 
and  jXj  are  the  real  and  imaginary  parts  of 
the  mechanical  impedance,  respectively.  At 
first  the  G-5  procedure  was  based  on  the  utili¬ 
zation  of  mechanical  impedance.  However,  it 
has  been  found  advantageous  to  use  the  mechani¬ 
cal  admittance  because  the  isolator- source  com¬ 
bination  produced  a  modified  source  whose  ad¬ 
mittance  is  the  sum  of  the  admittance  of  the 
source  and  of  the  isolator  [3]. 

The  success  of  future  dynamic  progress 
depends  on  the  utilization  of  mechanical  im¬ 
pedance  heads  and  force  generators  with  im¬ 
proved  recording  media  which  provide  accurate 


Frequency  (cps) 

Fie.  1  -  Plots  of  frafiilily  amplitude  and  f  ree  -  sun  rt  e  .im  pi  it  nrl  e 
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Fig.  2  -  Mechanical  admittance  of  supporting  structure 


measurements  of  force,  motion,  and  phase  rela¬ 
tionships.  The  ideal  instrumentation  system 
will  supply  automatic  plots  of  the  real.  A,,  and 
the  imaginary,  B, ,  parts  of  the  admittance  as  a 
function  of  frequency.  Presently,  no  equipment 
has  been  devised  to  automatically  plot  the  data 
as  presented  in  Fig.  2.  However,  some  systems 
are  being  developed  that  are  very  close  to  that 
stage.  One  system  [4]  being  developed  will  au¬ 
tomatically  record  ^e  following  data  as  a  func¬ 
tion  of  frequency:  (1)  the  ratio  of  the  accelera¬ 
tion  amplitude  to  the  force  amplitude,  and  (2)  the 
phase  angle  that  the  acceleration  lags  the  force. 
The  determination  of  the  admittance  amplitude, 
|Y,  i ,  can  be  accomplished  by  automatically 
plotting  the  ratio  of  (1)  on  special  charts  which 
show  both  the  admittance  amplitude,  |Y,1,  and 
the  ratio  of  (1)  as  a  function  of  frequency.  It 


then  becomes  a  simple  calculation  to  use  the 
admittance  amplitude,  |Y,1  ,  and  the  phase  angle, 
0 ,  to  determine  the  real.  A, ,  and  imaginary, 

B, ,  parts  of  the  admittance  [4]: 


A,(f)  = 

|Y,(f)|  sin  0 

(1) 

B,(f)  = 

|Y,(f)|  cos  0 

(2) 

Determination  of  the  Need  for  Isolation 

The  first  step  in  the  procedure  is  the  cal¬ 
culation  and  the  plotting  of  the  dimensionless 
admittances  A,  and  6,  as  a  function  of  fre¬ 
quency.  This  is  accomplished  by  using  the 
value  of  *  (100  pounds)  along  with  the  real  (A,) 
and  an  imaginary  ( i  B,)  values  of  the  original 
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data  given  in  Fig.  2.  These  data  are  inserted 
into  the  following  equation  from  the  procedure 
to  obtain  the  dimensionless  admittance  A,  and 

B.: 


The  value  of  2r;/g  is  1.62  x  10' 


A,(f) 


2v  Wf 
g 


A.(f) 


(3) 


B,(f) 


2v  Wf 
g 


B,(f)  . 


(4) 


The  values  of  A,(f)  and  B,{f)  for  this 
problem  are  shown  in  the  log- log  graph  of  Fig. 
3.  Distinction  between  the  positive  and  nega¬ 
tive  values  of  B,  is  made  by  plotting  the  posi¬ 
tive  values  as  a  solid  curve  and  the  negative 
values  as  a  dashed  curve.  Since  A  is  always 


positive,  the  A,  curve  is  displayed  as  a 
solid  line. 


The  next  step  in  the  procedure  is  the  cal¬ 
culation  and  plotting  of  the  maximum  allowable 
transfer  function,  H^.  The  G-5  procedure  de¬ 
termines  a  criterion  of  maximum  allowable 
transfer  function,  H„  ,  as  an  upper  limit  for  the 
vibration  of  the  equipment.  This  transfer  func¬ 
tion  is  given  as  [2] 


Hf(f)  I 

u.(f)  I 


for  the  sinusoidal  case,  and  is  given  as 


(5) 


(6) 


1  10  100  1000 

F rcqucncy  (cps ) 

Fii;.  ^  -  Normalized  admitlance  of  .supporlinL;  slruclure 
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for  the  random  case.  Where 

u,  is  the  complex  amplitude  of  sinus¬ 
oidal  input  displacement  expressing 
the  fragility  limit  of  the  equipment, 
inches. 

u,  is  the  complex  amplitude  of  sinus¬ 
oidal  output  displacement  of  the  un¬ 
loaded  source  vibration,  inches. 

S(Uf)  is  the  spectral  density  di^lacement 
which  the  equipment  can  withstand 
as  a  function  of  frequency,  in.^/cps. 

S(u,)  is  the  spectral  density  displacement 
of  the  unloaded  source,  in.*/cps. 


The  fragility  of  the  equipment  and  the  free 
vibration  of  the  source  may  be  expressed  also 
in  terms  of  velocity  or  acceleration  provided 
that  the  same  units  are  used  consistently  in  the 
above  equations.  The  data  from  Fig.  1  are  used 
to  compute  for  this  sample  numerical  solu¬ 
tion.  The  calculated  values  of  for  each  fre¬ 
quency  are  displayed  in  Fig.  4.  Cross-hatching 
lines  have  been  placed  above  the  curve  to 
show  that  the  curve  is  an  upper  limit.  The  lower 
frequency  limit  of  is  designated  f '  and  the 
upper  frequency  limit  is  designated  f“  .  The 
frequency  range  f'  to  f“  is  the  range  within 
which  all  succeeding  transmissibility  or  trans¬ 
fer  function  curves  are  to  be  plotted.  In  this 
problem  f'  is  5  cps  and  f"  is  1000  cps. 


Fig.  1  -  Plots  of  H  transfer  function 
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Another  transfer  function,  ,  is  calculated 
and  added  to  Fig.  4.  is  defined  to  be  the  ab¬ 
solute  value  of  the  ratio  of  the  velocity  of  the 
load  mass  to  the  free  velocity  of  the  source  when 
the  mass  is  attached  directly  to  the  source; 


Applying  the  relationships  of  Eqs.  (3)  and  (4), 
the  transfer  function  ll^f. becomes  [2]: 

=  [l  -  (10) 


Since 


0,(0,) 


^  1 

0,(oj)  1  +  .io,mY,(o,) 


(7) 


(8) 


ll„(tj)  can  be  equated  to: 

+  j^mo,A,(o,)j  j  .  (9) 


The  H„  curve  is  found  by  using  data  in  Fig.  3 
and  by  filling  in  Tables  1  and  2.  First,  the  peaks 
of  are  determined  by  inspection  of  Fig.  3. 

The  peaks  and  valleys  of  are  located  at  or 
near  the  frequencies  where  B,(  f )  equals  one. 
Therefore,  these  frequencies  occur  at  the  in¬ 
tersection  of  the  B,  positive  curve  with  the 
ordinate  value  of  1.0.  The  points  for  this  prob¬ 
lem  are  shown  in  Fig.  3.  The  lowest  frequency 
intersection  (29  cps)  is  labeled  as  and  the 
next  intersection  (48  cps)  is  labeled  fo^-  The 
other  intersections  need  not  be  labeled.  Since 
A,  is  a  relatively  small  value  at  foi  ,  the  fpi 


TABLE  1 

Determination  of  Frequencies  for  Transfer  Function  Calculations 


f'=^cps  f  =  1000  CPS 


Multiples  of  Center  Column 

PEAKS  and  VALLEYS 

Multiples  of  Center  Column 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.4 

2 

6 

50 

53 

9 

15 

20 

26 

foi  =  29 

32 

34 

43 

fo2  =  48 

53 

67 

196 

288 

1000 

5 

9 

12 

15 

17 

19 

24 

H 

28 

39 

50 

56 

62 

78 

112 

336 

_ _ 

1000 

TABLE  2 
Calculation  of 


f 

(cps) 

>■  © 

@ 

@ 

i^-®! 

® 

® 

. 

5 

0 

0.0200 

0.980 

0.960 

0 

0.960 

1.02 

9 

0 

0.066 

0.934 

0.872 

0 

0.872 

1.07 

15 

0.0122 

0.187 

0.813 

0.660 

0.000 

0.660 

1.23 

19 

0.0290 

0.315 

0.685 

0.469 

0.001 

0.470 

1.46 

24 

0.068 

0.560 

0.440 

0.194 

0.005 

0.199 

2.24 

r  29 

0.155 

_ 

0 

0 

6.43 

32 

0.280 

1.33 

0.33 

0.109 

0.078 

0.187 

2.31 

34 

0.440 

1.67 

0.67 

0.449 

0.194 

0.643 

1.25 

43 

2.50 

3.80 

2.80 

7.84 

6.25 

14.09 

0.266 

Jr  48 

9.0 

— 

0 

0.111 

53 

8.0 

-4.7 

5.7 

32.5 

64 

96.5 

0.102 

;  67 

4.4 

-3. >95 

4.95 

24.5 

19.4 

43.9 

0.151 

196 

0.39 

-2.4 

3.4 

11.6 

0.15 

11.8 

0.2911 

288 

0.068 

-2. 

3. 

9. 

0.00 

9.00 

0.33 

,  1000 

0.0200 

-2. 

3. 

9. 

0.00 

9.00 

0.33 
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intersection  represents  a  peak  in  the  curve 
and  since  is  relatively  high  at  f 02.  the 
intersection  represents  a  valley  [1].  The  meth¬ 
odology  of  the  G-5  procedure  is  followed  to  de¬ 
termine  the  minimum  number  of  points  neces¬ 
sary  to  properly  define  the  curve.  It  is 
shown  in  the  procedure  that  the  transfer  func¬ 
tion  curve  for  a  single  peak  or  valley  is  satis¬ 
factorily  determined  by  values  at  frequencies 
which  are  the  following  multiples  of  the  fre¬ 
quency  at  which  the  peak  or  valley  occurs:  0.1, 
0.3,  0.5,  0.7,  0.9,  1.0,  1.1,  1.4,  2,  6,  50.  These 
multiples  have  been  computed  for  this  sample 
problem  and  are  listed  in  Table  1.  The  values 
of  ffl,  (29  cps)  and  fgj  (48  cps)  are  entered 
in  the  center  column  of  Table  1,  with  the  small¬ 
est  value  being  listed  first.  Entries  were  made 
to  the  left  of  and  the  lowest  multiple  of  0.1 
was  deleted  and  replaced  by  the  value  of  f' 
which  is  to  equal  to  5  cps.  Then  entries  were 
made  to  the  right  of  f^j  and  the  highest  mul¬ 
tiple  was  deleted  and  replaced  by  f" .  Since 
the  ratio  of  foj/fgj  0.75,  the  fol¬ 

lowing  procedure  [1]  was  used  to  fill  in  the 
points  of  Table  1:  successive  entries  were 
made  —  1.1  to  the  right,  0.9  to  the  left, 
1.4  f(,i  to  the  right,  0.7  (qj  to  the  left.  At  this 
point  the  frequency  multiples  crossed  and  the 
last  entry  on  the  right  (41)  was  deleted. 

The  frequencies  determined  in  Table  1  for 
the  ll„  curve  are  used  in  Table  2  to  compute 
Eq.  (11)  from  the  procedure: 

M„(f)  =  ([1  - 

Table  2  which  is  taken  from  the  procedure  is  a 
convenient  form  for  efficiently  computing  values 
of  .  The  frequencies  found  in  Table  1  are 
entered  into  column  1  of  Table  2  with  the  lowest 
frequency  at  the  top  of  the  column.  The  check 
murks  to  the  left  of  column  1  indicate  frequen¬ 
cies  which  represent  actual  intersections.  For 
these  two  frequencies  (29  and  48  cps)  a  zero  has 
been  entered  in  column  4  and  a  horizontal  line 
has  been  entered  in  column  3  since  it  is  not 
needed.  Values  of  A,  and  B,  taken  from  Fig. 

3  are  listed  in  columns  2  and  3,  respectively. 
Negative  signs  have  been  included  wherever 
applicable  for  values  of  B,  .  Zeroes  have  been 
entered  for  values  of  A,  and  B,  which  are  less 
than  0.01.  The  computations  indicated  in  col¬ 
umns  4  through  8  are  straightforward  and  yield 
the  curve.  The  values  in  column  8  are  equal 
to  . 

The  data  for  this  example  are  displayed 
in  Fig.  4  for  comparison  with  the  curve. 

The  pe.ak  of  the  11^  curve  is  labeled  (29 
cps)  and  the  valley  is  labeled  fo2  (48  cps). 


Whereas  the  peak,  fgj ,  of  the  curve  exceeds 
the  upper  bound  of  and  exceeds  above 
150  cps,  vibration  isolation  is  required  and  this 
sample  problem  is  continued  to  determine  a  sat¬ 
isfactory  isolation  system. 


Determination  of  a  Minimum  Natural 
Frequency  for  the  Isolation  System 

The  definition  of  an  adequate  isolation  sys¬ 
tem  requires  the  determination  of  a  range  of 
acceptable  natural  frequencies.  The  lower 
limit  of  the  frequencies  corresponds  to  a  mini¬ 
mum  stiffness  for  the  isolators  and  is  a  function 
of  (1)  the  maximum  permitted  displacement, 
^max,  across  the  isolators  and  (2)  the  available 
sway  space  in  the  particular  installation  [1]. 

The  maximum  permitted  displacement,  of 

the  equipment  relative  to  its  static  position,  for 
this  problem  is  given  as  1  inch  and  the  largest 
peak  acceleration  of  the  supporting  structure, 
'•max  >  given  as  10  g.  This  information  is  used 
in  Eq.  (12)  to  find  [1]: 

fmin  =  =  9.9  cps.  (12) 


Determination  of  the  Region  of  Permissible 
Damping  and  Natural  Frequency  Parame¬ 
ters  for  the  First  Peak  in  the  Transfer 
Function  which  Considers  Isolators  and 
Structural  Admittance 

An  approximation  to  the  transfer  function, 
Hj,  of  the  desirable  isolation  system  is  pro¬ 
vided  by  the  family  of  classical  transmissibil- 
ity  curves  of  Fig.  5.  Figure  5,  in  the  form  of 
an  overlay,  is  placed  on  Fig.  4  of  this  problem 
to  establish  approximate  damping  values  for  the 
first  peak  in  the  transfer  function,  Hj  ,  of  the 
equipment-isolator  combination.  The  ordinate 
of  H  equal  to  one  for  the  overlay  is  juxtaposi- 
tioned  on  the  H  equals  one  ordinate  of  Fig.  4. 
With  this  juxtaposition  maintained,  the  overlay 
is  moved  laterally  on  Fig.  4.  The  frequency 
range  above  fo2  is  ignored  during  this  over¬ 
lay  estimation.  The  overlay  is  adjusted  so  that 
fn,  equals  f^jp.  The  largest  and  smallest 
values  of  in  Fig.  5  (in  the  range  from  0.01  to 
0.5)  for  which  all  of  the  corresponding  curve  in 
Fig.  5  lies  below  in  the  frequency  range 
below  fgj  are  noted.  The  extreme  values  of  /; 
at  the  frequency  are  plotted  on  Fig.  6. 

These  0.5  and  0.05,:-  are  indicated  as  points  1 
and  2,  respectively. 

Additional  usage  of  the  overlay  is  required 
to  further  define  the  approximate  region  of 
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(  I  Frequency  of  'o.  0^5 

I  first  peak 

Fig.  3  -  Family  of  curves  approximating  the  transmissibilily  achieved  at  low  frequen¬ 
cies  when  isolators  arc  placed  between  the  equipment  and  the  supporting  structure 
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TABLE  3 


Limiting  Values  of 

I - , 

i^nlfmax)  =  1/21 


r 


1 


Dnl 


© 

fnl 

(cps) 

® 

® 

1/2 

© 

®-  1 

© 

|/®/2 

Pnifmax) 

® 

© 

.  . 

10 

0.081 

12.4 

11.4 

1.68 

0 

<0.01 

15 

0.185 

5.40 

4.40 

1.05 

0.012 

<0.01 

18 

0.277 

3.61 

2.61 

0.81 

0.024 

0.012 

20 

0.353 

2.74 

1.74 

0.66 

0.035 

0.018 

25 

0.625 

1.60 

0.60 

0.387 

0.077 

0.040 

27 

0.799 

1.25 

0.25 

0.250 

0.130 

0.067 

permissible  damping  for  the  first  peak.  Figure 
5  is  adjusted  again  so  that  r„|  coincides  with 
a  frequency  that  is  5  cps  higher  than  or 
until  a  frequency  is  encountered  wherein  there 
is  a  noticeable  change  in  one  of  the  extreme 
values  of  p,,, ,  whichever  occurs  first.  Once 
more,  the  largest  and  smallest  values  of  are 
noted  and  plotted  on  Fig.  6.  These  values  are 
shown  by  points  3  and  4,  respectively,  on  Fig.  6. 
This  process  is  continued  until  a  frequency  is 
reached  wherein  no  transmissibility  curve,  for 
any  of  the  r  values,  falls  below  or  until 
f„,  =  0.9  fo, ,  whichever  occurs  first. 

Now  tliat  extreme  values  of  for  a  range 
of  f„,  values  have  been  determined  and  plotted 
in  Fig.  6,  the  next  step  is  to  correct  this  region 
by  using  Table  3  to  obtain  limiting  values  of  ,  . 

The  values  must  lie  between  a  maximum 
value  (m.ix)  and  a  minimum  value  /  „,  (min) 
in  order  that  the  damping  be  obtainable  by 
spring-dashpot  isolators.  Equations  (13),  (14), 
and  (15)  which  are  taken  from  the  procedure  are 
used  to  define  the  limits  for  the  maximum  and 
minimum  damping  at  each  distinct  value  of  f^,. 
Table  3  provides  a  convenient  and  efficient  form 
for  computing  Eqs.  (13)  and  (14).  The  maximum 
limit  at  each  f  ,  frequency  is  defined  bv 
Eq.  (13): 

Equation  (14)  gives  the  expression  for  a  di¬ 
mensionless  parameter  when  ■  is  at 
minimum 

2  •  „1  (min) 

D„,  .-.t  .„,imin)  -  j  (14) 

1  (■  4  .  (min) 


where 


D(Pi.  f/fj) 


(15) 


and  is  displayed  graphically  in  Fig.  7  which  is 
taken  from  the  procedure.  Figure  7  presents  a 
family  of  curves  of  D(/5.,  f  fj)  versus  f  f. 
for  constant  values  of  pj .  Thus,  having  estab¬ 
lished  A,(f„,)  values  which  are  equal  to  D„, 
at  (min) ,  it  is  easy  to  obtain  the  /•„!  (min) 
values  graphically  from  Fig.  7.  The  values  of 
f„,  (10,  15,  18,  etc.)  (or  this  problem  are  listed 
in  column  1  of  Table  3.  Corresponding  values 

which  are  taken  from 

Fig.  3  of  thir  problem  are  shown  in  columns  2 
and  6,  respectively,  of  Table  3.  Using  the  input 
data  of  column  2,  the  straightforward  calcula¬ 
tions  of  columns  3,  4,  and  5  are  completed  to 
yield  (m.nx)  values  which  are  equal  to 


of  Eq.  (13). 

The  input  data,  A^(  f„,) ,  of  column  6  are 
equivalent  to  values  of  Eq.  (14)  and  are  used  to 
enter  Fig.  7  to  obtain  the  values  of  (min)  in 
column  7.  The  ■  (min)  values  of  column  7 
are  obtained  from  Fig.  7  by  using  D  =  = 

A,(f„,)  and  f  f;  is  equal  to  1.  Therefore,  the 
ordinate  of  Fig.  7  is  entered  with  the  values  of 
from  column  6  of  Table  3  and  the 
abscissa  is  entered  only  at  f  f  j  equal  to  1. 

The  data  in  columns  5  and  7  are  compared  with 
the  extremes  previously  plotted  on  Fig.  6. 
The  points  which  are  outside  of  the  limits 
defined  by  columns  5  and  7  of  Table  3  must  be 
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Fig.  7  -  Plot  of  D(/3j,  f/fj)  versus  f/f^  for  constant  /3^ 
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changed  to  equal  the  value  at  the  nearest  limit. 
These  changes  are  shown  on  Fig.  6  where  the 
original  /8„,  points  at  25  and  27  cps  have  been 
crossed  out  and  lowered  to  match  the  (max) 
points  of  column  5  at  25  and  27  cps.  These 
newer  points  at  25  and  27  cps  are  0.387  and 
0.35,  respectively. 

The  points  in  Fig.  6  which  were  estab¬ 
lished  in  the  above  steps  are  connected  by 
dashed  lines  in  Fig.  6.  This  area  is  labeled 
the  "Region  of  Permissible  Parameters 
for  the  first  peak  of  H." 


Determination  of  the  Region  of  Permissible 
Damping  and  Natural  Frequency  (fi) 
Parameters  for  the  Equipment-Isolator 
System 

The  isolator/source  parameters  f„,  and 
/3„,  for  the  range  of  low  frequencies  from 
to  fo2  were  established  in  the  previous  section. 
Now  in  this  section,  these  same  parameters, 
f„i  and  ,  are  used  to  obtain  the  required 
isolator  parameters  f  j  and  /3j .  Data  from 
Section  4  are  entered  into  Table  4  to  compute 
Eqs.  (16),  (17),  (18),  and  (19)  which  are  also 
given  in  the  procedure.  First,  values  of  f„, 
and  their  corresponding  values  of  for  this 
sample  exercise  are  taken  from  Fig.  6  and  en¬ 
tered  into  columns  1  and  2  of  Table  4.  Two 
values  of  /3„,  are  entered  in  column  2  for  each 
fn,  frequency.  For  example.  Fig.  6  shows  /3„, 
equal  to  0.5  and  0.05  at  10  cps.  This  is  shown 
in  column  2  of  Table  4  wherein  a  double  entry 
is  made  for  the  /3„,  values  at  10  cps.  Two 
asterisks  are  noted  in  column  2  to  indicate  that 
two  values  of  /S„| ,  0.387  at  25  cps  and  0.25  at 
27  cps,  were  determined  in  Table  3  as  maxi¬ 
mum  values  of  , .  Column  3  of  Table  4  dis¬ 
plays  values  of  D„,  which  are  taken  from  Fig. 

7  of  the  procedure.  The  D„|  values  are  found 
by  entering  Fig.  7  with  f/f,  equal  to  1  and 
with  /3,  equal  to  /i3„,  as  listed  in  column  2  of 
Table  4.  Column  4  contains  values  of  A,(f„,) 
which  are  taken  from  Fig.  3  of  this  problem. 
The  value  of  A,  for  10  cps  is  less  than  0.01 
and,  in  agreement  with  the  procedure,  is  en¬ 
ter^  as  a  zero  in  column  4.  Column  5  is 
equivalent  to  column  3  minus  column  4  and 
thus,  column  5  satisfies  Eq.  (16).  For  column 
5  is  equal  to: 

-  A.(f„,).  (16) 

Column  6  of  Table  4  yields  D^,  values 

which  are  part  of  Eq.  (17): 

(17) 

*  “n  I 


Column  6  values  are  obtained  by  dividing  column 
3  values  by  twice  the  values  given  in  col¬ 
umn  2.  Values  of  i,(f„,)  which  are  taken  from 
Fig.  3  are  shown  in  column  7.  Finally,  the 

values  for  Eq.  (17)  are  established  in 
column  8  by  subtracting  column  7  values  from 
column  6  values.  This  siditraction  is  the  same 
as 


"nl 

24 1 


(17) 


At  this  point  in  this  section,  required  ad¬ 
mittance  parameters,  Aj  and  Bj ,  of  the  iso¬ 
lators  have  been  found  and  are  listed  in  columns 
5  and  8,  respectively.  Table  4  calculations  will 
be  continued  to  determine  fj  and  /3j  by  com¬ 
pleting  Eqs.  (18)  and  (19)  of  the  procedure: 


f.  =1 


/  Bi(fnl) 

-2  -  2 

4  (f„,)  +  B,  ff„,) 


(18) 


jL  ii_ 

2  ■ 


(19) 


The  and  values  needed  under 

the  radical  sign  of  Eq.  (18)  are  computed  and 
listed  in  columns  9  and  10,  respectively.  Then 
these  two  columns  are  added  and  tabulated  in 
column  11.  Following  this,  the  Bj  values  of 
column  8  are  divided  by  the  sums  of  A,*  and 
B|2  which  are  found  in  column  11.  These  quo¬ 
tients  are  placed  in  column  12  and  are  equiv¬ 
alent  to  the  quantities  under  the  radical  sign. 
Next,  the  square  roots  of  these  quantities  are 
found  and  listed  in  column  13.  The  values  in 
column  13  are  multiplied  by  values  taken 
from  column  1.  In  turn,  these  products  are 
listed  in  column  14  and  they  represent  the 
values  of  Eq.  (18)  which  are  equal  to  f  ^ .  The 
last  column,  15,  of  Table  4  is  determined  by 
Inserting  values  from  columns  14,  5,  1,  and  8 
into  Eq.  (19)  as  follows: 


-  1  !±  hiljill 

L  "Li  vpfu""  14)  »  (coluim  5) 
2  X  (column  1)  »  (column  8) 


(19) 


These  computations  complete  Table  4  and 
establish  the  required  parameters,  and  /3^, 
at  low  frequencies  for  the  isolators.  The  data 
from  columns  14  and  15  are  plotted  on  Fig.  6 
of  the  problem.  The  high  values  of  (1.701, 
1.25,  and  fall  outside  the  limits  of  the  figure 
and  are  simply  ignored.  A  horizontal  line  is 
drawn  at  =  0.5,  since  this  value  represents 
an  upper-boundary  limitation  to  practical 
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Additional  points  needed  to  define  the  region  of  permissible 


isolation  system  designs.  The  area  below  this 
practical  level  is  labeled  the  "Region  of  Per¬ 
missible  Parameters  fj.  /Sj  for  the  Isolation 
System." 


Investigation  of  the  Transfer  Function, 

Hj(f) ,  at  the  Higher  Frequencies 

Whereas  the  isolator  parameters  f  j  and 
/^i  found  in  the  above  section  have  ignored  fre¬ 
quencies  above  fgj,  it  becomes  necessary  to 
extend  this  investigation  to  include  the  higher 
frequencies.  In  this  section,  the  possibility  of 
the  isolator  transfer  function  exceeding  H„, ,  the 
maximum  allowable  transfer  function,  will  be 
determined.  The  area  where  trouble  may  be 
encountered  and  where  this  possibility  exists  is 
described  as  the  "suspect"  region.  The  trans¬ 
fer  function  of  an  isolator  system  is  given  in 
the  procedure  as; 

H.(f)  =  ([i  -  i,(f)  - 

+  [A.(f)  +  .  (20) 

An  upper  bound  is  established  for  the  Isolator 
transfer  function  when  the  term  (1-B,(f)  -  Bj(f)) 
is  equal  to  zero.  If  the  term  is  not  zero,  either 
in  a  positive  or  negative  case,  the  transfer  func¬ 
tion  is  smaller.  The  upper  bound  to  the  isolator 
transfer  function  is  given  by  Eq.  (21)  from 
the  procedure: 


Hj(f)  =  - - 

Ai(f)  A,(f) 

^ _ 1 _ 

-ii  0(^5;,  f/fj)  +  A.(f)  (21) 

^  i 

Table  5  has  been  set  up  to  facilitate  cal¬ 
culations  of  the  quantity  D(/v,  f  fj)  as  defined 
in  Eq.  (21): 


- 


The  determination  of  D(/5j,  in  Table  5 

is  accomplished  at  several  frequencies  to  check 
the  isolator  system  design  in  the  range  of  in¬ 
terest  above  .  Three  check  frequencies  are 
worthy  of  investigation  in  the  sample  problem. 
Sudden  reductions  in  occur  on  Fig.  4  at  65, 
150,  and  300  cps.  These  frequencies  are  en¬ 
tered  in  column  1  as  values  of  f^.  Correspond¬ 
ing  values  of  for  each  value  of  are  taken 
from  Fig.  4  also  and  entered  in  column  2  of 
Table  5.  Column  3  represents  the  reciprocals 
which  are  listed  in  column  2.  Values 
of  A,(f^)  are  taken  from  Fig.  3  and  are  tabu¬ 
lated  in  column  4.  The  values  of  column  3  minus 
column  4  or  the  values  of  1/H„,(  f ,.)  -  A,(  f  j)  are 
obtained  and  listed  in  column  5.  An  examination 
of  Fig.  6  shows  that  the  highest  frequency  per¬ 
mitted  by  the  solid  region  is  approximately  50 
cps.  This  value  of  {^  is  made  the  first  entry 


fc  [»i(n 


TABLE  5 

Determination  of  fp  /3.  Combinations  for  which  Satisfactory  Transfer  Function  Valves 

are  Assured  at  High  Frequencies 


(D 

® 

® 

® 

®  ^ 

r  ® 

® 

fc 

WJfc) 

1/® 

A.(f,) 

(D-(D 

f* 

fc  fi 

2 

(DX!) 

Region 
Fig.  7 

Satisfactory 

65 

0.4 

2.5 

1.65 

0.85 

50 

1.30 

0.503 

IV 

none 

40 

1.62 

0.334 

I 

>0.12 

30 

2.16 

0.182 

I 

>0.043 

20 

3.25 

0.084 

I 

>0.013 

10 

6.5 

0.020 

V 

all 

150 

0.2 

5.0 

0.165 

4.84 

50 

3.0 

"ITSST" 

IV 

none 

40 

3.75 

0.343 

n 

0.055  to  0.32 

30 

5.0 

0.194 

n 

0.02  to  0.5 

20 

7.5 

0.086 

V 

all 

35 

4.29 

0.263 

n 

0.024  to  0.42 

32.5 

4.62 

0.226 

n 

0.027  to  0.47 

300 

0.1 

10.0 

0.065 

9.94 

50 

6.0 

0.276 

n 

0.026  to  0.28 

40 

7.5 

0.187 

n 

0.013  to  0.34 

30 

10.0 

0.099 

V 

all 

_ 

L _ 

35 

8.58 

0.135 

m 

<0.43 

56 


in  column  6  of  Table  5.  This  value  (50)  is 
divided  into  the  value  of  (65)  and  the  result 
(1.30)  is  placed  in  column  7.  Column  5  is 
divided  by  the  square  of  the  value  in  column  7 
and  the  quotient  is  listed  in  column  8.  The 
values  in  column  8  are  equal  to  D  (/3. , 

The  first  calculation  of  is 

equal  to  0.503  for  the  condition  wherein  f 
equals  1.30  and  equals  50  cps.  Figure  7 
has  been  divided  into  five  regions.  For  this 
first  case,  Fig.  7  is  entered  at  .  equal  to 
1.30  and  D  equal  to  0.503.  This  falls  into  the 
region  which  is  marked  as  IV  in  Fig.  7.  The 
value  of  IV  is  entered  into  column  9. 


Information  in  column  10  is  found  by  using 
the  regions  of  column  9  according  to  the  follow¬ 
ing  scheme  [1]. 


Column  9 
Region  in  Fig.  7 

I 

n 

in 

IV 

V 


Column  10 
Satisfactory 

> - 

- to - 

< - 

none 

all 


variation  of  ±  20  percent  in  fraction  of  critical 
damping,  and  (4)  that  is  feasible  by  the  current 
state  of  the  art  in  the  production  of  vibration 
isolators.  The  nominal  value  of  the  trial  design 
is  shown  in  Fig.  6.  On  the  basis  of  the  selection 
rules,  the  selected  isolators  have  a  natural  fre¬ 
quency  of  20  ±  2  cps  and  a  fraction  of  critical 
damping  of  0.15  ±  0.03.  It  is  necessary  to  check 
this  trial  design  to  see  that  the  isolator  trans¬ 
fer  function  does  not  exceed  H„ ,  the  maximum 
allowable  transfer  function.  For  this  reason, 
the  isolator  transfer  function,  ,  is  calculated 
in  the  next  section. 


Calculation  of  Hj  the  Transfer  Function 
Achieved  with  the  Isolator  System 
Design 

The  check  of  the  trial  design  is  accom¬ 
plished  by  (1)  finding  the  peaks  and  valleys  of 
Hj ,  (2)  determining  the  frequencies  at  which  to 
calculate  H,  ,  (3)  calculating  the  values  of  the 
normalized  isolator  admittance  at  the  frequen¬ 
cies  found  in  (2)  and  adding  these  to  the  values 
of  normalized  structural  admittance,  (4)  using 
these  data  to  calculate  Hj  ,  and  (5)  plotting  Hj 
on  the  same  graph  as  H„  (Fig.  4)  tor  comparison. 


A  value  of  /3  is  found  on  Fig.  7  for  each  blank 
in  column  10.  Two  values  of  are  determined 
tor  region  H.  One  value  is  read  from  /3  curves 
which  are  traveling  upward  to  the  right,  positive 
slope;  the  other  from  p  curves  which  are  trav¬ 
eling  downward  to  the  right,  negative  slope. 

According  to  the  procedure,  the  region  IV 
provides  no  satisfactory  values  of  /8, .  Where¬ 
upon  the  above  steps  of  Table  5  are  repeated 
using  new  values  of  fj  in  increments  of  10  cps 
in  descending  order  until  the  f^/fj  versus  D 
point  lies  in  region  V  or  until  f  j  is  less  than 
f„i„.  Computations  for  Df/flj.  f^f,)  down  to 
10  cps  are  displayed  in  Table  5  for  this  sample 
solution.  The  data  in  columns  6  and  10  of 
Table  5  are  used  to  correct  the  solid  region  in 
Fig.  6  which  is  suspect  at  high  frequencies. 

The  suspect  region  is  indicated  by  the  slashed 
lines  in  the  high-frequency  area  of  Fig.  6. 


Choice  of  a  Trial  Isolator  Design 

The  next  logical  step  in  the  procedure  is  to 
select  a  resilient  element  having  frequency  and 
damping  parameters  within  the  range  of  the  ac¬ 
ceptable  values  of  Fig.  6.  A  trial  design  is 
selected  within  the  region  of  permissible  pa¬ 
rameters  (1)  that  avoids  the  su^ect  region, 

(2)  that  allows  a  variation  in  natura.'  frequency 
of  at  least  ±  10  percent,  (3)  that  permits  a 


The  isolator  transfer  function,  Hj ,  was 
previously  defined  by  Eq.  (20)  from  the  proce¬ 
dure.  The  data  from  Fig.  1  are  used  in  Tables 
1,6,  and  7  of  the  procedure  to  obtain  the  iso¬ 
lator  transfer  function  curve  for  this  numerical 
solution.  The  peaks  and  valleys  of  this  curve 
are  determined  by  using  the  overlay  of  Fig.  8 
on  the  data  of  Fig.  3.  Figure  8  provides  a  curve 
of  the  dimensionless  admittance  of  a  spring- 
mass  combination  with  a  natural  frequency  f 
The  frequencies  at  or  near  the  peaks  and  val¬ 
leys  of  Hj  are  found  at  the  intersections  of  the 
spring  part  of  the  admittance  of  the  isolator- 
equipment  combination  with  the  imaginary  part 
of  the  structural  admittance.  The  overlay  of 
Fig.  8  is  placed  on  Fig.  3  of  this  problem  with 
the  unity  ordinates  coinciding  and  the  f .  value 
of  Fig.  8  placed  at  the  numerical  value  of  f, 

(20  cps)  on  Fig.  3.  Then,  the  intersections  of 
the  overlay  curve  with  the  B,  curve  are  noted. 
Intersections  occur  near  17  and  56  cps.  These 
values  are  entered  into  Table  1  in  the  column 
labeled  peaks  and  valleys.  The  lowest  fre¬ 
quency  intersection,  17  cps,  indicates  the  pres¬ 
ence  of  a  peak  and  the  56  cps  represents  a 
valley. 

The  G-5  procedure  is  followed  to  deter¬ 
mine  the  minimum  number  of  points  necessary 
to  properly  define  the  H,  curve.  It  is  shown 
that  the  transfer  function  curve  for  a  single 
peak  or  valley  is  satisfactorily  determined  by 
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TABLE  6 

Calculation  of  the  Sum  of  the  Normalized 
Structural  and  Isolator  Admittances 


(D 

13 

0 

0 

0 

f 

f/f. 

D 

2 

(D@^2^i 

n 

D 

0+0 

0+0 

(cps) 

(Fig.  C) 

(2X-0 

Ai 

Bi 

. 

■ 

.. 

A.  +  A, 

B.+  Bi 

5 

0.250 

0.005 

0.060 

0 

0.020 

0.005 

0.080 

9 

0.450 

0.132 

0.027 

0.198 

0 

0.066 

0.027 

0.264 

12 

0.600 

0.175 

0.063 

0.351 

0 

0.118 

0.063 

0.469 

15 

0.750 

0.218 

0.123 

0.544 

0.012 

0.187 

0.135 

0.731 

17 

0.850 

0.242 

0.175 

0.687 

0.019 

0.243 

0.194 

0.930 

19 

0.950 

0.268 

0.241 

0.848 

0.029 

0.315 

0.270 

1.163 

24 

1.20 

0.320 

0.460 

1.28 

0.068 

0.560 

0.528 

1.84 

28 

1.40 

0.358 

0.701 

1.67 

1.28 

0.87 

0.829 

2.54 

39 

1.95 

0.435 

1.655 

2.83 

1.10 

2.65 

2.76 

5.48 

50 

2.50 

0.480 

3.00 

4.00 

9.7 

0 

12.7 

4.00 

56 

2.80 

0.490 

3.97 

4.57 

5.2 

-4.95 

9.17 

-0.38 

62 

3.10 

4.80 

5.18 

2.2 

-4.55 

7.0 

0.63 

78 

3.90 

0.490 

7.45 

6.38 

0.69 

-3.00 

8.14 

3.38 

112 

5.60 

0.430 

13.5 

8.02 

0.28 

-2.22 

13.8 

5.80 

336 

16.8 

0.192 

54.2 

10.75 

0.058 

-2.0 

54.3 

8.75 

1000 

50.0 

0.066 

165 

11.0 

0.02 

-2.0 

165 

9.0 

60 

3.00 

0.370 

3.33 

3.70 

2.90 

-4.8 

6.23 

-1.1 

200 

10.0 

0.300 

30.0 

10.0 

0.107 

-2.0 

30.1 

8.0 

TABLE  7 


Calculation  of  H, 


0 

f 

(cps) 

0 

B,.B, 

® 

l>-®) 

® 

© 

0+0 

H| 

5 

0.080 

0.920 

0.845 

0.845 

1.09 

9 

0.264 

0.736 

0.540 

1.36 

12 

0.469 

0.531 

0.282 

0.286 

1.85 

15 

0.731 

0.269 

0.072 

0.090 

3.33 

17 

0.930 

0.070 

0.005 

0.043 

4.83 

19 

1.163 

0.163 

0.027 

0.099 

3.18 

24 

1.84 

0.84 

0.707 

0.278 

1.01 

28 

2.54 

1.54 

2.37 

0.665 

0.573 

39 

2.76 

5.46 

4.46 

19.9 

7,63 

27.5 

0.190 

50 

12.7 

4.00 

3.00 

9.0 

161 

0.077 

56 

9.17 

-0.38 

1.38 

1.90 

85.9 

0.108 

62 

7.0 

6.3 

0.37 

0.137 

49.0 

49.1 

0.143 

78 

8.14 

3.38 

2.38 

5.65 

66.5 

72.2 

0.118 

112 

13.8 

5.80 

4.80 

23.0 

190 

0.069 

336 

54.3 

8.75 

7.75 

60.0 

0.018 

1000 

165 

9.0 

8.0 

64.0 

VwKil 

27,360 

60 

6.23 

-1.1 

2.1 

4.41 

38.9 

43.3 

0.152 

200 

8.0 

7.0 

49.0 

903 

952 

0.032 
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Normalized  Admittance 


$ 


Fig.  8  -  Plot  of  the  normalized  admittance  of  springmass 
combination  with  natural  frequency  f| 
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values  at  frequencies  which  are  the  following 
multiples  of  the  frequency  at  which  the  peak  or 
valley  occurs:  0.1,  0.3,  0.5,  0.7,  0.0,  1.0,  1.1, 

1.4,  2,  6,  50.  These  multiples  for  Hj  of  the 
trial  design  are  listed  in  Table  1.  The  values 
of  17  cps  and  56  cps  are  entered  in  the  center 
column  of  Table  1,  with  the  smallest  value 
being  listed  first.  Entries  were  made  to  the 
left  of  17  cps  until  f '  was  reached.  Then  en¬ 
tries  were  made  to  the  right  of  56  cps  and  the 
highest  multiple  was  deleted  and  replaced  by 
f*.  Since  the  ratio  of  17/56  is  less  than  0.75, 
the  following  procedure  was  used  to  fill  in  the 
Hj  points  of  Table  1:  successive  entries  were 
made  —  1.1  x  17  to  the  right,  0.9  x  56  to  the 
left,  1.4  X  17  to  the  right,  0.7  x  56  to  the  left, 

2  X  17  to  the  right,  0.5  x  56  to  the  left.  At  this 
point  the  frequency  multiples  crossed  and  the 
last  entry  on  the  right  (34)  was  deleted. 

The  frequencies  determined  in  Table  1  for 
the  H|  curve  are  used  in  Tables  6  and  7  to 
compute  Eq.  (20): 

Hjff)  =  ([l  -  B.(f)  - 

+  ^  (20) 

Table  6  provides  a  convenient  form  for  effi¬ 
ciently  computing  the  sum  of  the  normalized 
structural  and  isolator  admittances.  The  fre¬ 
quencies  for  Hj  found  in  Table  1  are  entered 
into  column  1  of  Table  6  with  the  lowest  fre¬ 
quency  at  the  top  of  the  column.  Using  20  cps 
for  f  j ,  the  values  of  f/fj  are  computed  and 
listed  in  column  2.  The  D(/3j,  f/fj)  values  of 
column  3  are  obtained  by  entering  Fig.  7  with 
the  (/f^  values  of  column  2  and  with  equal 
to  0.15.  The  D(/3,.  f/fj)  values  are  read  from 
the  ordinate  of  Fig.  7.  The  values  of  D  f*/f 
which  are  equal  to  Aj(f)  are  computed  by  mul¬ 
tiplying  the  square  of  column  2  values  by^the 
corresponding  column  3  values.  These  Aj(f) 
values  are  listed  in  column  4  of  Table  6.  Next, 
the  values  of  B|(f)  which  are  equal  to  fD(/3|, 
f/f ,  )/2/’, f ,  are  computed  by  dividing  the  product 
of  columns  2  and  3  by  2/3;.  The  values  of  Bj(f) 
are  tabulated  in  column  5.  Values  of  A,  and  B, 
are  taken  from  Fig.  3  and  are  entered  in  col¬ 
umns  6  and  7,  respectively.  Note  that  the  B, 
values  from  56  to  1000  cps  are  negative.  Col¬ 
umn  8  of  Table  6  is  the  sum  of  the  real  parts  of 
the  normalized  structural  and  isolator  admit¬ 
tances.  Thus,  the  values  in  column  8  are  ob¬ 
tained  by  adding  the  A;  values  in  column  4  to 
the  A,  values  in  column  6.  SimUarly,  column 
9  contains  the  sum  of  the  imaginary  parts  of  the 
normalized  structural  and  isolator  admittances. 
The  values  in  column  9  are  determined  by  adding 


the  B;  values  in  column  5  to  the  B,  values  of 
column  7. 

Now  the  admittance  data  complied  in  Table 
6  are  used  in  Table  7  to  establish  values  of 
Hi(f) .  The  frequencies  used  in  Table  6  are 
entered  into  column  1  of  Table  7.  The  admit¬ 
tance  values  of  (A,  +  Aj)  and  (B,  +  Bj)  are 
taken  from  Table  6  and  are  entered  into  col¬ 
umns  2  and  3,  respectively,  of  Table  7.  The 
absolute  values  of  |l-(B,  +  B;)|  for  this 
sample  exercise  are  computed  and  tabulated  in 
column  4  of  Table  7.  In  turn,  the  square  of 
these  values  is  calculated  and  listed  in  column  5. 
Column  6  contains  the  square  of  the  values  which 
are  listed  in  column  2  of  this  table.  Thus  the 
values  in  column  6  are  equivalent  to  the  square 
of  the  real  part  of  the  normalized  structural 
and  isolator  admittances.  A,  +  A; .  Values  from 
column  5  are  added  to  values  from  column  6  and 
entered  in  column  7.  These  sums  are  equal  to 
that  part  of  Eq-  (20)  which  lies  within  the  outer 
brackets.  Finally,  the  reciprocals  of  the  square 
roots  of  the  values  in  column  7  are  computed 
and  are  entered  in  column  8.  This  completes 
Table  7  for  the  values  in  column  8  are  equal  to 
H; ,  the  isolator  transfer  function.  The  data 
from  column  8  of  Table  7  are  plotted  on  Fig.  4. 
Since  these  data,  H; ,  are  less  than  the  maxi¬ 
mum  allowable  transfer  function,  H„,  at  all  fre¬ 
quencies,  it  is  evident  that  the  trial  design 
(^;  =  0.15  and  f;  =  20  cps)  provides  adequate 
protection  to  the  equipment. 


Statement  of  the  Isolation  System 
Design  Requirements 

It  was  concluded  in  the  previous  section 
that  since  H;  does  not  exceed  H„  at  any  fre¬ 
quency,  the  equipment  is  satisfactorily  isolated 
by  the  trial  design  selection.  Although  it  ap¬ 
pears  that  the  solution  has  been  completed,  it  is 
considered  worthwhile  to  continue  this  study  to 
further  define  the  isolation  system  and  to  es¬ 
tablish  an  envelope  of  permissible  transmis- 
sibility  curves  for  the  isolator-equipment 
combination. 

An  envelope  of  transmissibility  curves  is 
determined  for  all  the  permitted  isolation  sys¬ 
tem  designs  of  Fig.  6.  This  envelope  is  estab¬ 
lished  by  considering  designs  on  the  periphery 
of  the  accepted  design  region  (solid  line)  of  Fig. 
6.  The  overlay  of  Fig.  5,  which  contains  the 
family  of  classical  transmissibility  curves  is 
used  to  find  the  transmissibility  for  each  isola¬ 
tion  system  design  of  /S;  and  fj.  Correspond¬ 
ing  values  of  /3;  and  f ;  are  taken  from  the 
acceptable  design  region  of  Fig.  6.  Then 
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transmissibility  curves  for  these  /5j,  {^  points 
are  obtained  from  Fig.  5,  These  curves  are 
traced  from  Fig.  5  onto  a  three  decade  log-log 
graph.  The  frequency  f  ^  is  indicated  on  the 
abscissa  and  the  transmissibility  values  are 
plotted  on  the  ordinate.  The  envelope  of  trans¬ 
missibility  curves  for  the  permitted  Isolation 
system  design  applicable  to  this  example  is 
shown  in  Fig.  9. 

The  selected  isolators  have  a  natural  fre¬ 
quency  of  20  i:  2  cps  and  a  fraction  of  critical 
damping  of  0.15  ±  0.03.  These  two  ranges  are 
expanded  slightly  by  an  examination  of  Fig.  6. 
Using  the  same  nominal  values,  slightly  greater 
ranges  of  natural  frequency,  f  j ,  and  damping, 
/3.,  are  chosen  from  the  permissible  range  in 
Fig.  6.  A  range  of  20  ±  3  cps  is  selected  for 


the  natural  frequency  and  a  range  of  0.15  ±  0.035 
is  selected  for  the  damping. 

Equation  (22)  provides  an  alternate  method 
of  stating  the  isolator  specifications  [1]: 

g  ‘ 

(22) 

g  ‘ 

The  values  of  the  dynamic  stiffness,  kj ,  and  the 
coefficient  of  viscous  damping,  Cj  ,  in  Eq.  (22) 
are  divided  by  n ,  the  number  of  isolators,  to  ob¬ 
tain  individual  isolator  requirements.  Four 
isolators  are  required  for  this  numerical  exam¬ 
ple.  Individually,  the  isolators  have  a  nominal 
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dynamic  stiffness  of  1023  Ib/in.  for  the  nominal 
frequency  of  20  cps.  Similarly,  the  dynamic 
stiffness  values  for  the  limiting  frequencies  of 
17  and  23  cps  are  found  to  be  739  and  1352  Ib/in., 
respectively.  Thus  the  range  of  stiffness  for  the 
individual  isolators  runs  from  730  to  1352  Ib/in. 
with  a  nominal  value  of  1023  Ib/ln.  at  the  nomi¬ 
nal  frequency  of  20  cps.  The  values  of  the  co¬ 
efficient  of  viscous  damping  for  the  individual 
isolators  are  computed  for  (1)  the  nominal 
and  f>  values,  (2)  the  highest  ySj  and  fj  values, 
and  (3)  the  lowest  and  values.  For  the 
nominal  case  =  0.15;  f  j  =  20  cps),  the  indi¬ 
vidual  coefficient  of  viscous  damping  is  2.44 
lb/in./sec.  Likewise,  the  highest  case  (/3|  = 
0.185;  f  j  s  23  cps)  and  the  lowest  case  (/3^  s 
0.115;  =  17  cps)  values  are  3.46  and  1.59 

lb/in./sec,  respectively.  Obviously,  this  re¬ 
sults  in  a  range  for  the  coefficient  of  viscous 
damping  from  1.59  to  3.46  lb/in./sec  with  a  nom¬ 
inal  value  of  2.44  lb/in./sec  at  the  nominal  jSj 
and  fj  values. 


Rotational  Frequencies 


The  sample  problem  is  extended  to  inves¬ 
tigate  rocking  frequencies  by  basing  calculations 
on  zero  source  mechanical  admittance  [1]  and 
by  assuming  rigid  coupling  between  supports. 
The  rotational  natural  frequencies  for  this  ex¬ 
ample  are  determined  by  using  the  following 
equation  from  the  procedure: 
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Where  f„g,  f^y,  and  f„,  are  the  translational 
natural  frequencies  in  the  x ,  y ,  and  z  direc¬ 
tions,  respectively,  p„,  Py,  and  p,  are  the  radii 
of  gyration  of  the  equipment  about  the  x ,  y ,  and 
z  axes,  respectively;  and  a^,  and  a,  are 
distances  from  the  c.g.  of  the  equipment  to  the 
elastic  centers  of  the  individual  isolators  in  the 
X,  y,  and  z  directions,  respectively  [1]. 


A  satisfactory  isolator  design  has  been  es¬ 
tablished  in  the  previous  sections  and  the  iso¬ 
lator  characteristics  have  been  stated  in  terms 
of  (1)  the  natural  frequency  and  the  fraction  of 
critical  damping,  (2)  a  permissible  transmis- 
slbility  envelope,  and  (3)  the  stiffness  and  the 
equivalent  viscous  damping  coefficient  of  the 
isolation  system.  Since  the  procedure  is  ap¬ 
plicable  to  all  directions  individually,  only  cal¬ 
culations  for  one  of  the  XYZ  axes  are  given  in 
this  presentation.  This  sample  problem  has 
encompassed  the  isolation  of  equipment  which 
is  essentially  symmetrical  about  a  centrally 
located  c.g.  This  sample  problem  is  solved 
with  the  selection  of  four  isolators  which  are 
mounted  symmetrically  about  the  center-of- 
gravity  of  the  isolated  equipment. 

This  solution  has  been  accomplished  while 
assuming  that  only  translational  motion  exists 
in  each  of  the  three  principal  axes.  However, 
this  exclusive  translational  condition  seldom 
occurs  in  actual  practice.  It  is  indeed  a  rare 
design  which  results  in  translational  modes 
occurring  independently  without  rotational 
modes.  Nevertheless,  this  numerical  solution 
has  been  based  on  the  selection  of  the  most 
favorable  design  which  presents  negligible  cou¬ 
pling  and  rocking  vibrations.  Still,  realizing 
that  some  small  amount  of  rotational  motion 
may  be  induced  in  the  equipment,  it  is  consid¬ 
ered  desirable  and  appropriate  td  determine  the 
rotational  natural  frequencies  of  the  selected 
isolation  system. 


To  determine  the  three  rotational  natural 
frequencies  for  this  problem,  the  following  re¬ 
quired  data  are  given; 

f„,  =  20  cps,  f„y  =  20  cps,  and  f„,  =  20  cps. 

a,  2  12  inches.  By  s  12  inches,  a,  =  0, 

Equipnient  dimensions, 

20  X  20  X  20  inches. 


Then 

P*  ~  f’y  ~  ~  ~  *•  inchei 


^20*  X  12*  X  20*  X  12* 

g7I« 


41.59  cpi. 


A  last  final  criterion  from  the  procedure  is 
applied  to  these  rotational  natural  frequencies. 
These  rotational  frequencies  are  considered 
harmful  if  they  fall  within  0.56  to  1.5  of  any 
critical  frequency.  An  examination  of  Fig.  4 
shows  the  Hj  curve  to  have  a  salient  peak  at 
17  cps  which  tends  to  approach  the  H,,  curve. 
Then,  the  critical  range  for  the  rocking  fre¬ 
quencies  runs  from  0.56  x  17  to  1.5  x  17  cps  or 
from  9.5  cps  to  25.5  cps.  Since  the  rotational 
natural  frequencies  found  above  are  outside  of 
this  range,  the  isolators  appear  to  be  located  in 
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transmissibility  curves  for  these  .  f  j  points 
are  obtained  from  Fig.  5.  These  curves  are 
traced  from  Fig.  5  onto  a  three  decade  log-log 
graph.  The  frequency  f  j  is  indicated  on  the 
abscissa  and  the  transmissibility  values  are 
plotted  on  the  ordinate.  The  envelope  of  trans¬ 
missibility  curves  for  the  permitted  isolation 
system  design  applicable  to  this  example  is 
shown  in  Fig.  9. 

The  selected  isolators  have  a  natural  fre¬ 
quency  of  20  ±  2  cps  and  a  fraction  of  critical 
damping  of  0.15  ±  0.03.  These  two  ranges  are 
expanded  slightly  by  an  examination  of  Fig.  6. 
Using  the  same  nominal  values,  slightly  greater 
ranges  of  natural  frequency,  f  j ,  and  damping, 
/Sj,  are  chosen  from  the  permissible  range  in 
Fig.  6.  A  range  of  20  :i:  3  cps  is  selected  for 


the  natural  frequency  and  a  range  of  0.15  i  0.035 
is  selected  for  the  damping. 

Equation  (22)  provides  an  alternate  method 
of  stating  the  Isolator  specifications  [1]; 


*  e  ‘ 
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(22) 


The  values  of  the  dynamic  stiffness,  kj ,  and  the 
coefficient  of  viscous  damping,  Cj ,  in  Eq.  (22) 
are  divided  by  n ,  the  number  of  isolators,  to  ob¬ 
tain  individual  isolator  requirements.  Four 
isolators  are  required  for  this  numerical  exam¬ 
ple.  Individually,  the  isolators  have  a  nominal 
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Fig.  9  -  Envelope  of  maximum  allowable  transmissibilities 
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a  position  that  will  create  negligible  equip¬ 
ment  response  to  the  rocking  motions.  This 
completes  the  numerical  solution  of  the  isola¬ 
tion  design  which  is  adequate  to  protect  the 
airborne  equipment. 


Dr.  S.  Rubin  who  is  responsible  for  the  deriva¬ 
tion  of  the  procedure  and  who  supplied  most  of 
the  mathematics  and  the  sample  problem  used 
in  this  report. 


CONCLUSION 

Mechanical  impedance  and  admittance 
data  are  applicable  to  many  dynamic  problems 
in  various  fields  of  mechanical  engineering. 
However,  this  paper  has  been  confined  to  the 
specific  utilization  of  admittance  concepts  in 
the  reduction  of  vibration  transmission  from 
sources  (the  supporting  structure)  to  receivers 
(the  equipment).  In  this  regard,  this  paper  has 
attempted  to  lead  the  design  en^neer  through 
most  of  the  details  of  a  sample  vibration  isola¬ 
tion  problem  which  was  solved  by  the  procedure 
prepared  by  the  SAE  G-5  Committee  On  Aero¬ 
space  Shock  and  Vibration.  It  is  felt  that  the 
G-S  procedure  is  a  worthwhile  contribution 
which  will  be  invaluable  to  the  engineer  faced 
with  the  problem  of  designing  a  suitable  isola¬ 
tion  system  to  protect  delicate  airborne  equip¬ 
ment.  Obviously,  this  paper  has  presented  an 
excess  of  working  details  on  the  various  steps 
of  the  procedure.  This  intentional  overabun¬ 
dance  of  particulars  is  merited  since  experience 
with  the  G-5  procedure  has  shown  that  the  aver¬ 
age  engineer  encounters  difficulties  when  he 
first  attempts  to  apply  the  procedure.  The 
philosophy  has  been  followed  that  too  much  in¬ 
formation  is  better  than  insufficient  informa¬ 
tion.  As  this  paper  is  being  prepared  (18  July 
1961),  the  G-5  procedure  is  undergoing  revi¬ 
sions  which  are  not  reported  in  this  paper.  But 
the  essential  portions  of  the  procedure  have 
been  included  in  this  paper. 

It  is  predicted  that  mechanical  admittance 
measurements  will  become  commonplace  and 
that  the  application  of  mechanical  admittance 
concepts  will  prevail  within  the  next  decade. 
Future  dynamic  progress  will  depend  on  the  de¬ 
velopment  of  mechanical  admittance  measuring 
systems,  the  measurement  of  mechanical  ad¬ 
mittance  on  more  Aerospace  programs,  and  the 
application  of  these  data  to  optimize  structural 
and  equipment  designs. 
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NOMENCLATURE 

A  ■  Real  part  of  mechanical  admittance, 
ln./aec/lb 

A  =  Real  part  of  normalized  mechanical 
admittance  (Eq.  (3)),  dimensionless 

B  s  Imaginary  part  of  mechanical  admit¬ 
tance,  in. /sec/lb 

§  -  Imaginary  part  of  normalized  me¬ 
chanical  admittance  (Eq.  (4)),  dimen¬ 
sionless 

D  =  Function  defined  by  Eq.  (15),  dimen¬ 
sionless 

H  =  Transfer  function  defined  by  Eqs.  (5) 
and  (6),  dimensionless 

R  =  Real  part  of  mechanical  impedance, 
Ib/in./sec 

s  =  ^ectral  density  amplitude  of  random 
motion 

U  B  Complex  amplitude  of  sinusoidal 
motion 

v  =  u  =  velocity,  in. /sec 

w  =  Weight  of  equipment  s  mg,  pounds 

X  z  Imaginary  part  of  mechanical  imped¬ 
ance,  Ib/in./sec 

Y  B  Complex  mechanical  admittance, 
in. /sec/lb 

z  s  Complex  mechanical  impedance, 
Ib/in./sec 

a  B  Peak  acceleration,  g 

c  B  Coefficient  of  viscous  damping, 
lb/in./Bec 

f  B  Frequency  b  o/an ,  cps 

B  B  Acceleration  of  gravity,  386  in./sec^ 

i  =  j  » 

k  B  stiffness  of  linear  spring,  Ib/in. 
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tn  s  Mass  of  equipment  »  w/g ,  lb/in./sec^ 

u  =  Displacement,  inch 

/3  =  Fraction  of  critical  damping,  dimen¬ 
sionless 

u  =  Frequency  =  7vf ,  rad/sec 


SUBSCRIPTS 

o  -  Subscript  zero  Indicates  that  equip¬ 
ment  is  attached  to  source  without 
isolators 


[1]  SAE  G-5  Committee  on  Aerospace  Shock 
and  Vibration,  "Design  of  Vibration  Isola¬ 
tion  Systems." 

[2]  Dr.  Sheldon  Rubin,  "Theoretical  Back¬ 
ground  of  Procedure,"  Appendix  3.2  of 
Ref.  1. 


01.02  -  Used  with  f  to  indicate  frequencies  at 
which  first  peak  oi,  and  first  valley  02 
occur  on  the  curve 

f  -  Indicates  fragility  amplitude 

i  -  Indicates  that  quantities  refer  to 
isolators 

m  -  Used  with  H  to  denote  maximum  allow¬ 
able  transfer  function 


t  -  Indicates  quantities  of  source  structure. 


[3]  Dr.  Sheldon  Rubin,  "Design  Procedure  for 
Vibration  Isolation  on  Non-Rigid  Supporting 
Structures,"  SAE  National  Aeronautic  Meet¬ 
ing,  The  Ambassador  Hotel,  Los  Angeles, 
California,  October  5-9,  1959. 

[4]  H.  Himelblau  and  F.  Schloss,  "Measurement 
of  Mechanical  Admittance,"  Appendix  3.7 

of  Ref.  1. 


DISCUSSION 


Mr.  Stern  (General  Electric):  The  Isolators 
that  you  are  referring  to,  are  you  going  to  de¬ 
sign  and  make  these,  or  do  you  plan  on  buying 
commercial  isolators? 

Mr.  Mustain:  This  of  course  was  a  sam¬ 
ple  problem.  However,  we  have  ordered  from 
one  of  the  vendors,  under  a  little  pressure,  an 
isolator  unit  0. 18  damping  and  something  higher 
than  20  cycles  per  second. 

Mr.  Stern:  Generally  you  have  to  use  iso¬ 
lators  that  are  commercially  available,  and  the 
response  curves  are  usually  given  in  the  catalog 
literature  of  the  vendors.  This  really  brings 
up  one  of  the  objections  I  have  to  this.  Look  at 
the  isolators  that  are  used,  let's  say  if  you  were 
to  go  to  a  conventional  spring,  just  using  a  spring 
you  won't  get  any  damping.  If  you  use  a  mesh- 
type  isolator  to  get  damping  usually  you  are 
worried  about  icing  of  the  isolator  during  the 
humidity  test,  so  you  won't  isolate  that  way. 

The  re^Tular  rubber  isolators  get  hard  as  a 
rock  at  say  -65  degrees.  So  —  for  practically 
all  military  work  you  end  up  with  a  silicon  rub¬ 
ber  isolator  of  some  sort,  and  the  spring 


constants  are  far  from  linear.  If  anything  they 
are  pretty  close  to  a  tangent  function  or  a  hyper¬ 
bolic  sine  function.  So,  how  would  you  actually 
apply  this  technique  you've  evolved  to  an  isolator 
that  you  are  going  to  buy  from  a  catalog  when 
you  know  right  away  the  characteristics  are  far 
from  those  for  which  you've  set  up  in  this  plan? 

Mr,  Mustain:  Well  I  don't  think  you  can 
buy  it  off  the  shelf.  You  will  have  to  make  a 
special  specification,  and  this  is  the  second 
procedure  that  is  worked  out  in  this  G-5  Docu¬ 
ment;  How  to  specify  a  vibration  isolator. 

Mr.  Stern:  After  all  when  we  finally  ar¬ 
rived  at  the  silicon  rubber  isolator,  it  dMs 
work  fairly  well  compared  to  what  was  avail¬ 
able  before.  What  I  am  wondering  is  do  you 
really  think  people  are  going  to  pay  the  money 
and  take  the  time  to  design  isolators  when  we 
have  some  that  work  fairly  well  now  ? 

Mr.  Mustain:  Well,  as  I  said,  we're  forc¬ 
ing  the  manufacturer  to  give  us  0. 18  damping 
and  it  was  like  pulling  teeth  but  we  got  it  from 
him. 
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Mr.  Stern:  There  is  a  report  by  Mindlin, 
I'm  pretty  sure  you  are  familiar  with  it,  in 
which  he  has  handled  these  nonlinear  cases 
quite  well.  I  mean,  he  has  shown  how  it  is 
quite  easy  to  handle  them,  and  what  I  won¬ 
dered  was  why  you  haven't  made  any  attempt 
to  handle  these  cases  in  this  report.  Because 
these  are  the  isolators  people  use. 


Mr.  Mustain:  Well,  we  are  working  on  the 
problem.  Procedures  have  been  written  by 
Harry  Klein  of  Ramo-Wooldridge  that  will  cover 
the  specifications  to  vibration  manufacturers 
and  we  are  considering  all  these  points.  It  will 
probably  be  a  more  complex  specification  than 
any  that  has  ever  been  prepared  in  the  past. 


*  *  * 
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STRUCTURAL  RESPONSE  TO  DYNAMIC  LOAD* 


R.  M.  Mains 

General  Engineering  Laboratory 
General  Electric  Company 
Schenectady,  New  York 


This  paper  represents  an  attempt  to  generalize  the  problem  of  calcu¬ 
lating  responses  to  random  vibration  and  shock  into  a  set  of  simple 
principles,  which  are  sufficient  to  produce  numerical  solutions  to  prac¬ 
tical  problems.  Within  the  limitations  of  superposition  and  linearity, 
there  are  no  restrictions  on  the  recipes  given  so  that  the  methods  are 
general.  The  use  of  a  digital  computer  with  narrow  frequency  or  time 
intervals  leads  to  solutions  as  precise  as  desired,  yet  the  numerical 
summations  can  be  done  by  slide  rule  or  hand  computer  and  still  give 
acceptable  results.  The  related  subjects  of  load  definition  and  damage 
evaluation  are  discussed  so  that  the  prediction  of  structural  response 
is  placed  in  a  proper  frame  of  reference. 


UST  OF  SYMBOLS 
f(t)  3  a  time  function 

e  3  base  of  natural  logarithms 
t  3  time 
ui  3  frequency 
i  3  /T 

F(  iu)  3  the  complex  Fourier  integral  trans¬ 
form  of  f(t) 

F„  3  the  real  part  of  F(iu) 

F|„  3  the  imaginary  part  of  F(iw) 

M  3  mass  matrix 
C  3  damping  matrix 
K  3  stiffness  matrix 
X  3  motion  vector 
X,  X  3  time  derivatives  of  X 
P  3  operator,  d/dt 


Z(p)  3  Mp*  +  Cp  +  K 

X(  t )  3  time  function  of  X 

k(ia>)  3  the  Fourier  integral  transform  of  X(t) 

Y(  itu)  3  the  complex  frequency  response 
function 

a  3  the  real  (or  even,  or  cosine)  part  of 
the  response  vector 

b  3  the  imaginary  (or  odd,  or  sine)  part  of 
the  response  vector 

f  1  3  the  real  (or  even,  or  cosine)  part  of 
f{t) 

f]  3  the  imaginary  (or  odd,  or  sine)  part  of 
f(t) 

Yue  =  the  real  part  of  Y(ia>) 

Ym  3  the  imaginary  part  of  Yfiu) 

g(  iu)  3  complex  spectral  density  in  units  V 
cycle 

n  3  number  of  expected  cycles  at  a  given 
stress  level 


♦This  paper  was  not  presented  at  the  Symposium. 
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N  s  number  of  cycles  to  failure  at  a  given 
stress  level 


THE  PROBLEM  OF  DYNAMIC  DESIGN 

The  problem  of  dynamic  design  is  to  pro¬ 
duce  a  piece  of  equipment,  one  of  a  kind  or  in 
quantity,  which  performs  a  desired  set  of  func¬ 
tions  with  a  given  reliability  and  within  cost, 
size,  weight,  and  other  limitations.  Once  the 
functional  requirements  have  been  established 
and  the  economic  factors  have  received  due 
consideration,  the  problem  breaks  down  into 
three  essential  steps:  (1)  definition  of  the 
loads,  (2)  prediction  of  the  responses  of  the 
equipment  to  these  loads,  and  (3)  the  evalua¬ 
tion  of  the  resulting  performance.  In  the  gen¬ 
eral  sense,  loads  are  any  factors  which  by  their 
action  lead  to  malfunction  or  damage  of  the 
equipment,  while  responses  are  the  conse¬ 
quences  of  these  loads.  The  concern  of  this 
discussion  is  with  a  particular  variety  of 
loads  and  responses:  those  associated  with 
forces  or  motions  imposed  as  loads  upon  a 
mechanical  structure.  Such  specific  loads  are 
usually  either  static,  oscillatory,  or  transient 
(or  combinations  thereof)  and  are  referred  to 
as  steady  load,  vibration,  and  shock,  respec¬ 
tively.  For  any  problem  involving  the  produc¬ 
tion  and  use  of  several  pieces  of  the  same 
equipment  over  an  extended  period  of  time,  the 
loads  will  usually  vary  with  time,  with  the  par¬ 
ticular  (>iece  of  equipment,  and  with  the  indi¬ 
vidual  condition  of  use.  Consequently,  the 
loads  tend  to  be  random  in  character  and  to 
have  some  range  within  which  they  are  re¬ 
stricted  so  that  they  are  best  described  in 
terms  of  probabilities  and  statistics. 

The  response  of  a  linear  structure  to  a 
random  load  is  also  random,  but  of  a  modified 
character  [1].  The  way  in  which  the  response 
differs  in  character  of  randomness  from  the 
load  depends  upon  the  structure  itself  and  is 
predictable  for  linear  structures.  This  is 
covered  in  some  detail  below. 

The  development  of  malfunction  or  failure 
under  load  is  usually  subject  to  statistical 
variation,  even  for  steady  loads.  As  a  result, 
the  evaluation  of  performance  under  load  to 
determine  a  reliability  (or  its  inverse,  prob¬ 
ability  of  malfunction)  becomes  a  problem  of 
comparing  probable  malfunction  data  to  pro¬ 
duce  a  reliability  estimate.  This  amounts  to 
predicting  the  in-service  mortality  tables  in 
advance  of  actual  service. 


HOW  LOADS  ARE  DEFINED 

Since  we  are  concerned  with  force  or  mo¬ 
tion  loads  on  a  mechanical  structure,  the  am¬ 
plitude  variation  of  the  load  with  time  or  with 
frequency  must  be  described  in  order  to  define 
the  load.  Figure  1  shows  a  number  of  varia¬ 
tions  of  load,  together  with  the  names  by  which 
they  are  usually  called.  The  usual  processes 
for  determining  load  variations  are: 

(1)  By  direct  measurement  in  service  of  a 
prototype  or  model. 

(2)  By  deduction  from  in-service  measure¬ 
ment  of  similar  equipment. 

(3)  By  study  of  the  probable  environments, 
deduction  of  their  effect  upon  the  equipment,  and 
translation  of  this  information  into  loads. 

The  process  of  evolution  of  loads  informa¬ 
tion  is  usually  in  reverse  order  to  that  given 
above  because  the  prototype  is  never  available 
at  the  beginning  of  a  design.  Because  equipment 
reliability  is  a  function  of  the  entire  life  of  each 
individual  in  the  whole  population,  it  is  desirable 
to  describe  each  kind  of  load  in  terms  of: 

(4)  Root- mean- square  amplitude. 

(5)  Distribution  of  amplitudes  (histogram). 

(6)  Time  duration  or  total  design  life. 

(7)  Variations  in  (4)  and  (5)  with  time  or  at 
different  stages  of  design  life. 

For  all  loads  except  static  loads,  the  fre¬ 
quency  spectrum,  F(iu),  of  the  loads  is  needed 
rather  than  the  time  function,  f ( t ) ,  even  though 
load  measurements  are  usually  made  as  time 
functions.  To  change  from  the  time  function  to 
the  frequency  spectrum,  the  mathematical  proc¬ 
ess  used  is  ^e  determination  of  the  Fourier 
transform  (or  the  infinite  series  of  Fourier 
component  coefficients).  This  is  done  by  [2] 

• 

F{io;)  =  f  f{t)e’‘“*dt.  (1) 

Since  for  practical  problems,  f(t)  is  rarely  a 
convenient  analytic  function,  numerical  integra¬ 
tion  can  be  performed  by  substituting 

e'*“*  =  CO*  uit  -  i  sin  wt  ,  (2) 
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I  SINUSOIDAL 


FREQUENCY  ♦ 

Fig.  I  -  Different  kinds  of  load-time  functions 


whence 

•"re  =  ^  CO*  (3) 

•^iii  =  Z]  wt)  At  .  (4) 

t 

The  same  result  can  be  obtained  by  passing  the 
electrical  signal  representing  f(t)  through  a 
series  of  electrical  filters,  whose  characteris¬ 
tics  are  known,  and  deducing  from  the  filter  out¬ 
put  the  values  of  F(  iw) .  The  information  on 
phase  is  usually  lost  in  the  filtering  process, 
so  that  the  resulting  spectrum,  F(i<i>) ,  is  useful 
for  calculating  responses  to  random  oscillation, 
but  not  to  transients. 

It  sometimes  happens  that  the  histogram  of 
item  (5)  will  show  a  regular  distribution,  such 
as  Gaussian  or  Rayleigh,  and  the  variations  of 
item  v7)  will  not  occur.  In  such  a  case  the  func¬ 
tion  r(t)  may  be  what  is  called  a  "stationary 


time  function"  and  much  of  the  existing  mathe¬ 
matics  of  statistics  can  be  applied  [3],  In  other 
cases,  the  "stationary  time  function"  require¬ 
ments  may  be  fulfilled  for  short  intervals  in 
the  design  life,  or  may  be  assumed  to  be  ful¬ 
filled  so  that  advantage  can  be  taken  of  available 
statistical  theory.  This  is  primarily  helpful  in 
finding  representative  functions  and  testing  for 
correlations  and  fit  of  data. 

To  summarize  the  problem  of  definition  of 
loads,  consider  Fig.  2.  We  start  with  the  meas¬ 
ured  time  function  of  the  load,  fft),  filter  it  or 
transform  it  to  get  the  frequency  function  of  the 
load  F(iw).  The  frequency  functions  are  studied 
statistically  and  the  results  weighted  by  judgment 
to  account  for  the  degree  of  uncertainty,  the 
severity  of  effect  of  malfunction,  and  the  like. 

The  result  is  a  set  of  design  values  of  load, 
statistically  described,  the  general  character  of 
which  falls  in  one  or  more  of  the  categories  of 
Fig.  1. 
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Fig.  2  -  Flow  chart  for  load  definition 


HOW  RESPONSES  ARE  PREDICTED 


Substitute  Eqs.  (9)  and  (10)  into  (8)  and  get 


First,  consider  the  equations  of  motion  of 
a  linear  system  with  lumped  parameters, 

MX  +  cx  +  KX  =  f(t)  ,  (5) 

in  which 

M  s  mass  matrix 
c  =  damping  matrix 
K  =  stiffness  matrix 
X  =  motion  vector 
X  and  X  =  time  derivatives  of  X 

f ( t )  =  time  function  of  load,  a  vector. 

Now  let 


and  let 

Z(p)  =  Mp2  +  Cp  +  K  :  (7) 

then,  (5)  becomes 

Z(p)  X(t)  =  f(t)  ,  (8) 

in  which  X(t)  is  used  for  X  simply  to  empha¬ 
size  that  it  is  a  time  function.  By  Fourier's 
theorem  and  Eq.  (1)  [2] 

W 

f(*)  =  ^  f  F(ic.j)e‘“'da..  (9) 

•  • 

Define  a  displacement  vector  function  of  fre¬ 
quency,  x(  i-')  such  that 

^  f  e‘“'  .  (10) 

•i  m 


I  ‘  dciJ  =  ^  f  F(if^)  e“"‘dai.  (H) 

*•  ®  "•tZ 

Carry  Z(p)  inside  the  integral  and  get 

m 

^  f  Z(i<ij)  e““*  V(iai)  do) 

as 


■  ^  f  Ffioj)  e*"*  do)  .  (12) 

OD 

Equation  (12)  then  reduces  to 

Z(iaj)  y(iu,)  =  F(i^)  .  (13) 

Pre  multiply  by 

Z’  Viol) 

and  define 

Y(iw)  =  Z'^Cia)  ,  (14) 

whence 

Vfio;)  =  Y(iw)  Ffiw)  .  (15) 

In  Eq.  (15),  Yffiu)  is  the  calculated  or  meas¬ 
ured  mobility  matrix.  F(iM)  is  obtained  by 
transform  from  Eq.  (1);  then  v(i«)  can  be 
calculated.  Equation  (10)  then  serves  to  trans¬ 
form  y(i<*j)  into  X(t),  if  the  time  function  is 
needed. 

In  the  application  of  these  equations  to 
practical  problems,  two  separate  cases  occur; 

(a)  Preliminary  design,  for  which  calcu¬ 
lated  or  estimated  values  are  put  in  Eq.  (15) 
and  carried  through  to  a  solution. 

(b)  Final  design,  or  design  refinement,  for 
which  measured  data  can  go  into  Eq.  (15). 
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At  intermediate  stages  between  (a)  and  (b),  the 
input  for  Eq.  (15)  may  be  partly  calculated  and 
partly  measured. 

For  the  numerical  calculation  of  ¥(  iu) ,  it 
is  convenient  to  consider  Eq.  (5)  and  let 


X  =  a  cos  oit  +  b  sin  wt 
f(t)  =  fj  (cos  ut)  +  fj(ain  at)  . 


(16) 


When  these  substitutions  are  made  in  Eq.  (5), 
two  equations  result,  and  the  order  of  the  matrix 
solution  is  doubled, 

-  +  (iiCb  +  Ka  =  fj 

-  -  aCa  +  Kb  =  f,  , 


which  can  be  written  as 
■  M  +  K  +  ac 
-ac  -a^M  +  K 


a 

fi 

b 

^2 

(17) 


Equation  (17)  now  permits  the  computation  of 
the  vector 


for  any  desired  set  of  values  of  a,  with  the 
force  vector 


gives  the  first  row  and  column  of 


•re 


(19) 


and 


'lU 


(20) 


Similarly,  with  the  j  th  row  of  the  f  vector  of 
Eq.  (18)  taken  as  unity,  the  j  th  row  and  column 
of  Eqs.  (19)  and  (20)  result.  Even  for  as  few  as 
3  degrees  of  freedom,  the  bookkeeping  involved 
in  these  computations  is  staggering,  so  that  ma¬ 
chine  computation  is  advisable.  The  program¬ 
ming  for  the  machine,  however,  is  relatively 
simple,  since  only  a  simultaneous  equation  rou¬ 
tine  and  bookkeeping  are  involved. 

To  complete  the  set  of  relations  for  numeri¬ 
cal  computation,  recall  that 


=  cos  ojt  +  i  sin  at 


(21) 


taken  such  that  all  the  fj’s  are  zero,  and 
each  f  J  is  successively  unity  while  the  other 
f|’s  are  zero.  The  resulting  n's  give  the 
real,  or  in-phase,  components  of  Y(  ia) ,  and 
the  b's  give  the  imaginary,  or  out-of-phase, 
components  of  Y(  iw) .  For  example. 


-  uJC 


1 

"2 

0 

"3 

0 

"n 

_ 

bl 

ba 

0 

(18) 


and  apply  this  to  Eq.  (10)  to  get 

*  ^  [xrE**""*  +  Xii,  cosut]  Ao).  (22) 


The  foregoing  discussion  is  particularly 
interesting  because  it  leads  to  the  conclusion 
that  any  driving  function  which  can  be  expressed 
as  a  function  of  frequency,  Eqs.  (3)  and  (4),  to¬ 
gether  with  the  system  mobility,  Eqs.  (19)  and 
(20),  can  be  put  into  Eq^  (15)  to  give  the  response 
of  the  system  as  a  function  of  frequency.  If  the 
response  is  desired  as  a  function  of  time,  then 
Eq.  (22)  provides  the  means  for  doing  this  also. 
Of  course,  the  units  in  the  various  terms  of  the 
equations  must  be  compatible,  and  the  conditions 
of  linearity  and  superposability  must  apply. 

When  the  driving  functions  are  statistically  de¬ 
scribed,  then  the  various  equations  apply  to  the 
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mean- square  or  root-mean- square  values,  and 
Eq.  (15)  is  squared  to  accommodate  this.  For 
example, 

[x(ic.;)]^  =  X  Ga<^)  . 

in  which 

\  x(  iu)]  ^  =  mean  square  response  density, 
(units)  ^/cycle, 

[Y(iai)]*  =  square  of  |Y(iai)  1  ,  in  units 
suitable  to  change  G  into  x, 

G(  iu)  =  spectral  density  in  (units)  */ 
cycle, 

and  the  difference  between  Eqs.  (15a)  and  (15) 
is  essentially  only  that  between  the  customary 
description  of  a  continuous  random  spectral 
density  as  mean-square  units/cycle,  and  a  dis¬ 
crete  spectrum  as  just  units. 

The  foregoing  is  summarized  in  a  group  of 
charts  for  ready  reference.  Figs.  3,  4,  and  5. 
Figure  3  summarizes  and  illustrates  the  prin¬ 
ciples  involved.  Figure  4  shows  some  time- 
function  ■<  frequency-function  transforms  for 
some  common  driving  functions.  Figure  5  il¬ 
lustrates  what  can  be  done  to  solve  specific 
problems  numerically.  These  three  charts  pro¬ 
vide  the  essential  recipes  for  solving  numeri¬ 
cally  any  problem  to  which  the  limitations  of 
linearity  and  superposability  apply.  For  a  few 
(2  or  3)  degrees  of  freedom,  it  is  feasible  to  do 
the  work  by  hand  computation,  even  though  the 
bookkeeping  involved  is  tedious.  The  program¬ 
ming  for  machine  computation  is  simple,  since 
only  a  few  standard  subroutines  and  bookkeeping 
instructions  are  needed.  The  result  is  a  power¬ 
ful  tool  which  handles  all  structural  dynamics 
problems,  regardless  of  the  particular  name 
given  to  the  driving  function,  with  one  response 
equation  serving  all  cases,  Eq.  (15). 


HOW  PROBABLE  DAMAGE  IS  EVALUATED 

This  part  of  the  process  is  appreciably  less 
determinant  than  the  definition  of  loads  and  the 
prediction  of  responses  to  loads.  Investigations 
on  damage  accumulation  are  continuing,  and 
Refs.  3  and  4  contain  bibliographies  on  the 
subject.  Pages  365,  366,  and  389  in  Ref.  1  give 
numerical  examples  of  damage  accumulation 
calculations,  so  such  calculations  will  not  be 
repeated  here. 


The  principles  involved  in  damage  evalua¬ 
tion  are  relatively  simple  and  stralghtfoward, 
but  their  application  requires  careful  book¬ 
keeping.  First,  one  must  have  a  damage  factor, 
such  as  stress,  defined  statistically  so  that  a 
histogram  can  be  drawn  to  relate  stress  ampli¬ 
tude  to  the  number  of  occurrences.  Then  one 
must  have  data-relating  stress  amplitude  to  the 
number  of  cycles  required  to  produce  failure 
(or  malfunction).  Finally,  one  must  either  have 
data  or  make  some  assumption  about  how  the 
damage  accumulates  from  each  stress  cycle 
until  failure  occurs. 

For  example,  with  a  statistically  defined 
set  of  loads,  the  methods  discussed  above  can 
be  used  to  predict  responses  and  generate  a 
histogram  of  the  expected  stress  amplitudes. 
Next,  S-N  data  for  the  material  can  either  be 
found  in  the  literature  or  generated  by  test.  It 
is  then  common  to  assume  a  damage  accumu¬ 
lation  rule  such  as  the  linear  one  [5],  which 
states  that  a  proportion  of  total  damage  accrues 
at  each  stress  level,  ,  which  is  the  ratio  of 
the  number  of  expected  cycles  at  that  stress 
level,  n ,  to  the  number  of  cycles  required  to 
produce  failure  at  that  stress  level  in  the  virgin 
material,  N .  Tests  to  date  show  that  when 
these  portions  of  total  damage  add  to  something 
between  0.3  and  3.0,  failure  occurs.  In  symbols, 

^  =  0.3  to  3.0.  (23) 

This  is  not  a  very  satisfactory  situation 
for  stress  damage,  and  we  know  more  about 
stress  damage  accumulation  than  any  other 
damage  producing  factor.  Much  more  work 
needs  to  be  done  in  this  area. 


THE  EFFECT  OF  DAMPING  ON 
RANDOM  VIBRATION  AND  SHOCK 

First,  consider  the  two  cases  of  a  constant 
spectral  density  random  vibration,  and  an  im¬ 
pulse  shock,  for  both  of  which  the  frequency 
function  of  load  is  a  straight,  horizontal  line  as 
in  Fig.  4.  The  response  of  a  structure  to  these 
loads  is  then  governed  by  the  appropriate  fre¬ 
quency  response  function,  as  in  Fig.  3.  The  ef¬ 
fect  of  damping  on  the  response  can  then  be 
studied  by  considering  how  the  frequency  re¬ 
sponse  function  is  affected  by  damping.  It  is 
clear  that  a  decrease  in  any  part  of  the  fre¬ 
quency  response  function  decreases  the  re¬ 
sponse,  and  a  decrease  in  area  under  the 
frequency-response-function  curve  correspond¬ 
ingly  decreases  the  response. 
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Fig.  3  -  Summary  of  principles  for  response  calculation 
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Fig.  4(a)  -  Sample  time-function  to  frequency-f\inction  transforms 


FQW  SPECiriC  PWOeLEMS  -  (UNITS  OF  FUNCTIONS  MUST  BE  COMPSTIBLE I 


EO(IS)  CAN  BE  SOLVED  FREOUENCV- BY  •  FREOUENCT  TO  GET  X(<ii) 


AND  SBNLARLV  FOR  F,»|  t  ♦  FrC  »  V|M  ■  XlM 
EQ  (I ) ,  the  MRCCT  rauRCR  INTCSRAL  TRANSFORM  CAN  BE  SOLVED  NUMERtCAUV  BY 
Sc  (<•>)*?*<•)  COSalAl  ,  FM|(m)«  ^  I  (I  )(-SiN  •>!  )At 

FOR  EACH  SFSOFIC  FREOUENCV 

EO  (10) ,  THE  MVER8E  FOURIER  INTEGRAL  TRANSFORM  CAN  BE  SOLVED  NUMERICALLY  BY 

I  (I)  •  ^  3  [*RC  COS  -  X,M  ****  «•  ]  Am  T 

i  *(*Rt  SIN  .t  ♦  COS-l)  A. 

FOR  EACH  SPECIFIC  TIME 


FOR  RANDOM  ORMNB  FUNCTIONS  EXPRESSED  AS  A  MEAN  SQUARE  SPECTRAL  DENSTTV, 
G(l«).  USE  [Y(i«)]*  'Vm  Ym  SMC  MULTIPLY  FREQUENCY  •  BY  -  FREQUENCY  TO 
GET  THE  MEAN  SQUARE  RESPONSE  SPECTRAL  DENSITY  ,  [l(l»)]*: 


Fig,  5  -  Summary  for  numerical  calculation  of  response 


In  terms  of  the  (requency-by-frequency 
multiplication  of  load  function  by  frequency  re¬ 
sponse  function  to  get  the  response  function 
(Fig.  5),  the  response  is  reduced  by: 

1.  A  general  decrease  in  load  function 
level. 

2.  A  general  decrease  in  frequency- 
response-function  magnitudes. 

3.  A  shift  of  location  of  load  or  frequency- 
response-function  peaks  so  that  peaks  in  one 
curve  do  not  coincide  with  peaks  in  the  other 
curve. 

Now  a  general  increase  of  damping  can  ac¬ 
complish  item  2,  so  that  a  general  decrease  of 


response  is  obtained.  Discrete  increases  of 
damping  can  reduce  the  response  in  particular 
frequency  ranges  by  reducing  particular  peaks 
of  the  transfer  function.  Shifting  the  location 
of  peaks  in  the  frequency  reqmnse  function  is 
best  done  by  changes  in  stiffnesses  or  mass 
loadings. 

Consequently,  the  effect  of  damping  on  re¬ 
sponse  to  random  vibration  and  shock  is  to 
determine  the  height  of  peaks  in  the  frequency 
response  function  and  its  general  level.  In 
general,  an  increase  of  damping  reduces  re¬ 
sponse  and  a  decrease  of  damping  increases 
reqranse,  though  there  may  well  be  specific 
cases  for  which  this  generalization  does  not 
apply.  Reference  1,  Chapter  5,  summarizes 
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the  state  of  formal  knowledge  concerning 
structural  damping. 

For  any  specific  problem,  the  effect  of 
experimental  changes  in  damping  on  a  model 
or  prototype  transfer  function  can  be  meas¬ 
ured,  and  the  analytical  effects  of  changes  in 
damping  can  be  determined  from  Eq.  (17).  The 
qualitative  effects  are  as  just  stated. 


SUMMARY 

One  might  well  ask  what  all  of  the  foregoing 
means  in  terms  of  the  state  of  the  art  in  random 
vibration  and  shock.  What  has  been  discussed 
here,  together  with  the  references  in  the  bibli¬ 
ography,  leads  us  to  state  that: 

1.  We  can  with  certainty: 

a.  Measure  loads  and  describe  them 
statistically. 

b.  Measure  (on  prototypes  or  models) 
the  appropriate  transfer  functions. 

c.  Measure  (or  predict)  the  probable 
responses  to  the  loads  in  terms  of  known  dam¬ 
age  producing  factors. 


2.  We  can,  but  with  uncertainty: 

a.  Estimate  loads  and  devise  statisti¬ 
cal  descriptions  of  them  when  measurements 
are  not  feasible. 

b.  Devise  mathematical  models  to  cal¬ 
culate  transfer  functions  that  are  needed. 

c.  Hypothesize  damage  criteria  and  ap¬ 
ply  them  to  the  measured  or  calculated  re¬ 
sponses. 

d.  Judge  the  validity  of  the  process  and 
data  we  have  used,  and  form  an  opinion  of  the 
resulting  probability  of  damage. 

3.  We  need  more  information  on: 

a.  Damage  criteria  and  damage  accu¬ 
mulation. 

b.  The  validity  of  various  procedures 
for  devising  mathematical  models  to  represent 
physical  problems. 

c.  Other  damage  processes  than  stress 
fatigue,  which  may  be  susceptible  to  similar 
treatment. 


APPENDDC 
Worked  Examples 


Two  examples  have  been  worked,  one  to 
show  the  mathematics,  and  one  to  show  a 
numerical  calculation. 


EXAMPLE  1 

Consider  the  single-degree-of-freedom 
system  shown  in  Fig.  A-1.  The  equation  of 
nx>tion  is: 


mx  +  c*  +  kx  =  cx„  +  kx„. 

o  o 


Let 


X 


X 


o 


fl  cos  oit  >  b  sin  a)t 
A  cos  at  ^  B  sin  at 


}■ 


(A-l) 


(A-2) 


Substitute  (A-2)  in  (A-l),  and  ge't 

(-i<.’^nia  +  ka  +  a)cb)  cos  uit  ♦  (-oj^  nt>t  kb  -  ojcii)  sin  <4)t  = 

♦  (.*.’cB  +  kA)  cos  at  +  (-a'cA  +  kB)  sina)t.  (A-3) 
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Fig.  A-l  -  Single-degree-of- 
freedom  system 


(A-3)  can  be  written; 


whence 


-  ei)*m  4  k  4  a)c  ^ 

a 

oicB  4  kA 

Xbe  ■  ^be  7**  “  Fm  Yjk  I 

= 

.  (A-4) 

-  aie  -  a)*m4k 

b 

-  aieA*  kB 

Xju  =  FjjYjy  4  Fi^Ygj  J 

whence  and 


(A- 12) 


(-<ii*ni  +  k)k  +  (fttc)*  ^ 

(-w*m  +  k)*  +  (uc)^ 


>^‘>  =  57  f' 

im 


*  ^IM  oit)  dct> 


a 


(Vm  •incut  +  VjyCO»ejt)dcj.  (A- 13) 


(- cj^m  ■»  k)(- cjc)  ■»  <»)ck 
(-  Cl)*  in+  k)*  +  (cjc)* 


(-  cj*iii  -f  k)k  4-  (cjc)*  g 
(-  ci)*in  +  k)*  +  (ci)c)* 


(A-6) 


If  we  note  that  ^aa  is  an  even  function  while 
Fjm  is  odd,  that  is  even  while  Yj^  is  odd, 
it  follows  that  ;tK£  is  even  and  V|y  is  odd.  As 
a  consequence,  the  imaginary  part  of  (A- 13) 
vanishes  and  we  have 


An  examination  of  Eqs.  (A-5)  and  (A-6)  shows 
that  the  direct  frequency  response  function  is 


X(t) 


(-  cj*ni'fk)k  (CJC)* 

®  (-  cij*B  +k)*  +  (cue)*  (A-7) 


and  the  cross-frequency  response  function  is 


(- a)*iB  t  k)  ci)C  -  ci)ck 
(-  &)*!«  + k)*  +  (Ci)C)* 


(A-8) 


+  +  P||| Ypy)  sin  ci)t^  dci)  (A— 14) 

which  is  the  general  equation  for  the  time  func¬ 
tion  of  response  of  a  system  to  any  excitation 
which  has  a  proper  Fourier  integral  transform. 


EXAMPLE  2 


If  we  had  made  the  substitution 


at  Eq.  (A-2),  the  same  frequency  reqranse 
functions  would  result,  except  that  Y(.  would  be 
imaginary.  In  fact. 


Ybe  =  Yo  and  Y,£  =  Yc  •  (A-9) 


Now  with 


Consider  the  four-degree-of-freedom  sys¬ 
tem  shown  in  Fig.  A-2.  We  are  concerned  with 
the  stretch  of  the  q>ring,  kjj ,  which  connects 
masses  2  and  3,  under  a  step  shock  input  and 
under  a  random  vibration.  The  shock  load  is  a 
100  in. /sec  velocity  step  applied  to  the  founda¬ 
tion,  and  the  random  vibration  is  described  in 
Figs.  A-3  and  A-4. 

The  structure  is  characterized  by  the  mass 
matrix. 


r(iai)  =  I  dt  ,  (A-10) 


we  get 


■^XE 


=  f.« 


t)  COl  alt  dt 


.  (A-11) 


=  I  f(t)(-  iin  ait)  dt 

J 


M  » 


7.13 

0 

0 

0 


0 

4.38 

0 

0 


the  stiffness  matrix 


if40.0 

4^  6.62 

-  2.85 
4  5.54 


46.62 

47.59 

-0.445 

41.107 


0 

0 

1.13 

0 


0 

0 

0 

1.89 


-2.85 

-0.445 

41.057 

-0.395 


45.54 

41.107 

-0.395 

49.00 


.  (A- 15) 


xl0»,  (A- 16) 
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40 


*0 


y  y  ^  y  y  X  ^  ^  V  X  v 


joRIVINO  FUNCTION 


FOUNDATION 


notes:  each  sfnino  is  paralleled  by  a  viscous  damper 

MOTIONS  ARE  CONSTRAINED  TO  X  -  DIRECTION  ONLY 


Fig.  A-2  -  Four-degree-of-freedom  s/stem 
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FREQU£NCY-CPS 


Fig.  A-3  -  Random  vibration  spectral  density 


and  the  damping  matrix 

+267  +78  -50  +64 

+  78  +90  -12  +16.6 

C  * 

-  50  -12  +17.5  -  9.2 

+  64  +16.6  -  9.2  +65.4 

The  equations  of  motion  are  then 

Mx  +  Oi  +  Kx  =  f(t)  ,  (A- 18) 

and  we  substitute 


(A- 17) 

and 


fi 


f  2 


+49.31 
+14.872 
-  2.633 
+  15.252 


+359 
+  172.6 
-  53.7 
+136.8 


(A-21) 


X  =  a  cos  OiJt  +  b  sin  wt  (A- 19) 

to  get 

+  K  +  iiC  a  fj 

(A-20) 

-ciC  -CJ^M  +  K  b  fj 

Now  to  get  the  solution  for  driving  from  the 
foundation,  the  elements  of  the  f  j  vector  must 
be  the  sums  of  the  rows  of  the  K  matrix  and 
the  f  2  vector  the  sums  of  the  rows  of  the  C 
matrix  times  (-cj): 


This  produces  unit  response  at  zero  frequency. 
Equation  (A-21)  is  put  into  (A-20),  and  the  equa¬ 
tions  are  solved  simultaneously  for  the  a  (real) 
and  b  (imaginary)  vectors  for  a  range  of  values 
of  oj.  The  results  are  shown  in  Table  1  and 
Fig.  A- 5. 

To  compute  the  random  vibration  response 
frequency  spectrum,  we  use 

d(i&j)  *  i<ii)j  =  ^>f(i<i))j  ,  (A- 22) 

which  is  shown  in  Table  2  and  Fig.  A-6.  Since 
the  principal  response  is  in  the  frequency 
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Fig.  A- 4  -  Random  vibration  probability  plot 


TABLE  A<1 

Frequency  Reqranse  Functions 
(Sketch  in  spring  kjj  per  unit  foundation  motion) 


Freq. 

(cps) 

Ym 

(Y(i<u))® 

Freq. 

(cps) 

Yre 

Yim 

[Yfio.))* 

0 

0 

0 

0 

97.1 

>1.999 

-  2.150 

8.619 

\/2v 

-  1.43  X  10  ® 

+  1.326  X  10  * 

2x  10- 

100 

>0.911 

-  0.6955 

1.314 

10 

-  0.06001 

+  3.783  X  10-* 

0.00360 

no 

>2.124 

-  0.2039 

4.553 

20 

+  5.093  X  10  ® 

0.09048 

120 

>3.512 

-  0.6815 

12.80 

30 

0.04930 

1.236 

130 

-  4.924 

63.96 

40 

-  8.298 

+  2.100 

73.27 

134.8 

>1.489 

-10.804 

118.9 

42.9 

+  0.9642 

>40.09 

1608.1 

140 

-  5.831 

53.4 

45 

4 14. 246 

>  4.955 

227.5 

150 

-2.860 

-  1.494 

10.41 

50 

>  6.587 

>  0.4930 

43.63 

55 

^  6.504 

-  0.4172 

42.48 

180 

, 

7.69 

59.7 

^  3.206 

-13.22 

185.0 

250 

5.00 

-  0.5472 

0.4930 

400 

2.05 

^  1.119 

-  0.00580 

1.252 

700 

0.58 

^  1.612 

>  0.01065 

2.599 

lOOO 

0 

♦  2.753 

-  0.7177 

8.094 
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Fig.  A-5  -  Frequency  response  function 


TABLE  A-2 


Random  Vibration  Req>onse 


Freq. 

(cps) 

g(  ioj) 

gVcps 

1 

d(  iai) 

in.  */cps 

X  10-« 

V(i<u)* 

X(i^-)2 
in.* /cps 

X  10* 

Af 

(cps) 

X(ia))Af 
in.*  X  10*® 

rms 

in.  X  10'* 

0 

0.010 

0 

0 

0 

0.785 

0 

0 

l/2n 

0.013 

0.800 

2  X  10  - “> 

0 

5 

0 

0 

10 

0.100 

0.960 

0.00360 

3 

9.22 

28 

5.3 

20 

0.076 

45.7 

0.09048 

4 

10 

40 

6.3 

30 

0.088 

10.4 

1.236 

13 

10 

130 

11.4 

40 

0.142 

5.32 

73.27 

389 

6.45 

2510 

50.1 

42.9 

0.170 

4.80 

1608.1 

7720 

2.50 

19300 

139. 

45 

0.187 

4.38 

227.5 

998 

3.55 

3540 

59.5 

50 

0.214 

3.29 

43.63 

144 

5 

720 

26.6 

55 

0.227 

2.37 

42.48 

101 

4.85 

490 

22.1 

59.7 

0.188 

1.41 

185.0 

261 

5 

1305 

36.1 

65 

0.163 

0.875 

0.4930 

0 

7.65 

3 

1.7 

75 

0.132 

0.400 

1.252 

1 

10 

10 

3.1 

85 

0.108 

0.200 

2.599 

1 

10 

10 

3.1 

95 

0.098 

0.116 

8.094 

1 

6.05 

6 

2.5 

97.1 

0.097 

0.105 

8.619 

1 

2.50 

3 

1.7 

100 

0.097 

0.093 

1.314 

0 

6.45 

1 

1.0 

no 

0.102 

0.067 

4.553 

0 

10 

3 

1.7 

120 

0.112 

0.052 

12.80 

1 

10 

10 

3.1 

130 

0.124 

0.041 

63.96 

3 

7.40 

22 

4.7 

134.8 

0.130 

0.038 

118.9 

4 

5 

20 

4.5 

140 

0.136 

0.034 

53.4 

2 

7.60 

15 

3.9 

150 

0.145 

0.028 

10.41 

0 

20 

6 

2.4 

180 

0.172 

0.016 

7.69 

0 

50 

6 

2.4 

250 

0.135 

0.0033 

5.00 

0 

no 

2 

1.4 

400 

0.083 

0.0003 

2.05 

0 

225 

0 

0 

700 

0.030 

-- 

0.58 

0 

300 

0 

0 

1000 

0 

0 

0 

150 

0 

X  28180 

0 

X168.0 
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Fig.  A-6  -  Random  vibration  response 


range  around  42.9  cps,  a  fatigue  damage  evalua' 
tion  could  be  done  by  assuming  that  the  sum  of 
the  responses  (0.0281  in.^  mean-square  re¬ 
sponse,  0.168  in  root- mean- square  response) 
would  be  predominantly  42.9  cps.  The  damage 
evaluation  would  then  follow  the  process  of 
Ref.  1,  p.  389. 

To  compute  the  step  velocity  shock  re¬ 
sponse,  we  first  transform  a  unit  step  velocity 
Into  and  F,y;  thus 


F(M  =  l(t)e‘*’*dt.  (A-23) 


whence 


(A-24) 


We  then  get 

^  .  (A 

’‘’ill  =  ^IM^RE  J 

as  in  Table  3,  and  then  perform  the  inverse 
transform  to  get  x(  t ) , 


♦  ^im^be  ^  •  (A-26) 

The  resulting  time  function,  shown  in  Fig.  A-7, 
can  be  used  in  a  fatigue-damage  evaluation 
again  as  before. 


82 


fit)- 


TABLE  A-3 
Step  Shock  Response 


Freq. 

(cps) 

Fiu 

Vim 

’^FE 

Vjy 

0 

.  00 

0 

0 

0 

0 

1/27T 

-1 

- 

1.43  X  10- » 

+  1.33  X  10-» 

+1.33  X  10-9 

+1.43  X  10-3 

10 

-0.01592 

0.06001 

+  3.783  X  10-« 

+6.03  X  10-« 

+0.955  X  10-3 

20 

-0.00795 

- 

0.3008 

+  5.093  X  10-® 

+4.70  X  10- * 

+2.39  X  10-3 

30 

-0.00531 

- 

1.1105 

+  0.0493 

+2.62  X  10-9 

+5.90  X  10-3 

40 

-0.00398 

- 

8.298 

+  2.100 

+0.00837 

+0.0330 

42.9 

-0.00371 

+ 

0.9642 

+40.09 

+0.1486 

-0.00358 

45 

-0.00354 

+14.246 

+  4.955 

+0.1757 

-0.0504 

50 

-0.00318 

+ 

6.587 

+  0.4930 

+0.00157 

-0.0209 

55 

-0.00289 

+ 

6.504 

-  0.4172 

-0.00121 

-0.0188 

59.7 

-0.00267 

+ 

3.206 

-13.22 

-0.0353 

-0.00856 

65 

-0.00245 

- 

0.440 

-  0.5472 

-0.00134 

+0.001078 

75 

-0.00213 

+ 

1.119 

-  0.0058 

-1.24  X  10-s 

-0.002385 

85 

-0.00187 

+ 

1.612 

+  0.01065 

+1.99  X  10-s 

-0.003015 

95 

-0.00168 

+ 

2.753 

+  0.7177 

-1.205  X  10-3 

-0.00463 

97.1 

-0.00164 

+ 

1.999 

+  2.150 

-3.53  X  10-3 

-0.00398 

100 

-0.00159 

+ 

0.911 

+  0.6955 

-1.11  X  10-3 

-0.001448 

no 

-0.00145 

+ 

2.124 

-  0.2039 

-2.96  X  10-9 

-0.003224 

120 

-0.00133 

+ 

3.512 

-  0.6815 

-0.907  X  10-3 

-0.00467 

130 

-0.00122 

+ 

6.302 

-  4.924 

-6.00  X  10-3 

-0.00769 

134.8 

-0.00118 

+ 

1.489 

-10.80 

-1.275  X  10-3 

-0.001756 

140 

-0.00114 

. 

4.405 

-  5.831 

-6.65  X  10-3 

+0.00502 

150 

-0.00106 

- 

2.860 

-  1.494 

-1.58  X  10-3 

+0.00303 

180 

-0.00089 

2.50 

-  1.20 

-1.07  X  10-3 

+0.002225 

250 

-0.00064 

- 

2.00 

-  1.00 

-0.64  X  10-3 

+0.00128 

400 

-0.00040 

- 

1.30 

-  0.60 

-0.24  X  10-3 

+0.00052 

700 

-0.00023 

- 

0.40 

-  0.30 

-0.69  X  10-9 

+0.92  X  10-9 

1000 

-0.00016 

0 

0 

0 

0 
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